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An experiment has been completed 
which demonstrated quantum mechani- 
cal tunneling of electrons between two 
gold electrodes separated in vacuum. 
The tunneling current between the gold 
electrodes has been measured, for fixed 
voltages of 0.1 and 0.01 volts, as the 
electrode spacing was varied from a dis- 
tance of approximately 2.0 nm down to 
a point where the electrodes touched. 
Current changes of over five orders of 
magnitude were found for electrode 
spacing changes of approximately 1 .2 
nm. For the first time, these data enable 
one to deduce the work function of the 
electrodes in a tunneling experiment 
from experimental parameters indepen- 
dent of the tunneling device. Also ob- 
tained were current-voltage charac- 
teristics for fixed electrode spacings in 
the direct tunneling region where elec- 
trode spacings were less than 2.0 nm. 
An analysis is given which attempts to 
deduce an absolute electrode spacing 
and tunneling area from the nonlinear 
properties of the /-Fdata and the cur- 



rent versus spacing data. The analysis 
suggests that van der Waals and electro- 
static forces play a major role in deter- 
mining the I-V characteristics and that 
the tunneling area may be as small as 
10~ 16 m 2 . Along with a review of the 
theory of work functions and quantum 
mechanical tunneling, numerical calcula- 
tions of the tunneling current based on 
the free-electron model of the electrodes 
and the barrier, an image-potential re- 
duced barrier, and a WKB approxima- 
tion for the tunneling probability have 
been performed and compared with 
Simmons' theory and with the experi- 
mental results. 
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Preface 



The following paper is a reproduction of the au- 
thor's dissertation which was completed and sub- 
mitted to North Texas State University in August 
1978. A majority of the experimental work was 
performed at the National Bureau of Standards. 
Results from the experiments were first reported at 
the Spring 1978 APS meeting and later at a Na- 
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tional Science Foundation Workshop on Opportu- 
nities for Microstructures Science, Engineering, 
and Technology held in Airlie, VA, November 19- 
22, 1978. While a complete description of the ex- 
periments and results is available as a publication 
from University Microfilms, publication in the 
more accessible form represented here has been 
prompted by requests from colleagues and by the 
significant number of recent references to the 
work. In accordance with the requests, the Board 
of Editors of the Journal of Research has agreed to 
publish the document unchanged except for the ad- 
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dition of footnotes to identify recent work to 
which it is closely related. 

Renewed interest in these experiments has pri- 
marily resulted from the success and current wide- 
spread use of the scanning tunneling microscope 
(STM) introduced by Binning and Rohrer. 1 The 
findings of this work on vacuum tunneling should 
be of interest to many who are using or planning 
use of the STM. The current vs. voltage data, cur- 
rent vs. electrode spacing data, needs for tempera- 
ture stability, and observations about the potential 
influences of van der Waals forces in these static 
non-servoed tunneling experiments are important 
in understanding the servoed operation used for 
tunneling in the STM. While servoed operation 
overcomes much of the sensitivity of the STM to 
effects from vibration, thermal expansions, and ma- 
terials instabilities in the electrodes, these effects 
are proving to be significant factors in determining 
the stability and reproducibility of STM operation. 



1. Introduction 

According to quantum mechanics, a material 
particle can pass through a spatial region in which 
the particle's kinetic energy is less than its potential 
energy. The phenomenon is known as tunneling. A 
macroscopic example would be a ball which, when 
rolled toward a sharply peaked hill with insuffi- 
cient speed to get over the top, does not roll back, 
but appears on the other side of the hill without 
affecting the ball or the hill and continues on its 
way. The probability of such an event occurring is 
astronomically small at the macroscopic level. 
However, for electrons and other fundamental par- 
ticles and walls formed by nuclear and atomic 
forces, the probability while still small (one in 10 9 is 
not unusual) leads to major physical effects [l] 2 . 
This work pertains to the tunneling of electrons 
between two metal plates closely spaced in vac- 
uum. 

Tunneling is not a new idea; it is in fact almost as 
old as quantum mechanics. In less than two years 
after Schrodinger [2] had published his now classic 
papers on "Quantization as a Problem of Proper 
Values," Oppenheimer [3] had used the wave me- 



'Binnig, G.; H. Rohrer, Ch. Gerber, and E. Wiebel, Surface 
studies by scanning tunneling microscopy, Phys. Rev. Lett. 49 
(1982) 57; and Binnig, G., and H. Rohrer, Scanning tunneling 
microscopy, Helv. Phys. Acta 55 (1982) 726. 

2 Figures in brackets indicate literature references. 
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chanicai concept of tunneling to explain the au- 
toionization of atomic hydrogen, Fowler and 
Nordheim [4] had explained the main features of 
electron emission from cold metals in high external 
fields, and Gamow [5] had described the alpha par- 
ticle decay process in heavy nuclei. The efforts by 
Frenkel [6] and by Holm [7] to describe the anoma- 
lous temperature dependence of the resistance of 
electrical contacts were the earliest studies of elec- 
tron tunneling between metals. 

Since these early studies, electron tunneling in 
solids has remained a very active area of solid state 
physics. Discoveries over the last 15 years by Gi- 
aever [8] that tunneling could be used to measure 
the properties of superconductors, by Esaki [9] of 
interband tunneling in semi-conductors, by Joseph- 
son [10] of the unique characteristics of tunneling 
between superconductors, and Lambe and Jakelvic 
[11] that tunneling could be used as a spectroscopic 
tool have had a major influence on the develop- 
ment of physics, and all have important technologi- 
cal consequences. 

In the simplest form, an experiment to measure 
electron tunneling between two metals would be 
performed as shown in figure 1(a). Two parallel 
plates are separated by an insulator and connected 
to a variable voltage source. An ammeter and volt- 
meter are also included to measure the current 
flow and the applied voltage. 'Normally, a direct 
current would not flow through the circuit. How- 
ever, if the plates are spaced very close together 
(on the order of 10 atomic diameters), a tunneling 
current will flow. If the plates could be made suffi- 
ciently flat and parallel over an area of one square 
millimeter, the current flow would be several am- 
peres for an applied voltage of one-tenth volt. 

Conventional devices for conducting these ex- 
periments are fabricated by evaporating one elec- 



trode onto a glass substrate, forming an oxide layer 
with a fixed thickness of 1 to 3 nm on the electrode 
and then evaporating a second electrode over the 
oxide layer. These experiments have been widely 
used to measure the characteristics of Josephson 
junction devices, electrons and phonons in super- 
conductors, semiconductor surfaces and junctions, 
the material used for the insulating barrier, and in- 
tentionally introduced impurities at the barrier- 
electrode interface. While these experiments have 
made and continue to make valuable contributions 
to solid state physics and technology, they have 
been seriously restricted by the small number of 
electrode materials which tend to form natural ox- 
ides suitable for a tunneling barrier. In addition, 
interpretation of the experimental results has been, 
limited because of the difficulties in characterizing 
the atomic structure and composition of the insula- 
tor and the insulator-electrode interfaces. A major 
problem resulting from this lack of structure char- 
acterization is that the energy-momentum relation- 
ship of the electron while it is in the barrier cannot 
be modeled adequately. 

By making a change from these conventional 
devices to one which uses a vacuum barrier, mate- 
rial restrictions are removed, electron properties in 
the barrier are well known, structure of the elec- 
trode-barrier interface is more defined and, in addi- 
tion, the electrode spacing may be continuously 
changed. There are, however, some difficulties 
with vacuum barrier tunneling which must be 
weighed along with its advantages. 

A schematic energy diagram for a metal-vac- 
uum-metal tunneling structure is shown in figure 1 
(b). The energy diagram is based on Sommerfeld's 
electron theory of metals in which all the conduc- 
tion electrons of a metal are considered as being 
placed in a potential energy box and distributed 
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Figure 1-Schematic of a Tunneling Experiment; (b) is from 
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over the allowed energy values according to 
Fermi-Dirac statistics. At zero temperature, the 
electrons, therefore, fill the energy levels up to a 
maximum energy, the Fermi energy. The energy 
required to remove one of the electrons from this 
energy to a large distance from the metal is known 
as the metal's work function. As shown in the fig- 
ure, when a voltage is applied to one of the metals, 
its Fermi energy is lowered relative to the other; 
the energy difference is equal to the product of the 
magnitude of the charge on an electron and the 
applied voltage. 

It will be shown in section 2 that an approximate 
expression for the tunneling current, /, under these 
conditions is given by 



I = Cexp(-aSB U2 )A V, 



(1-D 



where Cs3.16xl0 13 A/mW, a = 10.25 (nm)" 1 
(eV) _1/2 , S is the electrode spacing, B is the average 
barrier height for an electron at the Fermi energy, 
A is the tunneling area, and V is the applied 
voltage. The important consequences of eq (1-1) 
for the vacuum tunneling experiment are that small 
changes in electrode spacing or mean barrier 
height have a very large effect on the tunneling 
current. A change in electrode spacing of 0. 1 nm 
produces approximately a factor of 10 change in 
the tunneling current for a mean barrier height of 5 
eV. Changes of about half this magnitude can be 
produced as a result of work function changes pro- 
duced by the formation of oxide or other films on 
the electrode surface. Thus great care must be ex- 
ercised in electrode preparation and in the mechan- 
ical and thermal design of the vacuum tunneling 
device. 

Young et al. [12] and Thompson and Hanrahan 
[13] have previously performed vacuum barrier 
tunneling experiments. Both of these experiments 
were primarily concerned with observing the tran- 
sition from Fowler-Nordheim field emission tun- 
neling to direct vacuum tunneling. Young et al. 
[12] did report a preliminary observation of a linear 
current-voltage characteristic for direct tunneling 
at a spacing of approximately 1.2 nm. Thompson 
and Hanrahan [13] have not reported any current- 
voltage characteristics. Uncontrolled fluctuations 
in electrode spacing limited the exploration of di- 
rect vacuum tunneling properties in both experi- 
ments. 

The objective of this work, which has evolved 
from that of Young et al., was to extend the previ- 
ous vacuum barrier investigations to a detailed ex- 
perimental study of the direct vacuum tunneling 
region. Such detailed experimental results would 



provide a means to test the present theoretical 
framework of electron tunneling between solids. 
Some of the unresolved questions in this area: 1) 
What is the relationship between the classical elec- 
trostatic image potential and the quantum mechani- 
cal exchange and correlation potential? 2) What is 
the transit time for a tunneling electron? and 3) 
How is the distance between the last atomic planes 
of the two electrodes related to the tunneling bar- 
rier thickness? 

Explicit goals of this work, which were partially 
achieved, were first, to measure the current 
through a vacuum tunneling device for fixed 
voltages of less than one-tenth volt as the electrode 
spacing was varied from a distance of approxi- 
mately 2.0 nm down to a point where the elec- 
trodes touched; second, to measure the 
current-voltage characteristics for fixed electrode 
spacing in the direct tunneling region; third, to 
compare the measured results with available theo- 
ries of electron tunneling; and finally, to establish a 
basis for using a vacuum tunneling device for use as 
a probe to measure the position of a surface to a 
precision of 0.1 nm, a precision 10 to 100 times 
greater than that of presently available techniques. 
Contributions to all these goals were obtained and 
will be described in later sections. The qualification 
on completely achieving the goals was given be- 
cause the quality of electrode surfaces and the pre- 
cision of electrode spacing measurements were not 
adequate in the area of the first three goals to re- 
solve completely the questions posed about the the- 
oretical framework. 

A discussion of some fundamental concepts re- 
lated to electron tunneling in solids, with particular 
emphasis on the origin of the work function and 
the barrier potential energy function, is given in 
section 2. Section 3 describes the requirements of 
the experimental apparatus for a vacuum tunneling 
experiment and the particular solution adopted for 
this work. Experimental results and an analysis of 
the relationship between the results and the numer- 
ical calculations of section 2 are given in section 4. 



2. Some Basic Concepts of Metal- Vacuum- 
Metal Tunneling Theory 

The objective of this section is to calculate, with 
the use of simple models, the current versus 
voltage (I-V) characteristics and, for a fixed 
voltage, the current versus spacing (/ cv — S) charac- 
teristics of the gold-vacuum-gold tunneling device 
used in the experiments of this work. The primary 
results of these calculations are tables and graphs of 
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the l-V and I cv —S characteristics which can be 
used for comparison with the measured data to be 
presented in section 4. 

While quite thorough treatments of electron tun- 
neling in conventional metal-insulator-metal junc- 
tions are available in the literature [14-23], no 
numerical results have previously been published 
which fully accounted for the image potential vari- 
ation with electrode spacing and applied bias 
voltage, both of which are essential for comparison 
with the metal-vacuum-metal tunneling data. Addi- 
tional motivations for performing the calculations 
were that 1) almost all of the available theoretical 
results have been seriously questioned in terms of 
numerical accuracy and overall validity of the ana- 
lytical approximations which were employed [24], 
and 2) even if one were to accept the validity of 
results which have been obtained, the effort to ex- 
trapolate and apply these approximate result^ to 
the gold-vacuum-gold tunneling device would be 
almost as great as carrying out the more accurate 
numerical calculations. 

The content of the following discussion is there- 
fore 1) a description of models of the electron 
properties of the tunneling device and of the tun- 
neling process, 2) a development based on these 
models to obtain analytical expressions for the tun- 
neling current which are, within the limitations of 
the models, "almost" exact, 3) a reduction of the 
exact expressions with commonly used approxima- 
tions to obtain analytical solutions for / vs. V and S 
in terms of experimental parameters, and 4) a de- 
scription of the procedure employed to carry out 
numerical computations of I vs. V and S along with 
a presentation of results from these computations. 
No attempt will be made to give either a review or 
a comprehensive survey of the literature that is re- 
lated to the topics of this chapter. Tutorial cover- 
ages of the theoretical and experimental aspects of 
electron tunneling into semiconductors, metals, 
and superconductors are given in the book edited 
by E. Burstein and S. Lundqvist [14]. A compre- 
hensive review and digest of the tunneling litera- 
ture up through 1969 has been given by C. B. Duke 
[15]. Recent development in the theories of tunnel- 
ing are covered in two papers of debate between 
T. E. Feuchtwang and C. Caroli et al. [21] and 
references contained therein. 

The Work Function and the Image Potential 

By its nature, the metal-vacuum-metal tunneling 
device can be completely represented by a model 
which accounts for the electron properties of two 
metal-vacuum interfaces and their interaction. Ne- 



glecting for the moment the interaction, the prob- 
lem is then to obtain a model that adequately 
describes a uniform metal and its interface to vac- 
uum. In this section properties of the metal-vacuum 
interface will be developed in terms of the general 
theory of metal surfaces formulated by Lang and 
Kohn; LKI [25], LKII [26], LKIII [27], Their the- 
ory will be used as a basis to describe the origins of 
the metal work function and to obtain a potential 
energy function for the metal-vacuum interface. 

Lang and Kohn's work is based on the density- 
functional theory of inhomogenous electron gas 
systems formulated by Hohenberg, Kohn, and 
Sham (HKS) [28]. The HKS formalism is a self- 
consistent-field method and as such its fundamental 
idea is to replace the actual interactions; electro- 
static, exchange, and correlation, between the elec- 
trons and the ion-cores by an interaction between 
one electron and an averaged self-consistent poten- 
tial produced by the other electrons and the ion- 
cores. It therefore fully accounts for the effects of 
electron exchange and correlation and is in princi- 
ple capable of providing exact results. Since the 
concepts of this formalism will be used so exten- 
sively in the following discussion of the work func- 
tion and the surface potential energy function, a 
brief summary of the HKS theory, based in part on 
a review by Hedin and Lundqvist [29], will now be 
given. 

The theory rests on two hypotheses. One states 
that the ground-state wavefunction of an interact- 
ing many-electron system, and thus all ground-state 
properties, are functionals of the electron density 
N(r). The other states that the ground state energy 
for a given external potential U(f) can be written as 

E*{N{r)] = J U(r}N(T)dT+F[N(ir)l (2-1) 

where F is a universal functional of iV, independent 
of the external potential U, and that E u is a mini- 
mum for the correct density function iV(r). 

With these hypotheses as a basis, HKS [28] fur- 
ther show that the total electronic ground-state en- 
ergy of a many-electron system in an external 
potential U{t) can be written in the following form: 



/ 



E u [N]=e \ U(r)N(r)dr 



+ 0/2) J | r _ r M drdi 



+ T S [N]+E K [N]. 



(2-2) 



The functional T S [N] is the kinetic energy of a non- 
interacting electron system, and E„£N] represents 
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the exchange and correlation energy. The first two 
terms in eq (2-2) represent the electron gas' elec- 
trostatic potential energy in the field C/(r) and the 
Coulomb self-energy of the electrons. All the com- 
plexities of the many-body interactions are con- 
tained in the functional jE, c [iV*]. 
An effective potential, C/ efr is then defined as 



UAN; T ] = U(r) + ef -^1-di 



+ UJ[N;tI 



where 



U«[Np} = 8E K [N]/8N(r)- 



(2-3) 



(2-4) 



U e( i[N;r] is then shown to play the role of an effec- 
tive single particle potential so that E a 
= T,[N] + U ca [N;r]. 8 denotes "variation of in the 
calculus of variations sense. This expression for the 
total ground-state energy combined with an as- 
sumed form of EJiN] enables one to solve the fol- 
lowing self-consistency problem to obtain an exact 
density distribution of the system of N interacting 
electrons: 



* 



[_ 2m +%rtAW/= ^ 



iV(r) = ||iMr)re(u.-e,), 



(2-5) 



(2-6) 



m is the electron mass and ft is Planck's constant 
divided by 2ir. The sum is over the N lowest-lying 
orthonormal solutions, $ i} of eq (2-5). Q(x) is unity 
for x <0, and zero for x>0 t u. is the chemical po- 
tential of the metal. The <?, are Lagrange multipliers 
introduced to fix the number of electrons as the 
states ypi are varied. Their identification as single 
particle energies is discussed later in this section. 
At finite temperatures, £/ sc is defined in terms of the 
system's free energy and the function is general- 
ized to the Fermi distribution function. 

The many-body problem has, therefore, been re- 
duced to a set of one-body problems coupled only 
through the effective one-body potential C/ erT [JV;r]. 
Solution of eqs (2-3 to 2-6) proceeds in the usual 
manner: one begins with an assumed iV(r) and 
E K [N(r)], constructs U en [N;r] from (2-3) and then 
finally calculates a new A^r) with the use of (2-5 
and 2-6). The procedure is then repeated until the 
new and old A^(r) differ by acceptable tolerances. 

Lang and Kohn (LKI) calculate the electronic 
structure at a metal-vacum interface by applying 
this formalism to the jellium model of a metal. In 



this model the positive ion-cores of the metal-lat- 
tice are replaced by a uniform, spatially fixed back- 
ground of charge with density 



=0 x>0. 



(2-7) 



The conduction or valence electrons with density 
N(t) are free to move and settle in an equilibrium 
configuration with respect to this fixed background 
of positive charge. N(r) is subject to the charge 
neutrality condition. 



/ 



[JV + (r)-JV(r)]</r=0. 



(2-8) 



For the tunneling experiment this model is a rea- 
sonable approximation since the ion-core lattice 
and its resultant electron band structure are not ex- 
pected to significantly affect the characteristics of 
the tunneling current. 

The electronic structure of an inhomogenous 
electron gas with the inhomogeneity being pro- 
duced by the background of eq (2-7) is calculated 
in LKI. They demonstrate by numerical calcula- 
tions that even though the electron density in the 
surface region of a metal varies quite rapidly, good 
agreement with experimental results is obtained if 
one uses the approximation 



i-/ 



EJLW=J e^N(T)) N{v)dr, 



(2-9) 



e K is the exchange and correlation energy per parti- 
cle of a uniform electron gas. By making use of 
Wigner's form of e w lN) 

<?«= -0.458/r s (A0-0.44/(r s (A0 + 7.8), (2-10) 
^(A0=[3/(47r fl ^)r, (2-H] 

and a U(r) determined by eq (2-7), the functions 
N(x), U a( ^x), and V{x), the electrostatic part of 
f/ e iT, are calculated for several values of the elec- 
tron gas density. In eqs (2-10 and 2-11), e IC is in 
units of Hartrees. (1 Hartree=27.2 eV), c B is the 
Bohr radius, N is the uniform background density 
in number/cm 3 , and r s is the Wigner-Seitz radius 
for a density N. Figure 2 is a graph constructed 
from a table given in LKI of these three functions 
for an electron density very close to that of elec- 
trons in gold, i.e., 5.76 X 10 22 electrons/cm 3 . This 
density corresponds to an r s of 3.03. 

Equations (2-3 and 2-4) in connection with the 
results plotted in figure 2 shows that for x less than 
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5.0 



Figure 2-Graph of Electron Den- 
sity, Effective Electron Poten- 
tial Energy, and Electrostatic 
Component of Potential Energy 
at a Metal- Vacuum Interface. 
The positive background den- 
sity corresponds to a Wigner- 
Seitz radius of 3. The graph is 
based on data from a table in N. 
D. Lang and W. Kohn, Phys. 
Rev. Bl (1970) 4555. 
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about —1.0 nm, C/ efr approaches a constant value 
within the bulk of the material as anticipated. 
Equation (2-5) then simplifies to the one particle 
Schrodinger equation 

[ ~£ V2+ Ve{ ~ "Hj^WOM^/- (2-12) 

Ve(— oo) is the electrostatic potential energy of an 
electron in the metal interior; ^(N) is the ex- 
change and correlation part of the total chemical 
potential jti(A r ) of an infinite uniform gas of density 
N . From the definitions given earlier for E XC [N] 
and e io , and the definition of the chemical potential, 



lx iC (N)=d((NeUN)ydn 



or 



(2-13) 
(2-14) 



=e K {N)+Nd(e^N))/dn. 

If the zero of energy is chosen such that 

V E (~oo)+ix lc (N a ) = 0, (2-15) 

then the equation describing the electron motion in 
the electrodes at distance greater than 1 nm from 
the surface is 



1m 



V^,=^,. 



(2-16) 



Therefore, for the case when iV(r) is a constant, the 
electrons in an interacting gas move as free parti- 
cles. 

Interpretation of the eigenvalues of eq (2-5) as 
the excitation energies of the elementary excita- 



tions of the electron gas is not valid for a general 
U eS j[N;r] [29-31]. However, arguments are pre- 
sented by Moruzzi et al. [31] that demonstrate that 
for most bulk simple metals this interpretation is a 
very good approximation. Slater [30] gives some 
numerical estimates of the energy errors involved 
in using this approximation for self-consistent cal- 
culations in atoms. Hedin and Lundqvist [29] con- 
sider the problem from the Green function 
standpoint and one may conclude from a compari- 
son of eq (2-5) with their equation for the excita- 
tion energies of quasielectrons that for uniform 
positive backgrounds the resulting energies and 
those of (2-5) would be the same. However, 
nonuniformities in the electron gas such as that 
generated by the pseudopotentials of the ion-cores 
and the discontinuity at a metal surface cause all 
the usual differences between one-electron theory 
and many-electron theory to become manifest, i.e., 
self-energy corrections, relaxation times, and dif- 
ferent energy-momentum relationships. 

The work from eqs (2-1 through 2-16) may ap- 
pear an overly complex way of obtaining this free- 
electron result. However, the reader should 
observe that in addition to the free-electron behav- 
ior, this approach has also, by eqs (2-12 and 2-15), 
very naturally given a means to determine the 
binding energy of the lowest energy state. The re- 
sults also clearly demonstrate that the binding en- 
ergy is due both to the electrostatic energy 
produced by the existence of a surface and to the 
energy resulting from exchange and correlation. 
Finally, the development through the calculated 
C4f f has produced a realistic potential energy func- 
tion at a metal surface. 
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The work function, WF, of a metal provides the 
necessary reference for relating V{— «s) to the vac- 
uum, energy level. It is defined as the minimum en- 
ergy required to separate an electron from the 
metal at zero temperature. Thus, 



WF=(^(+ «)+£„_,)_£■* 



(2-17) 



where V E {+ oo) is the total electrostatic potential 
energy far from the metal surface and E M is the 
ground state energy of the metal with M electrons. 
Arbitrary zeros of energy for Fand E do not affect 
the :esuh of eq (2-17). The chemical potential /x of 
a system is defined as E N —E N -i, thus eq (2-17) can 
also be written as: 



WF=Vz(+<»)-}x. 



(2-18) 



Since the electrostatic potential in the interior of a 
metal with and without a surface differs by 
V(— oo ), u. for a metal with a surface is 



ji=F E (- «)+£„+ n™CW.), 



<2-19) 



E P is the Fermi energy of a degenerate electron gas 
of density N a . WF is therefore 



fVF=lV £ (+°°)~Y E (-°o)] 
In a more compact form, 



(2-20) 



WF=LV E -y. ! (2-21) 

where AV E =V E (+ <»)— p>(— oo) and 

ji=£ F -4-p,«(JV a ). (2-22) 

In this form the WF is seen to consist of two com- 
ponents, a surface part AJK^, which is the rise in 
mean electrostatic potential energy across the 
metal surface and jX, a bulk part which is the intrin- 
sic bulk chemical potential of the electrons relative 
to the mean electrostatic potential energy in the 
bulk interior. LV E is due to the dipole layer created 
by the spilling-out of the electrons as shown by the 
distribution in figure 2. All many-body effects are 
contained in the exchange and correlation contri- 
butions to ji and in their effect on the barrier poten- 
tial &V E . 

Combining eq (2-20) with the zero energy con- 
vention of eq (2-15) shows that the lowest energy 
eigenvalue of eq (2-16) is at an energy of 
[AKe + u^W,)] less than P(+<»). Eigenstates of 
eq (2-16) are then according to the Pauli exclusion 
principle occupied at zero temperature from this 
energy up to the maximum value E F . A diagram of 
the electron energy as a function of x position in 
the region of a metal vacuum interface is given in 
figure 3. 

The work function of gold was estimated in 
LKU by again using eq (2-13) and Wigner's ex- 
pression for e lc (N) to obtain fi; AV E was calculated 




r 



V(+») 



VH-VW WQRI< 

vr - • v FUNCTION 



Figure 3-Energies pertaining to work function definition. Rela- 
tive energy scale was adjusted to give ratios of the quantities 
indicated which were that obtained by Lang and Kohn for 
gold. 
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from the double-layer solution of Poisson's equa- 
tion 



The result for a charge located at the position xi is 
to order (xi— x ) 3 given by 



A V E = (e 2 /£o) / x[N{x)-N + (x)]dx (2-23) V w (x) = const - ?7[4ire (x. -x )] 



and the density distributions obtained from self- 
consistent solutions of eqs (2-5 and 2-6). In eq (2- 
23), e is the electron charge and e is the free space 
dielectric constant. The results for gold were (all in 
units of electron volts) 



^=5.57,^= -6.75, AV E 
= 2.32, iVF = 3,.50, 



(2-24) 



Experimental values of the work function are [32] 
5,0+0.2 eV for polycrystalline samples; the uncer- 
tainty accounts for the range of values obtained by 
various experimental techniques. Attempts to 
achieve better agreement by incorporating lattice 
effects through a pseudopotential were not success- 
ful. 

It is not clear how the discrepancy between the 
theoretical and experimental work functions 
should be apportioned between u. KC and LV E . But 
regardless of this difficulty, Lang and Kohn's cal- 
culations have shown in a definite way the relative 
magnitudes of the two major contributing terms to 
the work function. In addition, this insight into 
how the work function energy is generated will 
later prove useful in interpreting the effects of ab- 
sorbed molecular species on the values of the work 
function. 

There remains one further significant question 
about the electron structure at a metal-vacuum in- 
terface: For purposes of a tunneling experiment, 
what is the potential energy of an electron 1 to 10 
nm outside a metal surface? A related question 
which must also be answered is: Where is the effec- 
tive metal surface? Both questions are addressed in 
LKIII. By solving eqs (2-5 and 2-6) self-consis- 
tently, an N(t) is calculated from a £/ eff determined 
by the sum of the energies due to the distribution of 
eq (2-7), a point charge outside the metal surface, 
and a weak external electric field F. The unper- 
turbed charge distribution produced by the charge 
density of eq (2-7) is then subtracted from this re- 
sult to obtain a screening charge density of N s (x). 
The potential energy of an electron outside the sur- 
face in the presence of a weak uniform electric 
field pointing out of the metal is then calculated as 
the interaction energy between the potential of the 
point charge, q , and its induced screening charge. 



-gE(x -xj. 



(2-25) 



The center of mass, x , of the induced screening 
charge is defined by 



x = JxN s (x)dx /§ N s (x)dx 



(2-26) 



For small fields, the center of mass of the induced 
screening charge is approximately the same as that 
for the point charge. The answer to the first ques- 
tion posed is, therefore, that the potential to a good 
approximation is the classical image potential re- 
ferred to a plane passing through jc . Relative to the 
edge x B of the uniform positive-charge back- 
ground, the position x a is determined by 
xo— Xb =0.076 nm for an r s = 3, i.e., gold's electron 
density. This value was obtained by interpolating 
between the values of x ~- x B given in LKIII for 
r s = 2 and r s —A. x and x B are located schematically 
in figure 4. 

In connection with the problem of determining 
the potential at a metal surface, field penetration 
effects are often considered to be a significant fac- 
tor. As pointed out by Lang and Kohn, the analogy 
here of penetration effects is just the fact that the 
induced charge distributions have finite widths. 

A definite answer for the second question is not 
as easily obtained. Locating x B of the positive back- 
ground model relative to the lattice of the metal is 
the problem. The approach taken in LKIII was to 
locate, somewhat arbitrarily, x B at a point one-half 
an interplanar spacing in front of the last plane of 
ions. A justification for this assignment can be ob- 
tained from the self-consistent calculations by Ap- 
plebaum and Hamann [34] for the Na (100) surface 
in which the lattice properties and atomic nature of 
the surface were treated at the same level of detail 
as is common in bulk calculations. They give a fig- 
ure which shows the effective potential and charge 
density in much the same manner as figure 2, with 
the ion-core locations shown. Taking the point at 
which the electron charge density has dropped to 
one-half the bulk value as the proper place to lo- 
cate the positive-background edge gives the x B 
value assumed by LKIII. 

In summary, the electron potential outside a gold 
surface is given by eq (2-25) with x located at a 
point one-half an interplanar spacing plus 0.076 nm 
in front of the last plane of ions (see fig. 4). This 
along with C/ efr is sufficient information to deter- 
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| Effective 

I Metal Surfaces 



c 



fc-1/2 d 



0.075— ?| |4— 



ELECTRODE 



ELECTRODE II 



Figure 4~Location of effective metal surfaces. Diagram illus- 
trates simple model of two identical gold electrodes with the 
atoms spaced at the nearest-neighbor distance in a gold crys- 
tal. The electrode spacing shown would correspond to a sepa- 
ration of the respective last ion planes of approximately 1.73 
nm; the nearest-neighbor distance was assumed to be 0.288 



x s =Xj+0.5d +0.076 nm+r, 



(2-28) 



where t is the distance from x to where C/' eff =£ F ; 
thus, / =0.105-0.076=0.029 nm. The positive x 
direction is out of the metal surface. 

At large, positive x , U x [N(t)] is thought to be- 
come the classical image potential [U xc (x)—>(e 2 / 
4n Co*)] [25,33]. However, no form of the 
exchange and correlation potential has as yet been 
put forth which has the expected asymptotic func- 
tional dependence on x . The transition of an elec- 
tron from inside a metal in the presence of many 
other electrons where energy changes result from 
electrons of like spin avoiding each other (ex- 
change) and electrons of either spin avoiding each 
other to minimize the coulomb energy (correlation) 
to a region in vacuum where at least the correla- 
tion is an insignificant factor has been extremely 
difficult to model. The slowly varying density ap- 
proximation for e KC [N] given in eq (2-9) leads to an 
effective exchange and correlation potential which 
vanishes exponentially as x goes to infinity (see 
LKI). A control calculation was performed in LKI 
by checking the self -consistent results against that 
obtained when the image potential (1/x) form was 
used for x values greater than the point where 
U c(i =E P . Differences in N(x) obtained with the 
two potentials were no greater than 0.3% for 
^=2.5. The more commonly used Slater X-a ap- 
proximation for the exchange and correlation po- 
tential varies as N(x) m . With reasonable N(x)'s it 
also vanishes much more rapidly than l/x as x goes 
to infinity [34]. 



mine the potential between two gold electrodes in 
vacuum. 

However, the effective spacing between the 
electrodes for tunneling is reduced even more than 
twice the sum just given since the tunneling dis- 
tance will, as shown later, be determined primarily 
by the distance between the points of the elec- 
trodes at which U e! r=E P . For gold (see fig. 2) the 
effective tunneling surface, by this convention, 
would be located 0. 105 nm plus one-half an inter- 
planar spacing in front of the last plane of ions. 
This is the distance obtained from figure 2 at which 
U c!f =0; x a in figure 2 is at the ordinate 0. For an 
interplanar spacing, d, the location of the gold sur- 
face for tunneling, x„ relative to the location of the 
last plane of ions, xi, is therefore 



X s =X t + (X B ~Xj) + (x -Xb) + 1 



(2-27) 



or, 



The Barrier Potential Energy Function 

In this section and the following section on the 
tunneling process, the tunneling device is assumed 
to be composed of two semi-infinite metal elec- 
trodes with plane parallel faces separated by a vac- 
uum space. The parallel faces are taken to be 
normal to the x axis with electrode II located fur- 
ther along the positive x direction than electrode I. 

Assume first that the electrodes are in equi- 
librium with each other. Then, the chemical poten- 
tials of the two electrodes will be equal [35] and eq 
(2-19) gives 

V l E {- oo)+E F +jil(Nl)=V P X- oo)+Eh l 

+t&(N?). (2-29) 

The superscripts denote the electrode labels. When 
the electrodes are not contacting this condition 
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combined with the definition of the work function 
in eq (2-20) requires that 

V V = V\{+ oo)-F£(+ «)= fFF 1 - WF n . (2-30) 

Thus, if the two electrodes have different work 
functions, the equilibrium condition is for there to 
exist an electrostatic potential difference between 
their surfaces. This potential difference is usually 
known as the Volta or contact potential [36]. Re- 
cent work by Bennett and Duke [37] has, however, 
shown that the name contact potential would be 
inappropriate since they conclude that when the 
metals are in contact, the junction potential, V h is 



V^El+vUNb-lEg+yUN?)]. 



(2-31) 



This result interpreted in terms of eq (2-29) means 
that when the electrodes are in contact, their back- 
ground electrostatic potentials, V{— oo), are con- 
strained to be equal; a very reasonable conclusion. 
Bennett and Duke's [37] arguments are not based 
on eq (2-29). Their results follow from solutions to 
Schrodinger's and Poisson's equations with appro- 
priate boundary conditions. 

The problem of describing the transition be- 
tween the conditions specified by eqs (2-30 and 2- 
3 1) has, unfortunately for use in a general vacuum 
tunneling configuration, not been investigated. For 
the particular configuration, i.e., with almost iden- 
tical gold electrodes, used in the experiments of 
this work the effects from both V } and V v were 
minimal. In the following discussion, V r and V } are 
assumed to be zero for the gold-vacuum-gold tun- 
neling experiment. 

Assuming equal Fermi energies and work func- 
tions in the two electrodes, the barrier potential 
energy function, V B (x) at equilibrium is, following 
the conventions of figure 3 given by 



V B (x)=E F +WF+V iP (x), 



(2-32) 



where V lf (x) is the image potential of a point 
charge between two conductors with plane parallel 
faces. The energy zero is assumed to be such that 
V\{— oo)+n4(flj) = 0. The classical image poten- 
tial (2-25) is employed because of the difficulties 
given earlier about the incorrect extrapolations of 
available forms of the exchange and correlation po- 
tentials. In vacuum, the classical double image po- 
tential energy function is [38] 

M*)=-5^-(?-+ 2 2 f , --), (2-33) 

8i7-£ v 2x ,i=i,„ rrs — x 1 ns / ' 



where s is the distance between the faces of the two 
conductors and x is the coordinate relative to an 
origin at the face of electrode I. To facilitate nu- 
merical calculations, a simpler approximate form of 
eq (2-33) obtained by Simmons [38] will be used. 
The approximate form is 



F IP (x)=-1.15(ln2)e 2 



8ire<yc(,y— x)' 



(2-34) 



With s and x expressed in nm and K ]P (x) in electron 
volts (eV), eq (2-34) reduces to 



V lf (x) = -0.2867 



x(s — x)' 



(2-35) 



Plots of V B (x) with the V lv of eq (2-35) and as- 
sumed electrode Fermi energies of 5.5 eV and 
work functions of 5.0 eV for s values ranging from 
0.2 nm to 5.0 nm are given in figures 5 and 6. Note 
that in these figures, the vacuum level is at an en- 
ergy of 10.5 eV. These graphs demonstrate that as 
the electrode spacing decreases, the maximum bar- 
rier height drops by more than one electron volt 
when s = 1.0 nm; when s =0.23 nm, the maximum 
barrier height has dropped to the Fermi energy of 
5.5 eV. This is coincidentally about equal to the 
nearest-neighbor distance for the gold lattice. 

For the barrier potential energy function plotted 
in figures 5 and 6 to correspond with the surface 
location described by eq (2-28) and the corrections 
to the classical image potential given in eq (2-25), 
the x =0 and x =s positions would correspond to 
the x described in these equations. While t, the 
distance from x to the x value at which U ea =E F , 
was found in LKIII to be 0.029 nm for gold, the 
value obtained from figure 5 is approximately 0.06 
nm. Neglecting for the moment the possible effects 
of nonzero t, the modifications of eqs (2-25 and 
2-28) mean that tunneling currents calculated from 
the barrier potential with a spacing parameter s 
should correspond to a true physical spacing, s T , 



s r =xj I —xf =s +d +0.152 nm. 



(2-36) 



The importance of eq (2-36) will be discussed in 
section 4 when an attempt is made to relate the 
calculated current values to measurements. 

When the positive terminal of a battery of poten- 
tial V is connected to electrode II and the negative 
terminal to electrode I, the mean electrostatic po- 
tential, for negative charges, of electrode II is low- 
ered by an energy of eV relative to that in 
electrode I. Since the E F and /a, c are fixed by the 
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Figure 5-Barrier potential energy 
vs. distance for electrode spac- 
ings of 0.2 nm to 1.0 tun. 
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Figure 6-Barrier potential energy 
vs. distance for electrode spac- 
ings of 1.0 nm to 5.0 nm. 
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electron density, the Fermi level is also lowered by 
an energy eV. Combining the barrier potential en- 
ergy function with this effect and assuming the 
field between the electrodes is uniform gives a po- 
tential energy function for the metal-vacuum-metal 
system which can be subdivided into three areas of 
interest: 



Region I: F T (x)=0 



*4 X 



x n >x 



Region II: V T (x)=-tV 

Region III: V B (x)=E P +WF~eVx/s 



+ V 1! (x)x 1 >x>x n , (2-37) 
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x l and x 11 are the values for which V B (x) is equal to 
and — eV, respectively. Figure 7 is a schematic 
illustration of eq (2-37) combined with the con- 
cepts of figure 3 concerning the work function. 
The areas of the diagram up to the Fermi energy 
represent eigenstates that are occupied when the 
electrodes are at zero temperature. 

The Tunneling Process 

From this point onward in the theoretical devel- 
opment no further attempt will be made to account 
for many-electron effects on the calculated tunnel- 
ing current. All the previous work to account for 
many-electron effects was required since, as 
shown, the dominant sources of the work function 
originate from electron-electron interactions and 
their effects on the electron distribution at a metal 
surface. A fully independent-electron modet-of the 
electrodes and of the tunneling process will be used 
in the remainder of this section. There were two 
major reasons for this decision. First, at the resolu- 
tions of J'V measurements anticipated for this ex- 
periment, the small changes in current due to 
electron-electron interactions, electron-phonon in- 
teractions, and self-energy changes produced by 
transition through the barrier would not be observ- 
able. The changes are usually less than one percent. 
Second, within the limitations of the independent 
electron model, the stationary state method of cur- 
rent-carrying states is a well-defined and conve- 
nient procedure for calculating current flow in 
quantum systems [39,40]. 

For an electron described by one of the solu- 
tions, iK(r), to the time-independent Schrodinger 
equation for the potential V r (x) 



[~^+ ^)]W')=E(k)#). (2-38) 



J(k,r)=-^ 0/f k Vi^-^VM 



(2-39) 



where * denotes the complex conjugate operation. 
The standard interpretation of eq (2-39) is that J-n 
is the probability per unit time that a charge — e 
will flow through the unit area with normal n. 
Since the electrons occupying the states i|* t (r) are 
assumed to not interfere with each other, the total 
current density is just the sum of the current den- 
sity contributions from each eigenstate of the 
metal-vacuum-metal system which is allowed by 
the Pauli exclusion principle. 

The non-equilibrium aspects of a biased junction 
are usually treated by assigning different thermal 
occupational probability functions to each metal 
with the reference energies offset by the applied 
voltage and differences in Fermi energies. Thus, 
with V\— oo)+u4 as the reference energy zero 
and a positive polarity voltage, V, connected to 
electrode II, the energies in the two electrodes 
would be related by the energy conservation equa- 
tion: 



E^Eu+Epj-EFu-eV. 



(2-40) 



the electron current density flow is given by [40] 



For a net current to flow, occupied states on one 
side of the barrier must exist at the same energy as 
empty states on the other. In figure 7, the applied 
voltage will establish, at zero temperature, empty 
states in electrode II opposite filled states with the 
same energy in electrode I. Therefore, in this in- 
stance, electrons would dominantly flow from 
electrode I to electrode II by tunneling through the 
barrier, i.e., via a current-carrying stationary state 
which existed in all three regions of the system. 
The current density is computed as the difference 
between the current density flowing from I to II 
and II to I. 

The thermal occupational probabilities are given 
by the Fermi-Dirac distribution 
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Figure 7-Energy Diagram for 
Metal-Vacuum-Metal Tunnel- 
ing. Energy and distance in the 
diagram are drawn to scale for 
the quantities represented. The 
small oscillations of the energy 
at the bottom of the occupied 
states represent the magnitude 
and spacing of the ion-core po- 
tentials. X\ and Xi are the classi- 
cal turning points for an 
electron near the Fermi energy. 
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/<i>)=(l + exp[^]r 



(2-41) 



where E is the energy eigenvalue, k e is Boltz- 
mann's constant and T is the system temperature. 
Equation (2-40) combined with eq (2-41) gives the 
result 



ME)=f(E) 
ME)=f(E+eV). 



(2-42) 



Again, taking the plane of the electrodes to be nor- 
mal to the x -direction, the current density, J^n, 
from I to II is then given by 



J^„=22 Jf^ry,^)!! -ME)]. 



(2-43) 



The factor of 2 is to account for the additional 
states produced by the two electron spin orienta- 
tions. Jf_.ii is the x component of J(k,r). Similarly, 






(2-44) 



Since V T {x) is a real potential, the solutions to 
Schrodinger's equation must obey the general re- 
ciprocity theorem of scattering processes. Blatt and 
Weisskopf [41] give a succinct statement of the the- 
orem: "This theorem states that the probability for 
a transition proceeding one way In time is equal to 
the probability for the same transition, but with the 
sense of time reversed," The sense of time in eqs 
(2-43 and 2-44) is manifest only in the direction of 
the wave vector, thus, Ji-.a(k,r)=J[^ a (—k,T) 
= Jf]„i(k,r), The expression for the net current den- 
sity is therefore 



J=2E J*(k,r) \f{E)-f{E +eV)]. 



(2-45) 



To calculate J x (k,r)=J*(k,i) (a subscript nota- 
tion is adopted hereafter), explicit solutions to eq 
(2-38) for all three regions of Vt{x) must be ob- 
tained. The solutions given here will be based on 
the results obtained by Merzbacher [42] with the 
use of the Wentzel, Kramers, Brillouin (WKB) ap- 
proximation technique for solving Schrodinger's 
equation. If periodic boundary conditions are im- 
posed in the y and z directions, the solutions to eq 
(2-38) may be written as 



^)=A ~ ]/1 exp(; V*W*,C*)> 



(2-46) 



where k^ and P are the wave vector and displace- 
ment, respectively, parallel to the electrode faces, 
k x is the wave vector component normal to the 



electrode faces, and A is the area of the tunneling 
device; M ks (x) is defined below. The eigenvalues 
E(k) are 



T-1A 



w^c*,+*-)=^ 



(2-47) 



Let jc l (JEi) and x R (E) be the classical turning points 
for the energy E at electrode I and electrode II, 
respectively, then with Merzbacher's conventions 

M kx (x)=Ak- 1 <\x)VJ + 

+Bk~ x! \x)W_ x£x L (E) 

M kx (x)=CK- in (x)Z_ 

M kx (x)=Fk- u \x)q + 

+ Gk~ in {x)q_ x^E^x, <2-48) 

where 

W±=texp[±iJ^ k s (x)dx] 

Z ± =exp[±( K(x)dx] (2-49) 

"L 

q ± =exp[±ij k x (x)dx] 
and 



*R 



K(x)=\^{V B (x)-{E{k)-^)}Y-. (2-50) 



2 k 2 - 
JEjl 

2m 



With these conventions, Bohm [43] then shows in 
detail that with the use of connection formulas 



^=[(1/2X26+1/(20)]^ 

where 



(2-51) 



e=exp(f K{x)dx). (2-52) 

If now the solution in region I is assumed to be 

Mi t (*) e ik * x +R e ~ ik ' x (2-53) 

with R being adjusted for normalization, then the 
solution in region II is, assuming G=0, i.e., no 
wave is incident from right, 

Ml (x)=2[29 + 1/(28)]^' e'^"--^> e'M (2-54) 

The conservation of k p and total energy in the tun- 
neling process was implicitly assumed in the calcu- 
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lations from eqs (2-46 to 2-54). Since M" (x) is the 
stationary state solution in region II and probability 
current is "conserved," it may be used as well as 
the solution in the other two regions to calculate 
the current density. Substituting eq (2-54) into eq 
(2-46) and then the result into eq (2-39) gives 

/,(k,r)=4[2e+ 1/(20)] - 2 eHk x /m . (2-55) 

Let 

Z>( J B,k,)=4[0+l/(2e)]- 2 , (2-56) 

then 

J x (k,r)=D(E, k p ) efikjm, (2-57) 

D(E,k P ) may be interpreted as a transmission coef- 
ficient for the part of the solution in region I into 
region II. 

Substituting eq (2-57) into eq (2-45) gives the 
complete expression for the current density / as 



J=2eWm 2 k x D(E, k p ) [f(E) 
-/(£+eV)]. 



(2-58) 



The summations are most conveniently carried out 
by changing to integrations over the variable E x 
defined as 



E x =H 2 k 2 x /2m =E(k)—n 2 k 2 /2m . 



(2-59) 



The conversion is performed with the following 
three relations: the first is the well known sum to 
integral conversion in k -space, the second one is 
based on the free-electron properties expressed in 
eq (2-47) and the third is the differential of eq 
(2-59): 



■!<■ 



, /c/(2tt) 3 



d 2 k p = (2irm/fi 2 )dE p 
kJk^im/^dE,, 



(2-60) 



With these relationships, eq (2-58) may be written 
in the following form: 

J=f o dE x D{E x )N{E x ) (2-61a) 

NiE^iAirme/h')^ \f(E) 

~f(E+ e V)]de p . (2-61b) 



In words this states that the net current density is 
equal to the sum of the number of electrons with 
the "jc-part of their energy" within dE x that are 
incident on the surface of the negatively biased 
electrode per unit time per unit area, N{E x )dE x , 
times the tunneling probability of an electron at 
this energy, D(E X ). 

This is the same form of the equations obtained 
by Duke [44] and as such is an acknowledgment 
that his work was frequently consulted in develop- 
ing eqs (2-61). Another form of eq (2-6 la) which is 
more useful for interpreting the energy aspects of 
the tunneling process is 



J=^fdE[f(E)-f(E+eV)] 
f2ir- 2 d%D(E,k f ). 



(2-61c) 



Equation (2-6 lc) facilitates the physical interpreta- 
tion of the origins of the tunneling current in terms 
of the overlapping of energy surface projections 
onto the tunneling plane and enables direct connec- 
tion with energy diagrams such as that given in 
figure 7. However, calculation of the tunneling 
current, particularly for the free-electron model of 
the electrodes, is most easily accomplished with 
the use of eq (2-6 la). 

Equation (2-6 lb) may be integrated by parts to 
obtain 

^(£^=(4 mek B T/h z ) In {FDl/FDll), 

FD\ = 1 + exp [(E F - E x )/k B T] 

FDll= 1 +exp [E F -E X ~eV)/k B T]. (2-61d) 

Note that the Fermi energies of the two electrodes 
have been assumed equal. Collecting all the defini- 
tions involved in eq (2-56) together for reference 
gives 

^(-EJ=4[2e-f-l/(2G)] 2 , 

e=exT?(fj L K(x)dx) t 

K(x)=[2m /K 2 )(V B (x) ~E x )f /2 , 
V B (x )=£ F +WF- eVx /s 

— (const)? [x(s —x)] ~ \ 

const = 1 . 1 5(ln2)e 2 /Sv€ Q 

=0.2867(eV nm). (2-61e) 

The accuracy of the WKB approximation for 
calculating the transmission coefficient has been 
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extensively investigated in the literature of this 
field. One of the most thorough checks was per- 
formed by Politzer [45], who compared values 
computed with the form 



-D(£)=exp[-2£ K(x)dx], 



(2-62) 



of the WKB transmission coefficient with those ob- 
tained from a numerical solution of Schrod inker's 
equation using essentially the same V T as given in 
eq (2-37). Her results show that the WKB result is 
valid at the fields and energies involved in most 
configurations of the vacuum tunneling experi- 
ment, i.e., fields less than 10 8 V/cm and energies 
much less than barrier height. 

The major difference between the exact results 
and the WKB approximation was that the exact 
D(E) was an oscillatory function of energy for en- 
ergies and fields which resulted in an effective 
bound state potential well above the top of the bar- 
rier, The amplitude of the oscillations in the exact 
D(E) was relatively small and such that 

0.5 <Log(z>(£) WKB /i?CE) NUM )«X8 (2-63) 

for all fields and for energies more than 0.5 eV 
from the top of the potential barrier. Near the 
Fermi energy, at which the major contribution to 
the tunneling current occurs, the ratio's mean value 
was about —0.2. The mean difference between the 
tunneling current for the two calculations of D(E) 
was not calculated by Politzer, and barrier heights 
applicable to the vacuum tunneling experiment 
were not investigated. 

Approximate Analytical Expressions for Tunneling 
Current 

In this section, the qualitative properties of the 
tunneling current will be derived from the zero 
temperature, low voltage solutions to eqs (2-61 a to 
e). The overall character of the tunneling current 
will be shown to be approximately ohmic with a 
small cubic dependence on voltage having impor- 
tance at higher voltages. At all voltages, the tun- 
neling current is found to be exponentially 
dependent on the electrode spacing. Finally, the 
zero temperature solutions are extended to finite 
temperatures and estimates of the effects of room 
temperature on the tunneling current are given. 

The tunneling characteristics determined by eqs 
(2-61) are most easily understood by first investi- 
gating the zero temperature behavior. At zero tem- 
perature 



N(E x )=A T zV, Q^Et4Er~eV 

=A T (E F -E !: ), E F -eV<E x <E f 

=0 E F $E X (2-64) 

A T is the "thermal" constant 4tt me/h 3 . Thus, eq 
(2-61a) takes the form 

J=A T [eVJ D{E x )dE x 

(2-65) 



+ J* F - e v {E ¥ ~E x )D{E x )dE x l 



For large G in eq (2-6 le), which for most barriers 
in tunneling experiments is a true condition, D{E X ) 
can be approximated by the form given in eq (2- 
62). Simmons [46] has proposed an average barrier 
height method to further simplify D(E X ). His ap- 
proach is to approximate the integral in eq (2-6 le) 
over K (x) with an average barrier height, B, above 
the Fermi level of the negatively biased electrode, 
i.e., he defines B as 

Jf»*RP 
etp V B &)dx, (2-66) 

where jc L f and Xkp are the classical turning points 
for E X =E F , and W^x^p— x LP is the width of the 
barrier at this energy. 

With this definition, Simmons [46] shows that for 
voltages of less than about one volt 



D(E x )=ex.p[-a1V(E F +B -E x f n ] 



(2-67) 



fl=(47r/h)(2m) 1/3 =10.24(eV)- 1/2 nm" 1 . By utiliz- 
ing this D(E X }, an approximate low voltage tunnel- 
ing current density can be calculated from eq 
(2-65). Let N=E F -E X sadp^aWB m , then eq (2- 
65) may be written as 

J=A T {-eVJ B exp[-p(\+n/B) m ]dN 

-J Nexp[-p(l +N/B) vl ]}dN. (2-68) 

For small voltage, N/B, will be much less than 
one, so the approximation (1+N/B) l/2 = (1+N/ 
2B) may be used. The integrals may then be evalu- 
ated to obtain the result: 



J =A r [(^e- p /(aWf) 

X ( 1 - exp[-a WeV/(2B ,/2 )])]. 



(2-69) 



Only a factor of exp( — a tVEp/B l/2 ) was neglected 
to obtain eq (2-68). If one takes just the leading 
term in the series for the exponential, the following 
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result for the tunneling current density at low bias 
voltages and zero temperature is obtained; 



/ L (7-=0)=/ OL exp(-p)F, 



(2-70) 



with 



y 3 =(7r 2 me/h 2 )=3.27 V~' nm" 1 , 



(2-73b) 



where J aL =[{2m) m /W]{e/hf 3 m , If B is ex- 
pressed in eV and W in nm, J 0L is 

/ 0L =3.16X \Q n B m W~ l (A/m : /V). {2-11) 

By evaluating eq (2-68) in a more exact manner, 
Simmons [46] derives a general formula for the cur- 
rent density at zero temperature. The formula and 
its associated relations are given below in what has 
been found to be a consistent algebraic form, which 
due to misprints, etc., is not correct in published 
forms. It is given here because it will be used for 
comparing with the numerical results of the next 
section. With the same definitions given in eqs (2- 
66 and 2-67), the formula is 

J(T=0)=J o (B exp(-aWB m ) 

-{B +eV)exp[-a^(£+eV) ,/2 ]), 

= 6.2 X \Q ll W~\A nmy<m 2 • eV). (2-72a) 

Based on eq (2-66), the analytic approximation for 
B including image potential effect is 

B = WF-(V/2S)(S l +S 2 )-[Q.2&6/(S 2 ~S 1 )] 

XMSi('S-Siy(S l (S-S 1 ))] (2-72b) 

S= electrode spacing in same sense as previously 
discussed. 

S i =0.3/ WF 

S^Sil-l.yQWFS + lQ-lVS^+S^. 

In eq (2-72b), the work function should be ex- 
pressed in electron volts and all distances involved 
in nanometers. 

Equation (2-72a) reduces quite easily to eq (2- 
70) by expanding both exponentials in a series and 
neglecting terms containing V 1 and higher orders. 
However, if one keeps all terms out through V i , it 
is found that the quadratic terms cancel to give the 
next higher order approximate result [19], For 
identical electrodes 



J(T = Q)=J ll dy+(J 3 S 2 /B)V 3 ], 



(2-73a) 



This value of j\ was taken from the results of 
Brinkman et ah [19] who expanded eqs (2-61) in 
powers of the voltage and obtained approximate 
expressions for the coefficients of the linear and 
quadratic terms of the conductance. They consid- 
ered eq (2-73b) to be accurate to within 10% for 
trapezoidal shaped barriers and where S is greater 
than 1 nm. An earlier value obtained by Simmons 
[47] for the same conditions was^ 3 = 1.15 V -1 nm -2 . 
No calculations of ji with image potentials correc- 
tions are available in the literature. In spite of the 
relatively unknown values for this constant, the 
qualitative behavior predicted by eq (2-73a) has 
been fairly well established in many different ex- 
periments [48,49,19]. 

In summary, the ideal tunneling characteristics 
are then as follows: 

1) Current density varies exponentially with 
electrode spacing, with an exponent of —10.25 X 
mean barrier height in electron volts X electrode 
spacing in nanometers. 

2) The current-voltage characteristic should be 
symmetric for identical electrodes. (This just 
serves as a check on the formalism or the experi- 
mental method since it must be true if the elec- 
trodes are truly identical.) 

3) For low voltages, the current-voltage charac- 
teristic is approximately ohmic. 

4) For voltages up to the work function, the non- 
ohmic conduction should be dominantly a cubic 
function of the applied voltage, with higher order 
terms becoming more important as the voltage is 
increased up to the work function. 

The extension of the previous results to finite 
temperatures is the last topic to be discussed in this 
section. If the approximations which led to eq (2- 
68) are substituted into eqs (2-61), the following 
result is obtained: 

J{T) =Ajk B T exp(-a WB m ) 

y-J E ln(FDI/FDir) 

exp( - a WN/2B m )dN, (2-74) 

where again N=E F ~Ej,. When k B T is less than 
2B U2 /aW, the integral may be evaluated by the 
techniques of Murphy and Good [50] to obtain 



187 



Journal of Research of the National Bureau of Standards 



a J{T)=A 1 Q- 2 \Qk B T/$m{iTQk II T)]Xexp 

(-p)[\ -exp(-0eV)] 5 (2-75) 



where Q={\/2)aWB- vl . For r approaching zero, 
the factor in the first bracket goes to one, and Sim- 
mons [46] shows that the remainder of the terms 
are equal to J(T=Q), Therefore, 



J(T)/J(T =0)=Qk B TMiL(jTQk 3 T). 



(2-76) 



Q is typically less than 4(eV) _1 . The argument of 
the sine is therefore much less than one at room 
temperature. The series expansion gives 

J(T)/J(T = 0)= l+(l/6)(i;Qk B Tf 

or in practical units (2-77) 

J(T)=J(T =0)(l + [3 X 10' "(eV/nmVJST 2 ) 

(WTf/B]). 

For W=\ nm, T=300 K, and 5=2 eV; 
J(T)= 1.013 J(T=0). Such small changes in the 
tunneling current would not be detectable with the 
present experiment No significant effect from con- 
ducting the experiment at room temperature was 
therefore anticipated. 



Numerical Calculations of Tunneling Current 

In order to facilitate the analysis of measured I- V 
and I cv —S tunneling data, the following calcula- 
tions, based on eqs (2-61) were performed. The re- 
sults are given in accompanying figures and tables. 

1) The product of D{E X )N{E X ) as a function of 
E x for four voltages; 0.001, 0.010, 0.10, and 1.0 V, 
at each of three electrode spacings; 5 = 0.3, 1.0, and 
2.0 nm. See figures 8, 9, and 10. 

2) For a simple model of the tunneling between 
spherical electrodes, described below, calculations 
of the tunneling current density and current were 
made for a fixed bias voltage of 0.1 V as a function 
of electrode spacing and work function. Computa- 
tion times were long for the current calculations, 
15 seconds per data point, so only a limited number 
were obtained. These results are given in table 1 
and Figures 1 1 through 17. The graphs in figures 11 
and 12 were at first plotted in a normalized form, 
but after finding that much information was lost 
because of the large orders of magnitude changes 
between graphs, the four separate graphs were 
chosen as the more useful representation. Figure 1 1 
shows the changes in curvature or lack thereof in 
the I-Vfoi a spacing of 2000 pm. Figure 12 shows 
the increase in curvature of the 1-V characteristic 
as a function of electrode spacing for a fixed work 
function of 4.5 eV. Figures 13-17 are graphs of the 
data in table 1 and enable a comparison of the cur- 
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Figure 9-Narmal-eitergy distribu- 
tion of tunneling electrons; elec- 
trode spacing of 1.0 nm. 
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rent and current-density so that the tunneling area 
may be estimated. 

3) The tunneling current density for a wide 
range of voltage and electrode spacing; 
V=N X 0.05 volts, N = l to 12 and S=(f)3+M 
X0.1)nm, M =0 to 27. See table 2. 



4) As a cross-check on the correctness of the 
numerical calculations and to determine the cor- 
rect form of Simmons' equations, eqs (2-72 and 2- 
77) were used to compute the tunneling current 
density as a function of electrode spacing and bias 
voltage. The results are given in table 3 and 
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Table 1. Tunneling Between Identical Spherical Electrodes; Current and Current Density Versus Electrode Spacing and Work Func- 
tion. 

Bias Voltage = 0.1 Volt 
Temperature = 300 Degrees Kelvin 
Fermi Energy == 5.5 Electron Volts 





— WF= 


= 4.0 eV— 


— WF = 


4.5 eV— 


— WF=S.O eV— 


— WF=S.S eV— 


Electrode 




Current 




Current 




Current 




Current 


Spacing 


Current 


Density 


Current 


Density 


Current 


Density 


Current 


Density 


(nm) 


(A) 


(A/m 2 ) 


(A) 


(A/m 2 ) 


(A) 


(A/m 2 ) 


(A) 


(A/m 2 ) 


.3 


.256+0! 


.924+13 


.155+01 


.627+13 


.955+00 


.424+13 


.601+00 


.287+13 


.4 


.420+00 


.168+13 


.210+00 


.918+12 


.110+00 


.515+12 


.601-01 


.297+12 


.5 


.610-01 


.253+12 


.262-01 


.117+12 


.120-01 


.567+11 


,576-02 


.287+11 


.6 


.852-02 


.357+11 


.319-02 


.143 + 1! 


.128-02 


.609+10 


.546-03 


.273 + 10 


.7 


.118-02 


.495+10 


.384-03 


.173 + 10 


.136-03 


.648 + 09 


.514-04 


.258+09 


.8 


.161-03 


.680+09 


.460-04 


.207+09 


.144-04 


.685+08 


.481-05 


.242+08 


.9 


.218-04 


.926+08 


.547-05 


.247+08 


.151-05 


.721+07 


.449-06 


.226 + 07 


1.0 


.295-05 


.126+08 


.649-06 


.294 + 07 


.158-06 


.756+06 


.418-07 


.211+06 


1.1 


.398-06 


.169+07 


.766-07 


.347+06 


.165-07 


.790+05 


.388-08 


.196 + 05 


1.2 


.534-07 


.228+06 


.902-08 


.410+05 


.171-08 


.823+04 


.360-09 


.181+04 


1.3 


.716-08 


.306+05 


.106-08 


.482+04 


.178-09 


.856+03 


.333-10 


.168+03 


1.4 


.956-09 


.409+04 


.124-09 


.566+03 


.185-10 


.889+02 


.307-11 


.155+02 


1.5 


.127-09 


.545+03 


.146-10 


.663+02 


.191-11 


.921+01 


.283-12 


.143+01 


1.6 


.170-10 


.727+02 


.170-11 


.775 + 01 


.198-12 


.954+00 


.261-13 


.132+00 


1.7 


.226-11 


.966+01 


.199-12 


.906+00 


.205-13 


.986-01 


.240-14 


.122-01 


l.S 


.299-12 


.128+01 


.232-13 


.106+00 


.211-14 


.102-01 


.221-15 


.112-02 


1.9 


.397-13 


.170+00 


.270-14 


.123-01 


.218-15 


.105-02 


.203-16 


.103-03 


2.0 


.526-14 


.226-01 


.314-15 


.143-02 


.225-16 


.108-03 


.187-17 


.946-05 



demonstrate that his theory gives values which are 
surprisingly close to those of the numerical calcula- 
tions 3 . 

The product D(E X )N(E X ) is known as the nor- 
mal-energy distribution [51], and is the distribution 
of tunneling electrons in energy associated with the 
component of momentum normal to the electrode 
faces. The normal-energy distributions were calcu- 
lated and plotted to give insight into the thermal 
distributions and barrier transmission aspects of the 
tunneling process. Interpretation of figures 8, 9, 
and 10 in terms of the Pauli exclusion principle and 
the energy diagram of figure 7 is difficult. For a 
total energy distribution, the graphs should have a 
sharp peak at the Fermi level of the negatively bi- 
ased electrode and have a width equal to the bias 
voltage with some small amount of thermal smear- 
ing at both edges. Plots of the total energy distribu- 
tion for tunneling in Al-GaSe-Au structures are 



3 The negative current density values for an electrode spacing 
of 0.3 nm result from the approximations used in Simmons' the- 
ory. In his theory the current density is expressed as the differ- 
ence between two terms (see the equations in the computer 
program listing in Appendix I) whose magnitudes become ap- 
proximately equal at electrode spacings of 0.3 nm and work 
functions of 5 eV. 



given by Kurtin et al. [23]. The utility of the nor- 
mal-energy distribution may be seen by recalling 
that N(E X ) for zero temperature is given by eq (2- 
64); N(E X ) at room temperature is only slightly dif- 
ferent with small spreading at E X =E F and 
rounding at E X —E F — eV. Combining this trun- 
cated triangular shape with the graphs of figures 8, 
9, and 10 enables one to readily see the exponential 
changes in the distribution due to changes in the 
tunneling probability as E F —E X and (E F —eV)—E x 
are changed. Other observations from figures 8, 9, 
and 10 are: 1) that the sources of non-linearity of 
the distributions for different voltages are apparent 
in the increasingly unequal heights as the electrode 
spacing increases, and 2) that the thermal smearing 
of about 25 mV at room temperature results in the 
distributions for 1 and 10 mV being almost identi- 
cal. 

The other figures are reasonably self-explanatory 
and, while no further discussion of them will be 
given here, they will be referred to in the analysis 
of experimental data. 

The computer programs used to obtain the data 
described in items 1 through 4 are given in Ap- 
pendix I, along with annotations which explain the 
major mathematical difficulties encountered while 
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Figure lla-Theoretical current-voltage characteristic; E F 
(shown as FE) and WF are in eV; D is in pm. 
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Figure lid-Theoretical current-voltage characteristic; Ff 
(shown as FE) and WF are in eV; D is in pm. Repeat of figure 
1 la. with IFF of 4.0 eV. 
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Figure lib-Theoretical current-voltage characteristic; Fp 
(shown as FE) and WFare in eV; D is in pm. Repeat of figure 
11a with WFof 5.0 eV. 
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(shown as FE) and WF are in e V; D is in pm. 



_H CURRENT US. VOLTAGE T-388. E^E.S UF-4.S S-a3«a 




CURRENT US. VOLTAGE T-388. FE-S.S UF-4.5 D-E889 



— i — r — i — ■ — i — r — i ' — i — i — r— i — — i — i — i — i— 

e.i *. a a.a e.4 e.s 8.6 

UOLTASE llfOLTB* 




9.9 



"n — : — — [— 
8,8 8.1 8.2 8.3 8.4 8.5 0.6 

UOLTACE W0LTS* 

Figure llc-Theoretical current-voltage characteristic; E F 
(shown as FE) and WF are in eV; D is in pm. Repeat of figure Figure 12b-ThearetLcal current- voltage characteristic. Repeat 

11a with WF of 4,5 eV. of figure 12a with an electrode spacing of 2000 pm. 



191 



Journal of Research of the National Bureau of Standards 



_, g CURRENT US. UOLTACE 7-393. E F -&.5 UF-4.S D-E5 




CURRENT VS. DISTANCE T-390. E-5.5 U-.ie 



~T 1 1 ■ — 

e,a 0.3 s.< 

UCLTASE tUQLTSI 



Figure 12c-Theoretical current-voltage characteristic. Repeat 
of figure 12a with an electrode spacing of 2500 pm. 
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Figure 14-Theoretical log current-spacing characteristics for 
four work functions; bias voltage of 0.1 volt. Work functions 
are in units of eV. 
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Figure 12d-Theoretical current-voltage characteristic. Repeat 
of figure 12a with an electrode spacing of 3000 pm. 
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Figure 15-Theoretical current-voltage characteristics; work 
function of 4.5 eV, electrode spacings of 2.0, 2. 1, and 2.2 nm. 
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Figure ]3-Theoretical log current-spacing characteristics for 
four work functions; bias voltage of 0.01 volt. Work functions Figure lfr-Theoreiical log current density-electrode spacing 

are in units of eV, characteristics; work functions of 4.0, 4.5, 5.0, and 5.5 eV. 
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Figure ll-Theoretical current density-Mohage characteristics.; 
work function of 4.5 eV and electrode spacings of 2,0, 2.1, 
and 2.2 nm. 



bringing the programs up to a correctly operating 
state. There are, however, several additional points 
about the physics involved which should be de- 
scribed here. To calculate all the constants in- 
volved, the following values of the fundamental 
constants were used: 

m = mass of the electron = 9.11 XlO -3 ' kg 

e = charge on. the electron = 1.60 X 1CT 19 C 

fc B =Boltzmann's constant= 1.38 X \0~ 2i 
J/K=86.2 yeV/K 

h = Planck's constant =6. 625 X 10 _3+ J s. 

With these values, 



Table 2. Tunneling Between Identical Electrodes: Current Density Versus Electrode Spacing and Voltage (Current Density in Units 
of A/m 2 ). 



Temperature = 300 Degrees Kelvin 
Work Function = 5.0 Electron Volts 
Fermi Energy = 5.5 Electron Volts 













Voltage (Volte) 












Electrode 


























Spacing 


























(run) 


.OS 


.10 


.15 


.20 


.25 


.30 


.35 


.40 


.45 


.50 


.55 


,60 


.3 


.212+13 


.424+13 


.636+13 


,849+13 


.106+14 


.128+14 


.149+14 


,171 + 14 


,192+14 


.214+14 


.236+14 


,258+14 


.4 


.257 + 12 


,515+12 


.773 + 12 


.103+13 


.129 + 13 


,155 + 13 


.182+13 


.208 + 13 


.235 + 13 


.262+13 


.289+13 


.317+13 


.5 


.283 + 11 


.567 -rU 


.852+11 


.114+12 


.143 + 12 


.172+12 


.201 + 12 


.231 + 12 


.260+12 


.291 + 12 


.321 + 12 


,353 + 12 


.6 


,304+10 


.609+10 


.916+10 


.122+11 


.154 + 11 


,1854-11 


.217+11 


.249+11 


.282+11 


.315-11 


.349+11 


.383 + 11 


.7 


.324 + 09 


.648+09 


.975 +09 


.130+10 


.164+10 


.197+10 


.232 + 10 


.266-,- 10 


.302 + 10 


.338+10 


.375 + 10 


.413—10 


.8 


,342+08 


.685 + 08 


.103+09 


.138+09 


.173+09 


.209+09 


.246 +D9 


.283 + 09 


.322+09 


.361 + 09 


.402+09 


.444+09 


.9 


.360+07 


.721+07 


.109 + 08 


.14S+08 


.183+08 


.221 + 08 


.260+08 


.300+08 


.342+08 


.384+08 


.428+08 


.474+08 


1.0 


.377 + 06 


.756+06 


.114+07 


.153 + 07 


,192+07 


.233 + 07 


.274+07 


.317 + 07 


.362+07 


.408 + 07 


.456+07 


,506+07 


1.1 


.394+05 


.790+05 


,119+06 


.160+06 


.201+06 


.244+06 


.288+06 


.334+06 


.382+06 


,432+06 


,484+06 


.539+06 


1.2 


.410+04 


.823+04 


.124+05 


.167+05 


.211+05 


.256+05 


.303+05 


.351+05 


.403+05 


.456+05 


.513+05 


.573+05 


1.3 


.426+03 


.856+03 


.129+04 


.174 + 04 


.220-04 


,267+04 


.317+04 


.369+04 


.424+04 


.482 + 04 


.544+04 


.610+04 


1.4 


.442+02 


.889+02 


.134+03 


.181+03 


.229+03 


.279+03 


.332-^03 


.387 + 03 


.446-03 


.50B + 03 


.576-^03 


.648 + 03 


J. 5 


.458+01 


.921+01 


.139+02 


.188+02 


.238+02 


.291+02 


.347+02 


.406 +02 


.469+02 


.536+02 


.609+02 


.688+02 


1.6 


.474+00 


.954+00 


.144+01 


.195+01 


.248+01 


.303+01 


.362+01 


.425+01 


.492+01 


.565+01 


.644+01 


.731+01 


1.7 


.490-01 


,986-01 


.149 + 00 


.202+00 


.257+00 


.315+00 


.378+00 


.444+00 


.517 + 00 


.596+00 


.682+00 


.777+00 


1.8 


.506-02 


.102-01 


.155-01 


-209-01 


.267-01 


.328-01 


,394-01 


.465-01 


.543-01 


.628-01 


.721-01 


,825-01 


1.9 


.522-03 


.105-02 


.160-02 


.217-02 


.277-02 


.341-02 


.411-02 


.486-02 


,570-02 


,661-02 


.763-02 


.877-02 


2.0 


.538-04 


.108-03 


.165-03 


.224-03 


.287-03 


.354-03 


.428-03 


.509-03 


.598-03 


.697-03 


.808-03 


.933-03 


2.1 


.554-05 


.112-04 


.170-04 


.232—04 


.297-04 


.368-04 


,446-04 


,532-04 


,628-04 


.735-04 


.856-04 


.993-04 


2.2 


.570-06 


.115-05 


.175-05 


.239-05 


.308-05 


.382-05 


.465-05 


.556-05 


.659-05 


.775-05 


.906-05 


.106-04 


2.3 


.586-07 


.118-06 


,131-06 


.247-06 


.319-06 


.397-06 


.484-06 


.582-06 


.692-06 


.317-06 


,960-06 


.112-05 


2.4 


.603 -0B 


.122-07 


.186-07 


.255-07 


.330-07 


.412-07 


.504-07 


.608-07 


,726-07 


.862-07 


.102-06 


.120-06 


2.5 


.619-09 


.125-08 


.192-08 


.263-08 


.341-08 


.428-08 


.525-08 


.636-08 


,763-08 


,909-08 


.108-07 


.128-07 


2.6 


.636- 10 


.129-09 


.198-09 


,272-09 


.353-09 


.444-09 


.547-09 


.665-09 


,802-09 


.960-09 


.115-08 


.136-08 


2.7 


.654-11 


.133-10 


.204- 10 


.281-10 


.366-10 


.461-10 


.570-10 


.696-10 


.842-10 


.101-09 


.122-09 


.145-09 


2.8 


.671 - 12 


.136-11 


.210-11 


,289-11 


.378-11 


.479-11 


.594-11 


.728-11 


.886-11 


,107-10 


.129-10 


.155-10 


2.9 


.689-11 


.140-12 


.216- 12 


,299-12 


.391-12 


.497-12 


.619-12 


.762-12 


,931-12 


.113-11 


.137-11 


,166-11 


3.0 


,707-14 


,144-13 


.222-13 


,308-13 


.405-13 


.516-13 


,646-13 


.798- 13 


.980-13 


.120-12 


.146-12 


.177-12 
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Table 3. Tunneling Between Identical Electrodes; Current Density Versus Electrode Spacing and Voltage. Calculated from Simmons 
Theory (Corrected). (Current Density in Units of A/m 2 ). 

Temperature = 300 Degrees Kelvin 
Work Function = 5.0 Electron Volts 

Voltage (Volts) 



Spacing 






















<nm) 


,05 


.10 


.15 


.20 


.25 


.30 


.35 


.40 


AS 


.50 


.3 - 


-.704+12 


-.142+13 


-.216+13 


-.292+ 13 


-.371 + 13 


-.452+13 


-.538+13 


-.627+13 


—721 + 13 


-.821 + 13 


.4 


.101 + 12 


.204+12 


.308+12 


.413 + 12 


.520+12 


.628 + 12 


.738+12 


.849+ 12 


.962+12 


.10B+13 


.5 


.182+11 


.367+11 


.556+11 


.749+11 


.946+11 


.115 + 12 


.135 + 12 


.156+12 


.178 + 12 


.200 + 12 


.6 


.236 + 10 


.478 + 10 


.725+10 


.978+10 


.124+11 


.150+11 


.178+11 


.206+11 


.235+11 


.265 + 11 


.7 


.279 + 09 


.564+09 


.856+09 


.116+10 


.146+10 


.178+10 


.211 + 10 


.245 + 10 


.280+10 


.316+10 


.& 


.314 + 08 


.636+08 


.966+03 


.131-1-09 


.166 + 09 


.202 + 09 


.239+09 


.278+09 


.319 + 09 


.361+09 


.9 


.345 + 07 


.698+07 


.106+03 


.[44-r08 


.183 + 08 


.223+08 


.265+08 


.308+08 


.354+08 


.401+08 


1.0 


.372 + 06 


.755+06 


.115+07 


.156-07 


.198 + 07 


.242 + 07 


.288+07 


.336+07 


.386+07 


.439 + 07 


1.1 


.398 + 05 


,,807+05 


.123+06 


.167+06 


.212 + 06 


.260+06 


.310+06 


.362+06 


.417+06 


.476 + 06 


1.2 


.421+04 


.855+04 


.130+05 


.177 + 05 


.226 + 05 


.277 + 05 


.331+05 


.387 + 05 


.447+05 


.511+05 


1.3 


.444+03 


.901+03 


.138+04 


.187 + 04 


.239 + 04 


.293 + 04 


.351+04 


.412+04 


,477+04 


.547 + 04 


1.4 


.465+02 


.946+02 


.144+03 


.197 + 03 


.251 + 03 


.309 + 03 


.371+03 


.437+03 


.507+03 


.583 + 03 


1.5 


.486+01 


.988+01 


.151+02 


.206+02 


.264+02 


.325+02 


.391+02 


.461+02 


.537+02 


.619+02 


1.6 


.506+00 


.103+01 


.158 + 01 


.215 + 01 


.276+01 


.341+01 


.411+01 


.486+01 


.568+01 


.656+01 


1.7 


.526-01 


.107+00 


.164+00 


.224 + 00 


.288-00 


.357 + 00 


.431+00 


.511+00 


.599+00 


.695+00 


1.8 


.545-02 


.111-01 


.170-01 


.233-01 


.300-01 


.373-01 


.451-01 


.537-01 


.631-01 


.735-01 


1.9 


.565-03 


.115-02 


.177-02 


.242-02 


.313-02 


.389-02 


.472-02 


.563-02 


.664-02 


.776-02 


2.0 


.584-04 


.119-03 


.183-03 


.251-03 


.325-03 


.405-03 


.493-03 


.590-03 


.698-03 


.819-03 


2.1 


.603-05 


.123-04 


.189-04 


.260-04 


.337-04 


.422-04 


.515-04 


.618-04 


.734-04 


,865-04 


2.2 


.621-06 


.127-05 


.196-05 


.270-05 


.350-05 


.439-05 


.537-05 


.647-05 


.771-05 


.912-05 


2.3 


.640-07 


.131-06 


.202-06 


.279-06 


.363-06 


.456-06 


.560-06 


.677-06 


.810-06 


.962-06 


2.4 


.659-08 


.135-07 


.208-07 


.288-07 


.376-07 


.474-07 


.584-07 


.708-07 


.851-07 


.101-06 


2.5 


.678-09 


.139-08 


.215-08 


.298-08 


.389-08 


.492-08 


.608-08 


,741-08 


.894-08 


.107-07 


2.6 


.697-10 


.H3-D9 


.221-09 


.307-09 


.403-09 


.511-09 


.634-09 


.775-09 


.939-09 


.113-08 


2.7 


.716-11 


.147-10 


.228-10 


.317-10 


,417-10 


.530-10 


.660-10 


.810-10 


.986-10 


.119-09 


2.8 


.735-12 


.151-11 


.235-1] 


.327-11 


.431-11 


.550-11 


.687-11 


.847-11 


.104-10 


.126-10 


2.9 


.755-13 


.155-12 


.242-12 


.338-12 


.446-12 


.571-12 


.716-12 


.886-12 


.109-11 


.133-11 


3.0 


.774-14 


.159-13 


.249-13 


.348-13 


.461-13 


.592-13 


.746-13 


.927-13 


.114-12 


.140-12 



N{E X )= 1.3585 X I0 ,0 r In [ ] A/eV K m 2 



(2-78) 



The values of a = 10.25(eV)- |/2 (nm)~' and Ll5e 2 
ln2/(S7reo)= 0.2867 (eV nm) were also obtained 
with these values for the fundamental constants. 

Due to the nature of N(E X ), which was ex- 
plained earlier, eq (2-61 a) was decomposed into 
three parts: 

£" F -c-eV 

J\ =A TksT&zV/ksT) / D{E x )dE x ], 

(2-79a) 

MFDl/FDYtydE*, 
jF ' £ - eV (2-79b) 



/•E F +Wk B T 



with JszJi+Ji+J}. c was chosen as 40kgT and 
ETOP is the top of the barrier potential function. 
Neglecting small changes in the maximum barrier 
height with applied bias voltages, ETOP is given 
by 



ETOP= [0.575 e 2 ln 2/(8ire<,)]5" 



(2-80) 



(2-79c) 



J-3 was not included in the calculations unless 
ETOP was less than (E F +80k B T). 
The form for D{E X ) used in the calculations was 

i>(^)=(l+exp[2j[ L K0c)dx])- 1 , (2-81) 

K(x) is defined in eq (2-6 le). This form may be 
approximately obtained from eq (2-6 le) by factor- 
ing 4Q 2 from the second term and performing a 
series expansion. The result is £>(!£,)=( 1/2 4- 2 ) -1 . 
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Equation (2-81) has been obtained by Miller and 
Good [52] by analytic continuation and was shown 
to hold for all energies, including energies greater 
than the barrier height; one determines complex 
turning points in this instance. The Miller and 
Good form of the WKB transmission enabled the 
calculations to proceed smoothly up to the top of 
the barrier, where x l =Xr and D(E X )=0.5. The 
value for the exact transmission coefficient of a 
particle with energy at the top of a rectangular 
shaped barrier is 0,2. Since one would expect the 
image potential rounded barrier to be less reflec- 
tive, the 0.5 value is at least intuitively correct. To 
avoid the complexities of determining complex 
turning points, a simple linear interpolation be- 
tween 0.5 and 1 was used from ETOP to the vac- 
uum level. 

Current Calculations for Spherical Electrodes* 

Current flow was obtained from the current den- 
sity calculations by assuming that in the almost pla- 
nar electrode structure, the current flow was 
parallel to the axis of the two spheres. The incre- 
mental tunneling areas are then rings around this 
axis of radius r with spacing greater than the clos- 
est spacing between the electrodes. The current 
flow in this approximation is 



-J 



J(x)dA, 



(2-82) 



J(x) is the current density for an electrode spacing 
x, Xo is the on-axis electrode spacing and x m is the 
spacing between points off-axis at which J(x) has 
dropped to a negligible value. For r much smaller 
than the sphere radius R, 



'The high lateral resolution implied by the imaging of the Si 
(111) surface obtained by G. Binnig, H. Roarer, Ch. Gerber, 
and E. Weibel (Phys. Rev. Lett. 50 (1983) 120) has stimulated 
much theoretical investigation, of the effect of electrode shape 
on the tunneling current and effective tunneling area. The 
strong exponential dependence of the tunneling current on elec- 
trode spacing tends to confine the tunneling current to very 
small regions of one of the electrodes if it is sharpened. Some of 
the investigations of the shape of the electrodes on the tunneling 
current and tunneling area are: N. M. Miskovsky, S. J. Shep- 
herd, P. H. Cutler, A. A. Lucas, and T. E. Sullivan, Effect of 
Geometry and Multiple Image Interactions on Tunneling and 
I-V Characteristics of MVM Point-Contact Junctions, App. 
Phys. Lett. 37 (1980) 189; J. Tersoff and D. R. Hamann, Theory 
and Application for the Scanning Tunneling Microscope, Phys. 
Rev. Lett. SO (1983) 1998; N. Garcia, C. Ocal, and F. Flores, 
Model Theory from Scanning Tunneling Microscopy: Applica- 
tion to Au (110X1X2), Phys. Rev. Lett. SO (1983) 2002; T. E. 
Feuchtwang, P. H. Cutler, N. M. Miskovsky, Theory of Vac- 
uum Tunneling Microscopy, Phys. Lett. 99A (1983) 67; A. 
Baratoff, Theory of Scanning Tunneling Microscopy - Methods 
and Approximations, Physica 127B (1984) 143. 



r(jc>=[J£(jc-jeo)] 1/I . (2-83) 

With dA =%Tirdr, eq (2-82) is therefore 

I^ttR f "j(x)dx. (2-84) 

The 10" 11 factor employed in the programs in Ap- 
pendix I results from the use of x in nm and R in 
mm. 

3. Experimental Apparatus and Method 

The rule-of-thumb expression for the tunneling 
current /, at small voltage, V 



I = Cexp(-aSB h2 )A V, 



(3-D 



and its ramifications are the major factors which 
determine all aspects of the apparatus design, con- 
struction, and use of apparatus for metal-vacuum- 
metal tunneling. C is a constant, A is the device 
tunneling area, a is 10.25 (eV) _1/2 (mil)" 1 , S is the 
electrode spacing, and B is the mean barrier height. 
Some insight into the experimental problems which 
had to be addressed may be obtained by observing 
that according to eq (3-1) a change in electrode 
spacing of 0.1 nm, about one-half an atom diame- 
ter, will generate a factor of 10 change in the tun- 
neling current, assuming a barrier height of 5 eV. 
In a metal-vacuum-metal tunneling device, an elec- 
trode spacing stability of one pm is therefore re- 
quired to achieve a noise level of 2% or less in the 
tunneling current for a fixed voltage. This corre- 
sponds to a signal- to -noise ratio of only 34 dB. 

Mechanical and Thermal Requirements 

The electrode spacing stability is in turn deter- 
mined by the mechanical and thermal properties of 
the electrodes and the assembly into which they 
are placed. As representative values, let the effec- 
tive length, L, for thermally generated changes in 
electrode spacing be 5 mm and the expansion coef- 
ficient, k T , of the materials involved be 10 ppm/K. 
Then, the required temperature stability, dT, to 
achieve a 1 pm stability, dL , in electrode spacing is 



dT = dL {L At) " ' = 20 /nK 



(3-2) 



The laboratory in which the experiment was to be 
performed had a temperature stability of ±0,5 K 
over a 24-hour period. However, when local 
sources of heat such as that of the body heat of the 
experimenter are taken into account, the local vari- 
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ations in environmental temperature are more rep- 
resented by ± 1 K. An approximate attenuation of 
I0 _i of the environmental temperature variations 
must, therefore, be provided by a temperature con- 
trol system. 

Requirements for mechanical stability were de- 
termined by measuring the amplitude and fre- 
quency composition of building vibrations at the 
planned laboratory location of the experimental ap- 
paratus. The transducer used for these measure- 
ments was a Texas Instruments S-36 very low 
frequency seismic detector whose sensitivity was 
1.97 V/cm/s. The transducer's response was ap- 
proximately 10% greater at its resonance Fre- 
quency of 2 Hz and dropped by about 15% at 1 Hz. 
Figure 18 shows the velocity amplitude spectrum 
of the building vibrations in the National Bureau of 
Standards Reverberation Chamber. Measurements 
were performed here because it was one of the qui- 
etest potential locations for conducting the vacuum 
tunneling experiment. In addition, it was the only 
laboratory which allowed convenient access to the 
Fourier transform electronics necessary for obtain- 
ing real-time vibration spectra. Vibration spectra 
from this location were compared with other possi- 
ble locations for the experiment with the use of 
oscillographs. These oscillographs showed that the 
vibrations present on a vibration table with a reso- 
nant frequency of 1 Hz and located in a basement 
laboratory were similar in amplitude and frequency 
content to those measured in the well-isolated NBS 
Reverberation Chamber. The NBS Reverberation 
Chamber is a room-within-a-room structure with 
an inside room volume of about 400 ra 3 . The inner 
room is supported on springs to produce a room- 
spring resonance frequency of 3.5 Hz. 

Depending on the frequency, the vibration am- 
plitudes which must be accounted for, even with a 



vibration isolation table, are therefore between 4 
and 25 nm. The amplitudes are deduced from fig- 
ure 8 with the 25 nm amplitude corresponding to 
frequencies of about 3.5 Hz, and the 4 nm to a fre- 
quency of about 35 Hz. Figure 18 only shows the 
vertical component of vibration displacements; 
measured horizontal components were similar in 
amplitude, but with major contributions existing 
only at frequencies near 35 Hz and below 10 Hz. 
The mechanical design of the electrode holder 
must, therefore, attenuate these vibrations by about 
10~ 4 to produce the 1 pm stability. 

A second factor which must be accounted for in 
designing for mechanical stability is acoustical 
pressure variations in the air around the experimen- 
tal apparatus. No measurements of the magnitude 
and frequency content of the acoustical pressure 
variations were made since reliable estimates could 
be obtained from the literature. Hunter [53] gives a 
graph of the noise amplitude as a function of fre- 
quency for an average room. For a mean noise 
level of 40 dB referenced to 2 X 10 _s N/m 3 , the am- 
plitude at 20 Hz was 40 dB and dropped almost 
linearly to -5 dB at 5000 Hz. Within the labora- 
tory used for this experiment, there were no dis- 
tinctive noise sources, so this description of the 
acoustical environment was taken as representa- 
tive. 



Choosing an Electrode Material and Preparing the 
Electrodes 

Variations in the mean barrier height also pro- 
duce large changes in the tunneling current. For 
small changes in B, dB, eq (3-1) may be used to 
obtain 



dl/l=~ 1/2 aS B 1/2 dB 



(3-3) 



Figure 18-Measured velocity am- 
plitude spectrum of building vi- 
brations in the National Bureau 
af Standards Reverberation 
Chamber. 
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With 5 = 1 nm, B=4 eV, and dB=0.1 eV, dl/I 
— —0.26, B must be stable to 0.008 eV to achieve a 
34 dB signal to noise ratio. 

Changes in work function and thus barrier 
height are produced by the adsorption and 
chemisorption of any kind of foreign atom on a 
metal surface [27,32]. The change in the work func- 
tion WF, WFC, is a local microscopic effect, but in 
terms of a spatially averaged change, it can for- 
mally be related to an effective dipole moment p 
per adsorbed particle via the equation [32,54] 



WFC=4pN, 



(3-4) 



where N is the number of adsorbed atoms or 
molecules per unit area. Adsorption is viewed as 
changing the [V(+ oa)-V(-ao)] term of the WF 
described in section 2. Changes in WF range from 
+0.08 eV for a 1/3 monolayer of selenium on 
nickel to —3.1 eV for a monolayer of cesium on 
tungsten (100) [32,54]. Larger increases in the WF 
are also found, but usually no greater than 1 eV. 
WFC's due to more common adsorbates such as 
oxygen, nitrogen, and hydrogen are typically about 
1 eV, and may be of either sign depending on the 
substrate-adsorbate combination. 

This demand for stable clean electrode surfaces 
must be carefully balanced with the opposing need 
for high temperature stability. There are two major 
reasons for the conflict between these two require- 
ments. First, to achieve atomically clean metal sur- 
faces, such that the work function would not be 
changed by any contaminants, some form of in-situ 
vacuum cleaning process must be employed. The 
most common and reliable methods for in-situ 
cleaning all require large energy inputs either to 
simply evaporate or decompose all surface and 
near surface contaminants, e.g., high temperature 
flashing of W or high temperature cycling in oxi- 
dizing and reducing atmospheres, e.g., Ni in 2 and 
H 2 . Second, the time for a contaminating layer to 
form on a clean surface is comparable with esti- 
mated times for an electrode assembly to thermally 
equilibrate to the stabilities demanded by the vac- 
uum tunneling experiment. 

The time, T, in which an initially clean surface 
will become covered with a monolayer of adsor- 
bate when the surface is exposed to a gas at pres- 
sure P Torr is (see Appendix III) 



Tzs2xl0 6 /P seconds. 



(3-5) 



Equation (3-5) is only a rough estimate, but it is 
useful in determining the vacuum conditions neces- 
sary for this particular application. A minimum 



time for a system that has been raised to 1300 K to 
obtain a 10 uiC stability is about 24 hours or 
8.6 X 10 5 seconds, Thus, maintaining a clean surface 
over this period of time would require that the par- 
tial pressures of all possible adsorbates be less than 
10" n Torr. Vacuum pressures at this level are at 
the limits of modern technology and are achieved 
with only cryopumping or by very careful use of 
vacuum materials. 

As a compromise between the demands of tem- 
perature stability and the requirement for clean sur- 
faces, gold was chosen as the electrode material 
because none of the common adsorbates, including 
molecular oxygen, have been found to chemisorb 
at room temperature and it could be conveniently 
formed into electrodes with clean, smooth surfaces 
prior to installation into the temperature control 
and micropositioning assembly. 

The approach taken was to fabricate the gold 
electrodes in one UHV system, to store and trans- 
port the fabricated electrodes in desiccated con- 
tainers, and finally to assemble the electrodes into 
an electrode assembly in the laboratory atmo- 
sphere. As discussed in the next paragraph, this 
procedure certainly does not result in an ideally 
clean surface. However, after the electrodes are 
placed in a UHV environment, the surfaces should 
be altered by at most a monolayer of hydrocarbon 
and require no further processing before the vac- 
uum tunneling experiment could reliably be per- 
formed at a stable temperature. 

The most definite study of gold's resistance to 
adsorbing molecular species was performed by 
Chesters and Somorjai [55]. They found that of the 
species most likely to be encountered in a labora- 
tory environment, oxides, sulphides, hydrocarbons, 
nitrides, and water vapor, the only ones which 
were observed to adsorb on a gold (111) crystal 
face or a stepped crystal face were heavy hydro- 
carbons. Their investigations were with the crys- 
tal's temperature less than 500 °C and partial 
pressures less than 10 _s Torr. Naphthalene ad- 
sorbed on both types of gold surfaces and the ad- 
sorption behavior indicated dissociative 
chemisorption. Surface steps did not play a critical 
role in the adsorption properties of gold, in con- 
trast to that of the similar chemical material, plat- 
inum. Earlier work by Trapnell [56] using 
polycrystalline gold films showed that reversible 
molecular adsorption of CO, C2H4, C 2 H 2 , and other 
light hydrocarbons did occur with low heats of ad- 
sorption (0.4 eV/molecule) for high pressures. 
These low heats of adsorption allow the adsorption 
process to be reversible and Trapnell [56] shows 
that at low pressures and room temperature, the 
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light hydrocarbons desorb with no measurable, i.e., 
less than 2% coverage at pressures less than 10~* 
Torr. 

Work function changes, WFC, in gold due to 
adsorption have received very little study, in part 
because gold is the most common reference elec- 
trode for measurements of the work function [32]. 
In particular, no work to determine the WFC's due 
to the adsorption of heavy hydrocarbons was 
found in the literature. However, WFC's in plat- 
inum have been studied for a wide variety of 
molecular species. While the type and behavior of 
molecular species which adsorb on gold surfaces is 
markedly different from that on platinum surfaces, 
the JVFC's observed in the studies of adsorption on 
platinum at least set an upper bound on the WFC's 
to be expected for adsorption on gold. The values 
of the WFC's of the platinum do reflect the effec- 
tive dipole moments of the adsorbed molecules on 
a dominantly s -electron like surface and should 
give a valid estimate of the dipole moment of the 
molecule when it is adsorbed on gold. The major 



difference between platinum and gold from the 
standpoint of chemical properties is that platinum 
has an unfilled 5d electronic subshell. For plat- 
inum, the WFC's were measured by Gland and So- 
morjai [57] for a wide range of heavy 
hydrocarbons. WFC's were found to be between 
— 1.0 and —2.0 eV depending on the type of adsor- 
bate and the degree of coverage of the platinum 
surface. 

The gold electrodes were formed by melting the 
end of a 1.25 mm gold wire in UHV with electron 
bombardment. The apparatus used for forming the 
electrodes is shown and described in figure 19. 
Gold wire with 99.99% purity was used. Follow- 
ing degreasing with a trichloroethylene, acetone, 
ethanol sequence, the wires were etched in boiling 
aqua-regia, i.e., 3.5 parts HC1 and 1 part HN0 3 ; 
both acids were fresh ACS grade chemicals. After 
rinsing the wires thoroughly in running tap water, 
they were ultrasonically cleaned in de-ionized, fil- 
tered water and finally blown dry with hot air. The 
wires were then installed in the electron bombard- 
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Figure 19-Apparatus for forming 
the gold electrodes. 
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ment apparatus and vacuum pumping was begun. 

The vacuum system used for the melting opera- 
tion utilized a sublimation-ion pump combination 
which maintained pressures of 2 to 4X10" 8 Ton- 
while the gold was molten. With the 1,4 kV indi- 
cated in figure 19, only about 15 mA of bombard- 
ment current was required to melt the end of a gold 
wire. 

The stainless steel mass-thermal resistance is an 
essential component of the melting apparatus. This 
component combined with the nichrome wire pro- 
vides sufficient thermal resistance to allow an equi- 
librium between input and dissipated power to be 
obtained by convenient manual adjustment of the 
filament heating current. Without these compo- 
nents, the equilibrium is difficult to obtain; the 
nonequilibrium result in this case is that the molten 
region of the wire moves up the length of the wire, 
forms a sphere and drops off the wire. Aftef sev- 
eral of these occurrences, the configuration shown 
in figure 19 was developed to avoid this problem. 
With this design, a molten sphere could be main- 
tained on the end of the wire with the support of 
the liquid's surface tension for 10 to 20 minutes. 
These long periods in the molten state in vacuum 
provided a means to evaporate many contaminants 
from the surface and were sufficient for the surface 
tension forces to determine the microscopic surface 
topography. 

After melting, the electrode surfaces were in- 
spected with an optical microscope at a magnifica- 
tion of 300X. The process from etching in aqua 
regia to melting was then repeated until no surface 
contamination was visible. For the pair of elec- 
trodes used in the experiment, the process was re- 
peated four times. The resulting surfaces were very 
smooth with no visible topography or contami- 
nants. The electrodes were finally examined with a 
scanning electron microscope (SEM) after the ex- 
periment was completed. The SEM micrographs of 
the electrode surfaces are shown in figures 20 and 
21. The surfaces were sufficiently smooth that in 
both the electron and optical microscopes good fo- 
cus was difficult to obtain. 

The anticipated electrode surfaces for this exper- 
iment were, therefore, clean gold surfaces altered 
by at most a monolayer of adsorbed heavy hydro- 
carbons to produce work functions which could 
range from gold's work function, 5.0 eV, down to a 
minimum of 3.0 eV. 

Temperature Control System 

In order to provide the required stability of ± 1 
pm for the interelectrode spacing, the design goal 










; j#. 







Figure 20- Scanning-electron- microscope photomicrograph of 
gold electrode surface; view normal to surface. 
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Figure 21-Scanning-electron-microscope photomicrograph of 
gold electrode; surface projection. 



for the temperature control system was to produce 
a temperature stability of ±20 jiK., Other workers 
[58,59] have shown that by careful design, temper- 
ature stabilities at this level can be achieved with 
the use of a two-stage control system. This ap- 
proach was therefore adopted. The first and outer 
stage was designed to operate at a temperature of 5 
K above the laboratory temperature with a stability 
of ± I mK. The second and inner stage was de- 
signed to operate at a temperature 0.5 K above the 
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outer stage's temperature with the specified design 
goal stability. 

The apparatus used for the combined purposes of 
temperature control, vibration isolation, and acous- 
tical isolation is shown and described in figure 22. 
Vibration and acoustical isolation aspects of the ap- 
paratus are described in a later section. Several fea- 
tures of this apparatus should be emphasized; 1) 
Bifilar windings were used for both control heaters 
to insure that stray magnetic fields in the region of 
the electrode assembly would be minimal. 2) Radi- 
ation baffling was provided at the vacuum pump- 
ing part of the inner heater enclosure to minimize 
the radiant coupling between the outer heater's 
fluctuations and the electrode assembly. The baffle 
design insures that any radiation entering the 
pumping port must be scattered at least five times 
from the surface of the inner heater enclosure. 3) 
All joints in the polyurethane insulation and open- 
ings for the various vacuum parts were taped with 



building construction duct tape to minimize con- 
vection losses. 4) Electrical leads passing from the 
laboratory through the first stage and from the first 
stage through the second stage to the electrode as- 
sembly were in both cases solidly heat-sunk to the 
enclosure of each respective stage. For the outer 
stage, thermal contact for each electrical lead was 
obtained by soldering one point of the lead to a 
copper post which was mounted to the outer en- 
closure with the same high thermal conductivity 
epoxy used to mount the heater wires; see figure 
22. Thermal contact for leads passing through the 
inner stage enclosure was obtained by using copper 
feed-throughs in the base of the enclosure, item 3 in 
figure 22, which were electrically insulated with 
beryllium oxide washes. Total contact area for 
each lead was 3 cm 2 . 

A block diagram of the temperature control sys- 
tem employed for both control stages is shown in 
figure 23. The basic design approach of both stages 



Figure 22-Apparatus for vibration isolation, acous- 
tical isolation, and temperature control, consist- 
ing of the following components: I) Electrode 
holder and micropositioning assembly. These are 
described in figure 25. 2) Inner stage temperature 
control enclosure (6061 Aluminum) and heater. 
The heater wires shown are only schematic, ac- 
tual wires were of #30 nichrome in a bifilar con- 
figuration which almost completely covered the 
enclosure's cylindrical surface. The radiation baf- 
fle is the cross-hatched piece located immediately 
above the electrode assembly. 3) Aluminum (6061 
alloy) base for inner enclosure 4) Alumina stand- 
offs. Three were used to provide thermal isola- 
tion between the two temperature control stages. 
5) Outer stage temperature control enclosure, 
heater and vacuum chamber walls. The cylinder 
wall thickness was 3.2 mm, top and bottom thick- 
nesses were 1.9 cm. Heater wires were #31 cop- 
per spaced approximately as indicated in the 
drawing. An epoxy coating approximately 3.0 
mm thick with a thermal conduction coefficient 
of 0.01 (cal cm)/(sec cm : *C) [60], was used to 
bond the heater wires to all surfaces of the outer 
enclosure and to provide acoustical damping. 6) g 
cm thick polyurethane foam. This foam provided 
both thermal and acoustical isolation of the sys- 
tem enclosure and 7) Epoxy standoffs, three, 
were used to provide thermal insulation of outer 
enclosure and vacuum system from the vibration 
isolation table. These standoffs were bonded to 
the vibration isolation table top and permitted 
rigid coupling of the entire apparatus to the isola- 
tion table top. And 8) Granite isolation table. The 
table was 61 cmx25 cm and had an approximate 
mass of 409 kg. It was supported with four 
damped steel springs from a simple aluminum 
frame. The resonant frequency of the system was 
less than 1 Hz for vertical motion. 
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T (environment) 

Figure 23-Block diagram of temperature control electronics. 

was to use the output of a bridge to drive a high 
gain amplifier and compensation network to ac- 
count for the lags of the detectors and stage- 
masses, then to close the control loop by 
incorporating a thermistor detector as one element 
of the bridge. A full discussion of the electronics of 
both control stages is given in Appendix V. There 
it is shown that one of the most important parame- 
ters of the control loop is the time response of the 
thermistors with respect to changes in power input 
to the heaters. For this reason, the thermistors 
were mounted in a manner to allow the closest pos- 
sible coupling between the heaters and the thermis- 
tors. 

The principal difference between the two stages 
of temperature control was that a dc bridge and a 
sensor amplifier were used for the outer stage 
while an ac excited bridge and lock-in amplifier 
were used for the inner stage. The compensation 



network, power amplifier and heaters were tai- 
lored for the requirements of each stage. Values for 
the heat capacities, thermal resistances, and masses 
of each stage are given in Appendix V. 

Both controllers were designed to allow conve- 
nient and accurate changes in operating tempera- 
ture. The outer stage's setpoint could be changed 
in 2 mK increments with the use of a switched re- 
sistor network. Setpoints for the inner stages were 
changed with a seven-digit ratio transformer which 
permitted 10 u-K increment changes. 

Performance of the outer control stages was 
checked by using the inner stage detector system as 
a temperature monitor. When this stage's compen- 
sation network was adjusted as described in Ap- 
pendix V, temperature deviations from the 
setpoint, measured over a 10-hour period, were no 
greater than ±0.5 mK. The only check on the in- 
ner control stage's performance was to monitor the 
stability of the interelectrode spacing. The results 
of this check are given in section 4 where it is 
shown that no changes in electrode spacing due to 
temperature variations were observed. A major in- 
dication of the sensitivity of the inner stage con- 
troller was that a setpoint change of 30 juK would 
cause an observable change in the output of this 
stage's sensor amplifier. The response of both 
stages to setpoint changes is described in Appendix 
V. 



Vacuum System for the Tunneling Experiment 

A schematic diagram of the vacuum system used 
for the tunneling experiment is given in figure 24. 
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The vacuum chamber and pumping system were 
constructed from 304 stainless steel with standard 
UHV copper gasket flanges employed for all con- 
nections between system components. The system 
was pumped from atmospheric pressure down to 
the operating pressure using only the three pumps 
shown in figure 24 in order to avoid exposing the 
gold electrodes to any unnecessary hydrocarbon 
gases. Base pressure of the vacuum chamber as de- 
termined by a crosscheck of the ion-pump current 
and the Bayard-Alpert gauge was 5xlO~ 9 to 
1 X 10~ 8 Torr after several weeks of pumping. The 



Bayard-Alpert gauge was not operated during the 
critical temperature control parts of the experi- 
ment. 



Electrode Holder and Micropositioning Assembly 

The final design of this critical part of the vac- 
uum tunneling experiment is shown in figure 25. A 
photograph of the electrode assembly, along with 
the temperature control heaters and vacuum sys- 
tem, is given in figure 26. The assembly as shown 




Figure 25-Electrode holder and micropositioning assembly. (1) Housing for electrode holders. This housing is basically an aluminum 
(6061 alloy) block 8.9 cmX6.35 cmx6.35 cm with a 3.17 cm hole bored and reamed to accept the electrode holders. The diameters 
of the electrode holders were machined to match the hole in the housing to a clearance no greater than 0.01 mm. (2) Gold electrodes. 
In the experiment, both gold wires were inserted into their respective holders as far as the spherical ends would allow. (3) OFHC 
copper mount for the differential-screw mechanism and main body of left electrode holder. (4) 0.025 mm thick polyimide film [63] 
to provide electrical insulation between the electrode holders and the housing. (5) Diaphragm flexures. The thickness of both flexures 
was 0.38 mm. (6) Threads to allow prestressing of the piezoelectric element-diaphragm flexure system. (7) Prestressing nut. (8) 
Hardened steel washer to prevent motion transfer between prestressing nut and the piezoelectric element. (9) Electrical connections 
to gold electrodes. (10) Clamp-holder for left electrode, OFHC copper. (11) Threads to allow mounting of the differential-screw 
mechanism. (12) Hardened steel end caps to distribute stress uniformly over piezoelectric element's surface. (13) Clamp-holder for 
right electrode. (14) Support for outer diaphragm flexure, OFHC copper. (15) PZ displacement element. 
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in figure 26 provides a system which is compatible 
with ultra-high vacuum, has mechanical rigidity to 
maintain electrode spacings to stabilities of approx- 
imately 1 pm, and permits the electrodes to be as- 
sembled into a tunneling structure on a laboratory 
scale. In addition, the assembly is capable of pro- 
ducing interelectrode spacings of 1 to 2 nm and has 
the facility to change the electrode spacing in a 
rapid and convenient manner. 

The basic approach of this design is that a preci- 
sion-differential screw mechanism is inserted at 
point 11 in figure 25 to position the left electrode 
clamp, item 10, which contains one gold electrode 
until the electrode spacing is less than the distance 
which can be spanned by the piezoelectric (PZ) 
element. At this point, all components of the as- 
sembly are clamped rigidly in place, the differen- 
tial-screw mechanism is removed, and finarfine 
tuning of the electrode spacing is achieved by tem- 
perature adjustment and piezoelectrically gener- 
ated displacements of the second electrode holder. 

With reference to figure 25, the left electrode 
holder is composed of the electrode clamp, and a 
differential-screw-mount. Both were made of 
OFHC copper to obtain minimum thermal resis- 



tance between the electrode and the holder hous- 
ing and minimum electrical resistance between the 
electrode and its electrical connection. The elec- 
trode clamp was essentially two half-cylinders con- 
nected together by a 0.25 mm X 0.5 mm flexure 
along the length of the clamp and located next to 
the hole for the electrode. A separate screw is pro- 
vided for clamping the electrode into place. The 
mount for the differential-screw had a similar form. 
All other clamping is performed with the elec- 
trode-holder housing which is cut only in the 
shaded area of the drawing. 

Each vernier division of the differential-screw- 
mechanism corresponded to an output displace- 
ment of about 25 nm. This resolution was such that 
approximately one-quarter turn of the differential- 
screw's input covered the PZ element's 500 nm 
range. 

The right electrode holder consisted of an elec- 
trode clamp mounted to a double diaphragm-flex- 
ure support to insure axial motion of the electrode, 
a PZ element, and a mechanism to prestress the 
driving element-flexure combination so that micro- 
scopic motion losses would be minimized. End 
caps for the PZ element and the prestressing mech- 
anism were made from hardened steel. The di- 




Figure 26-Electrode assembly, temperature control heaters, and vacuum system. 
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aphragm-flexures and piezoelectric element were 
prestressed with a force of approximately 890 N 
(200 pounds). All other parts of the electrode 
holder, except the 304 stainless steel fasteners, were 
machined from OFHC copper. The PZ element 
was electrically insulated from the rest of the 
holder with four sapphire spheres. 

The PZ element consisted of six ceramic discs, 
each 1.27 cm in diameter by 0.25 cm thick. The 
ceramic was a barium-titanate-zirconate composi- 
tion which has very low hysteresis, linear charac- 
teristics [61]. The relevant expansion coefficient for 
the material is d 3} = 225 pm/V for fields along the 
poled direction of the ceramic; here perpendicular 
to the disc faces. The individual disks and end caps 
were bonded into one element with the use of a 
low viscosity, electrically conductive epoxy. The 
epoxy's electrical resistivity was less than 0.002 
ohm-cm [51]. After baking in air at 100 °C for one 
hour, the outgassing rate was sufficiently low that 
the 10~ l0 Torr base pressure of a bell jar system 
with a 50 1/s ion pump was unaffected by the pres- 
ence of the epoxy bonded piezoelectric element. 

Displacement characteristics of the completely 
assembled micropositioning system were measured 
before and after the tunneling experiments de- 
scribed in section 4. Before the experiments, the 
displacements were measured with a precise but 
uncalibrated mechanical gauge to be (.56 ±.06) 
nm/V. After the experiment, the displacements 
were measured with a Hewlett-Packard 5526A 
laser interferometer measurement system in combi- 
nation with a 10565B remote interferometer and 
retroreflector. This measurement system enabled 
displacements to be measured with an accuracy of 



±25 nm. For both expansion and contraction of 
the piezoelectric element, the displacements were 
measured as (0.47 ±.03) nm/V. Thus, taking into 
account the uncalibrated state of the mechanical 
gauge, the displacement characteristics of the mi- 
cropositioning system were stable throughout the 
experiment. 

The need for 1 pm electrode spacing stability 
combined with the 470 pm/V sensitivity of the PZ 
element requires that the high voltage power sup- 
ply used to drive the PZ element have a very low 
voltage noise level, less than 2 mV. The power sup- 
ply used for this purpose had a noise level of less 
than 75 uV peak-to-peak and the capability of be- 
ing switched in 1 mV increments from — 1 100 volts 
to +1100 volts [64]. To enable continuous voltage 
sweeps over ±10 volts, a ramp generator was 
added in series with this power supply. The circuit 
diagram of the ramp generator and a filter circuit 
to suppress switching transients are shown in figure 
27. 

Two other experimental techniques for holding 
and positioning the tunneling electrodes were thor- 
oughly explored before finally adopting the elec- 
trode holder and micropositioning assembly design 
just described. The technique first explored was to 
fabricate the tunneling structure as a microelec- 
tronic device with the use of photolithographic 
methods. This approach was not adopted because 
the technology to which the author had access 
could not meet the dimensional requirements of the 
tunneling device. However, as a byproduct of this 
work, the approach may result in an electronic 
device with a very high transconductance; values 
up to 1 mho appear feasible. A patent application 
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on the concept is being submitted to the National 
Bureau of Standards Patent Review Committee. 
The concept and a prototype device are described 
in Appendix IV. 5 

Lever amplification of the piezoelectric ele- 
ment's length changes was also explored as a possi- 
ble means to bridge the gap between laboratory 
scale distances and the microscopic dimensions re- 
quired by the tunneling experiment. While this 
technique proved to have excessive vibration levels 
for the tunneling experiment, it has characteristics 
which are very useful for scanning-electron and 
optical microscope stages [65]. 

Vibration and Acoustical Isolation 

In the previous metal-vacuum-metal tunneling 
experiment [12] similar to this work and in all ear- 
lier attempts made in connection with this work, 
vibrations of the tunneling device structure at the 1 
nm level have limited a successful execution of the 
experiment. Every effort was made, therefore, to 
carefully examine the mechanical design of each 
part of the apparatus to insure that all known 
sources of vibrational energy were attenuated to 
levels which would not result in detrimental 
changes in the electrode spacing. The vibrational 
properties of what was thought to be the most 
likely contributors to changes in electrode spacing 
will now be described in terms of one-dimensional 
models. 

Starting with the electrode holder and micropo- 
sitioning assembly, the most likely contributor is 
vibrational motion of the right electrode holder as 
the entire assembly is vibrated along the assembly's 
axis. This motion could be produced either from 
acoustical pressure variations on the bottom of the 
vacuum chamber or from building vibrations driv- 
ing the isolation table top and, thus, the entire vac- 
uum chamber. As a model of this motion, let a 
mass, m , be connected by a spring with force con- 
stant k, and a viscous damper characterized by the 
constant C, to a vibrating wall whose displacement 
is given by x a =A Q sin at. Then the equation of 
motion of in in terms of x, measured relative to the 
wall is 



x =R*[(l-Rif + 4yiRi]-w Xa 



(3-7) 



tn (x+x a )-i- Cx+kx=0. 



(3-6) 



The steady state solution to eq (3-6) with the as- 
sumed form of x a is 



S A variation on the concept proposed here has recently been 
developed by J. Morelaud and J, W, Elkin, Electron Tunneling 
Eiperiments Using Nb-Sn 'Break' Junctions, J. Appl. Phys. 58 
(1985) 3888. They use a bending beam to produce and vary the 
spacing between electrodes, 



R =td/<o Q , 6i =(k/m) 1/2 , and V= l/2C(mk)~ 1/2 . If at 
is much less than a>o, then 



=R*x 



(3-8) 



Thus, if oq is the lowest resonant mode fre- 
quency of the micropositioning assembly for axial 
vibrations and the full amplitude of the isolation 
table top's horizontal motion is transmitted along 
the axis to the electrode holder housing, / (=<»o/ 
2tt) must be greater than 6 kHz. This follows with 
the use of an w of (27r)30 radians/s, an x„ of 40 nm, 
and a requirement that x be less than 1 pm. The 
only mechanical element with a resonant frequency 
near this value is the diaphragm flexures, which 
were designed to have a resonant frequency of 10 
kHz, when unstressed. When stressed, the fre- 
quency would be even higher and motion would 
lien be determined almost entirely by the pre- 
stressing mechanism. 

For motions normal to the axis of the assembly 
all resonant frequencies are much higher than 10 
kHz. The relatively large gold wires were used to 
avoid low resonant frequencies in those directions. 
For example, even with a very high modulus of 
elasticity material such as tungsten, the transverse 
vibrational resonant frequency of a cantilever only 
2 mm long and .025 mm in diameter is 4 kHz [66]. 

To minimize the affects of any acoustically gen- 
erated vibrations, the electrode holder and micro- 
positioning assembly are only connected 
mechanically to the laboratory environment by the 
bottom plate of the vacuum chamber. This plate is 
heavily damped both by the epoxy layer for attach- 
ing the heater wires and by the polyurethane foam. 

The four lowest eigenfrequencies of the bottom 
plate are approximately: 884 Hz, 1.85 kHz, and 3.0 
kHz, and 3.5 kHz [66]. These neglect any mass 
loaded decreases in the frequencies due to the mass 
of the tunneling assembly, which is supported by 
the plate and any increases in the frequencies due 
to the epoxy damping. The lowest of these frequen- 
cies again is the most likely source of transmitting 
acoustical pressure variations to the electrode as- 
sembly. As a safety factor for design, a 40 dB white 
noise spectrum was assumed to exist up to one kHz, 
i.e., pressure variations, P a , with amplitudes of 
2X 10~ 3 N/m 2 . For frequencies much less than the 
eigenfrequencies, the vibrational amplitudes, Y, of 
the plate at the center of the plate would be [55] 



r=3P u (l-s 3 )«V(4(e) ! 



(3-9) 



205 



Journal of Research of the National Bureau of Standards 



s=Poisson's ratio=0.33, a = radius of the plate 
= .127 m, t= thickness of the plate =1.9 cm, and Q 
= modulus of elasticity of the plate = 1. 86 X10 10 
N/m 2 . With the indicated P a , Y is only 8.7 pm. For 
frequencies on resonance, the amplitudes of vibra- 
tion should be no larger than 10 times this value 
with the epoxy and polyurethane damping [67]. 
Since the assumed pressure variations were larger 
than expected and did not account for any attenua- 
tion due to the polyurethane foam, no further stiff- 
ening of the plate was performed. 



Tunneling Electronics 

The electronics used to measure the conduction 
properties of the tunneling devices are shown in 
figures 28 and 29. Requirements of the electronics 
were that it permit one to test devices whose 
impedances could range from \0~ 2 to 10 1! ohms 
with dc test voltages of typically 10 mV. The de- 
sign shown in figure 28 utilized an operational am- 
plifier in a follower configuration to insure that, 
other than voltage drops in the current leads, the 



Figure 28-Block diagram of tun- 
neling electronics. 
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Figure 29-Schematics for compo- 
nents of the tunneling electron- 
ics. (A) Ramp generator for 
sweeping bias voltage. (B) Man- 
ual sweep and bias supply. (C) 
Logarithmic amplifier. 
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output impedance of the voltage source would be 
less than 0.005 ohm. The operational amplifier was 
a 356H integrated circuit device. The four-wire 
lead system was carried through to the connection 
points on the electrode holders shown in figure 25. 

With all systems of the tunneling experiment in 
operation, i.e., vacuum system, temperature con- 
trol system, high voltage power supply for the mi- 
cropositioning assembly, and the tunneling 
electronics, the total system was carefully de- 
bugged for ground loop problems and 60 Hz 
pickup. The noise current in the tunneling elec- 
tronics after this process was about 10" " A. This 
value was measured with a 10 s ohm test resistor 
substituted for the tunneling device and obtained 
from both the /- F characteristics of the resistor and 
the output of the logarithmic current to voltage 
converter. The major source of this noise was the 
operational amplifier and could have been reduced 
with a more elementary resistive voltage genera- 
tor. However, the level was adequate for this ex- 
periment, and the alternative was a greatly 
increased source resistance. 

Both dc voltmeters employed with the tunneling 
electronics had floating-guarded inputs with dc 
common-mode rejection ratios of greater than 140 
dB. The voltmeter used to measure the bias voltage 
applied to the tunneling device had a gain of 10 
with readout resolution to 10 u,V, The voltmeter 
for measuring the tunneling current, in the I-V 
configuration of the electronics, had variable gains 
up to 10* with a meter noise of about 150 nV p-to- 

P- 

The logarithmic current-to-voltage amplifier 
output was trimmed according to the manufac- 
turer's specifications. After warmup, the offset was 
less than 100 u-V. Specifications for the amplifier 
modules were that the output voltage would be a 
linear function of the logarithm of the input current 
over a range from 10~ 9 to 10 -3 A to within 1%. 
Measurements with test resistors showed that, in 
addition, the nonlinearity was only 5% for the pos- 
itive polarity module and less than 2% for the neg- 
ative polarity module for currents of 10 -10 A. To 
reduce the acoustical noise level of the chart 
recorder when the tunneling device formed an 
open circuit, a lowpass filter with a cutoff fre- 
quency of 32 Hz (see figure 29c) was installed. 

In operation, the electronics were changed from 
the Log / cv — S configuration shown in figure 28 to 
the configuration for I-V measurements by first 
shorting the input of the Log module with the 1 
ohm resistor, disconnecting the Log module and 
the piezoelectric sweep generator and finally con- 
necting the dc voltmeter for current detection. 



Changing back to the Log I QV — S configuration was 
performed in a similar manner. 



4. Experimental Results and Analysis 

Taking into account the extreme delicacy of the 
interelectrode spacing adjustments, the apparatus 
described in section 3 has proved to be very reli- 
able and to give reproducible results. The data 
which will be presented were obtained with the 
same essential characteristics in two separate ex- 
periments with data being obtained over a period 
of several weeks in both instances. 

After several unsuccessful attempts at getting the 
apparatus to perform as anticipated, a standard se- 
quence of operations was evolved which permitted 
the experiment to be carried out in a consistent 
manner. Following formation of the electrodes and 
their installation into the holder and microposition- 
ing assembly, as described in section 3, the assem- 
bly was placed on a temperature-controlled plate 
of copper. The temperature of the plate was main- 
tained at (300 ±.05) K which was the planned oper- 
ating temperature for the experiment The 
electrode spacing was then observed with a low 
power stereo-microscope and adjusted with the 
differential-screw mechanism until the electrodes 
appeared as shown in figure 30. 

Concurrent with this adjustment the current 
flow through the electrodes for an applied voltage 
of 10 mV was monitored with the tunneling elec- 
tronics. Final adjustment of the left electrode in the 
sense of figure 25, was such that with the PZ ele- 
ment at the maximum extension possible in air the 
electrodes would just make contact. Since the ex- 
tension of the PZ element in air was only 0.7 of 
that possible in vacuum, slight misadjustments 
could be corrected after installation in the tempera- 
ture-controlled vacuum system. Misadjustments 
out of range of the PZ element, either in contrac- 
tion or expansion, could be corrected with small, 
±0.5 K, changes in the operating temperature. Use 
of the temperature-controlled copper plate was re- 
quired to overcome the expansion and contraction 
generated by unknown changes in temperature re- 
sulting from body heat during the adjustment pro- 
cess. 

All parts of the electrode holder and microposi- 
tioning assembly were firmly clamped in place at 
this time. The assembly was then allowed to equili- 
brate with "hands-off." (Plastic gloves over cotton 
gloves were actually used at all times during the 
assembly to provide thermal insulation and to in- 
sure the cleanliness required for UHV.) After con- 
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Figure 30-Gold electrodes mounted in electron assembly. All 
the apparent structure and cuts in the electrodes are images 
formed by the mirror like properties of the electrode surfaces. 



firming that an open and a short could be produced 
by contracting and expanding the PZ element, the 
electrode assembly was installed in the vacuum sys- 
tem and pumped to a pressure of approximately 
1 X 10 8 Torr. This usually required several weeks 
with the small ion pump which was used on the 
vacuum system. 

Determining the PZ element voltage and operat- 
ing temperature at which the electrodes were posi- 
tioned to within 1 to 2 nm of each other was 
accomplished by searching for the abrupt change 
in tunneling current as the electrode spacing was 
swept through the tunneling region. With the cir- 
cuit shown in figure 27 the sweep rate for a range 
of 5 nm was approximately 0. 1 nm/s. The search 
was greatly facilitated with the logarithmic current 
amplifier. In fact, all attempts at finding the re- 
quired PZ voltage and temperature settings with- 
out the use of a logarithmic amplifier were 
unsuccessful. After completing this procedure to 
locate the tunneling region the first time, succes- 
sive searches were not difficult. 



Data and Analysis 

The first sets of Log I cv —S data, while exhibiting 
the character of many orders of magnitude change 
in current over 1 to 2 nm change in electrode spac- 
ing, contained much structure akin to noise; see fig- 
ure 31. However, some of the structure may be due 
to the oscillations in the tunneling probability pre- 
dicted by numerical calculations based on an exact 
solution to Schrodinger's equation [45]. Regardless 
of the possibilities, the exploration of the phe- 
nomenon was prevented by an accidental electrical 
discharge as an attempt was being made to measure 
the constant-current characteristics which were 
measured in both previous vacuum tunneling ex- 
periments [12]. 

The series of measurements which was made 
subsequent to this accidental discharge of the con- 
stant-current power supply will be referred to as 
the first experiment. Representative Log I cv —S 
data obtained during this experiment are given in 
figures 32a and 32b. A comparison of figures 32 
and 3 1 shows that much of the structure present in 
figure 31 was probably due to fluctuating surface 
conditions. The larger scale oscillations in figure 3 1 
were not found in any subsequent Log I cv —S mea- 
surements even though they were consistently 
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Figure 31-One of the first measured log current-spacing charac- 
teristics. Curve A was measured as the PZ element was ex- 
panded then curve B as it contracted. The hysteresis is a 
combination of hysteresis in the PZ element and mechanical 
relaxation of the electrodes. Cycling gave reasonable repro- 
duction of both parts of the graph. 

present in all measurements prior to the accidental 
discharge. Exploration of the meaning of the oscil- 
lations and the conditions necessary for their obser- 
vation was left for future experiments. 6 

The overall shape and features of the data shown 
in figure 32 were reproduced about 50 times over a 
period of several weeks. Small changes in average 
slope and local details were observed as the elec- 
trodes were brought into contact during repeated 
closing and opening of the electrode spacing. All 
the Log I ev —S data had three characteristic re- 
gions. Expressed in terms of the tunneling current, 
the First of these regions was the initial transition 
into the current range of the Log amplifier. This 
region extended up to currents of approximately 
1CT 10 A and covered the non-linear part of the Log 
amplifier's response. The direct tunneling region 
usually extended from currents of 10 _ia A up to 
10~ 5 A for bias voltages of 10 mV. The final region 
occurred as the spacing of the electrodes was such 
that multiple portions of the electrode surface, 
other than the previously dominant part, began to 
make significant contributions to the total current 



^Subsequent investigations of this phenomenon have proved 
very interesting. R. S. Becker, J. A. Golovchenko, and B. S. 
Swartzentruber (Electron Interferometry at Crystal Surfaces, 
Phys. Rev. Lett. 55 (1935) 987) and G. Binnig, K. H. Frank, H. 
Fuchs, N. Garcia, B. Reihl, H. Rohrer, F. Salvan, and A. R. 
Williams (Tunneling Spectroscopy and Inverse Photoemission'. 
Image and Field States, Phys. Rev. Lett. 55 (19S5) 991) demon- 
strated oscillations in the tunneling voltage as a function of elec- 
trode spacing For a constant tunnel current Figure 31 is a graph 
of tunneling current as a function of electrode spacing for a 
constant bias voltage. The scale of the structure here is similar 
to that shown in the two references cited. 



and as the initially dominant portions made con- 
tact. 7 

A second attempt at obtaining constant-current 
measurements was made with a different power 
supply. In this attempt the transient discharge oc- 
curred again and was sufficient to microscopically 
weld the electrodes. Efforts to melt the bridge re- 
sulted in broken lead wires. The entire tunneling 
system was therefore disassembled, new lead wires 
installed, repositioned as described earlier, reassem- 
bled and all of the set-up operation repeated. Once 
the vacuum was at the 10 -8 Torr level, the elec- 
trodes were cleaned with controlled capacitive dis- 
charges. A representative Log I C \~S graph from 
this second experiment is given in figure 33. Five 
discharges of both polarities obtained from a 0.47 
ufd capacitor charged to 14 V were used for the 
cleaning and/or shaping of the electrode surfaces. 
The surface condition of the electrodes for the data 
of figure 33 is obviously different from that of fig- 
ure 32, but the average slopes in the direct tunnel- 
ing region of the two sets of data are approximately 
equal. The data in figure 33 were also reproduced a 
large number of times over several weeks. 

The electrode surfaces were examined in detail 
in a SEM after the completion of both experiments. 
No evidence of the discharges could be found. 

To confirm that the shape and magnitude of the 
Log I ov —S graphs were not being significantly af- 
fected by the time response of the Log amplifier, a 
very slow scan of the PZ element's voltage was 
performed. The data shown in figure 34 were ob- 
tained at a scan rate of 0.01 nm/s. One may con- 
clude from a comparison of figures 33 and 34 that 
the shape of the Log curves was not affected. 

A Log 7 CV versus time measurement was also per- 
formed to check for temperature and any other in- 
stabilities in the experiment. Results are shown in 
figure 35. None of the changes in current appears 
to be the result of temperature variations since the 
transitions occur so rapidly. Changes such as the 
ones shown in figure 35 could be due to either sud- 
den changes in surface conditions generated by the 
relatively large electric fields involved or to mo- 



'There have been two subsequent observations of the expo- 
nential variation of tunneling current as a function of electrode 
spacing for a constant bias voltage: G. Binnig, H. Rahrer, Ch. 
Gerber, and E. Weibel, Tunneling Through a Controllable Vac- 
uum. Gap, Appl. Phys Lett. 19 (1982) 178; and C. F. Quate, 
Low-Temperature Vacuum Tunneling Spectroscopy, Appl. 
Phys. Lett, 45 (1984) 1240. The capability for making such ob- 
servations is implicit in the work by Becker, et al. cited in foot- 
note (d) and in the work by Moreland, et al., cited in footnote 
(g), but no explicit graphs of current versus distance or equiva- 
lent are given. 
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Figure 32a-Measured log current- 
spacing characteristics. The 
sloping dashed line corresponds 
to a work function af approxi- 
mately 3.0 eV. The 100 mV 
curve was measured immedi- 
ately following the ID mV 
curve. Horizontal placement of 
the two curves was determined 
by the zeros of the PZ element's 
ramp generator. 
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Figure 33-Measured log current-spacing characteristics, bias 
voltage =10 mV. 
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Figure 32b-Measwed log current-spacing characteristics, bias 
voltage =10 mV. Curve A was measured as the PZ element 
expanded, then curve B. 

tion of the electrode surfaces resulting from lattice 
relaxation from the strains produced by the me- 
chanical stress of van der Waals forces and electro- 
static forces. The magnitude of these mechanical 
stresses is discussed in the next section. 

The work functions of the electrodes may be ob- 
tained directly from the Log J CT — S data given in 
figures 32 and 33. The slope of these curves in the 
direct tunneiing region may be related to the work 
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Figure 34-Measured log current-spacing characteristics; ex- 
panded spacing scale, bias voltage= 10 mV. 
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CHANGE IN CURRENT 
J CORRESPONDS TO ELECTRODE SPACING CHANGE OF 0,3 NM 




Figure 35-Log current vs. time for 
a fixed bias voltage of 10 mV. 
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function with the use of eq (2-70). By taking the 
ratio at two different spacings, one obtains 



QnI l ~lnJ 2 )/(.W l -W 2 )m~aB 



1/2 



where a = 10.25 (ev)- 1/2 (nm)" 1 , ^=barrier 
widths at the Fermi level, and B= average barrier 
height. The equation only serves to qualitatively 
illustrate the method relating the slope to the WF 
since B is a function of WF. The exact relation be- 
tween the slope of Log / cv — 5 data and the WF can 
be obtained from the calculations of section 2, s The 
sloping dashed line in figure 32 was drawn with the 
use of tables 1 and 2 and figure 13. Based on the 
discussion of the changes in WF with adsorbed hy- 
drocarbons, the WF value of 3.G eV is very reason- 
able. Note that this WF value was obtained without 
the use of any assumed property of the device or 
the electrode structure. The only input information 
was the piezoelectric element expansion coefficient 
and measured currents. This is the first time that 
such measurements have been performed. All pre- 
vious measurements employed separate tunneling 
devices with thickness distributed over the values 
allowed by the fabrication process. Determinations 
of the tunneling barrier thickness were also made 
from capacitance measurements and therefore 
were not independent of the device characteristics. 
Current-voltage measurements of the tunneling 
device were made during the second experiment, 
i.e., condition B, The results for an approximate 



l A. recent examination of many of the topics discussed in sec- 
tion 2 and of the validity of measuring work functions from the 
slope of (he In lev versus electrode spacing data has been per- 
formed by M. C. Payne and J. C. Inkson, Measurement of Work 
Functions by Tunneling and the Effect of the Image Potential, 
Surface Science 159 (1985) 485. In addition to the In J cv versus 
electrode spacing calculations given in section 2, they also use a 
formula derived by Simmons to compute the log of the tunnel- 
ing voltage as a function of the electrode spacing for a constant 
tunnel current density. 



electrode spacing of 1.2 nm are shown in figure 36. 
The approximate spacing of the S curve was ob- 
tained by determining the PZ element voltage at 
which touching occurred, assumed to be at the first 
major drop in slope of the Log I CV ~S curve, then 
contracting the PZ element the desired distance. 
The other curves were obtained by changing the 
applied PZ voltages by the amount to generate 0.1 
nm changes in electrode spacing, i.e., 0.2 V. 

Many attempts were made to obtain a continu- 
ous series of I-V graphs with electrode spacings 
varying from 2 nm to that when the electrodes 
made contact. In all attempts, the state of the elec- 
trodes did not remain sufficiently stable for the 
measurements to be performed. A series of mea- 
surements at two to three electrode spacings such 
as shown in figure 36 was obtained four times. 9 

Relating the I-V date to the theoretical calcula- 
tions of section 2 requires an estimate of the tunnel- 
ing area and the absolute electrode spacing. (The 
deduction of WF's from the Log 7 SV — S data only 
required knowledge of changes in electrode spac- 
ing.) Ideally a consistent approach to obtaining 
both the area and the absolute spacing would be 
first to deduce an absolute spacing with the use of 
eq (2-73) which relates the curvature of the I-V 



'The author is not aware of any other published results which 
show the variation of metal- vacuum-metal tunneling I-V curves 
with small changes in electrode spacing such as is shewn in 
figure 36. The data given by J. Moreland, S. Alexander, M. 
Cox, R, Sonnenfeld, and P. K. Hansma (Squeezable Electron 
Tunneling Junctions, AppL. Phys. Lett. 43 (1983) 387) for silver- 
vacuum-silver tunneling exhibit similar properties in that their 
curves have an ohmic region near zero bias with non-linear be- 
havior at larger bias voltages i.e., greater than about 0.2 V. 1-V 
data for structures with lead electrodes are also given which 
demonstrate behavior consistent with superconducting elec- 
trodes. /- V data is also given by S. A. Elrod, A. L. de Lozanne, 
and C F. Quate for small bias voltages applied to a W/vacuum/ 
NbjSn structuie, in the paper Low-temperature Vacuum Tun- 
neling Microscopy, Appl. Phys. Lett. 45 (1984) 1240. 
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Figure 36a-Measured current- 
voltage characteristics for three 
electrode spacings. S was ap- 
proximately 1.2 nm. The S, 
S — 0,1, and S— 0.3 nm curves 
were drawn manually from the 
A, B, and C curves of the "raw 
data" of figure 36b, respec- 
tively. The uncertainty in S is 
much larger than that of the dif- 
ference between the indicated 
spacings. The indicated noise 
level was a function of voltage, 
increasing with increasing 
voltage. The source of the noise 
was apparently changing sur- 
face condition. 



Figure 36b-Measured current- 
voltage characteristics for three 
electrode spacings. Curves A 
and B were measured with elec- 
trode spacings differing by 0.1 
nm; curve C differs in spacing 
by 0.2 nm from curve B. 
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characteristics to the barrier height and spacing. 
Since the barrier height can be obtained from the 
slope of the Log I m — S data, the absolute spacing 
would then be the only undetermined parameter. 
Once it is obtained, the current density calculations 
could then in combination with the corresponding 



Log Icv—S data and/or the I-V data be utilized 
tocalculate an average tunneling conduction area. 
Let G{V)=dl{V)/dV , then carrying out this 
procedure for the 5—0.1 nm curve of figure 36 
gives G(+0.5)/G(0) =21.42 and G(-0.5)/ 
G(0)= 12.23. The corresponding S parameters are 
Si =25. 13 B nm 2 and Si = 15.0 B nm 2 . With the B 
values calculated from the Log 7 CV — S data in fig- 
ure 32, both values are greater than 6 nm. Based on 
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the procedure outlined earlier for obtaining these 
data, S should have been Less than 2 nm. A possible 
explanation of this discrepancy in terms of mechan- 
ically generated strains is given in the next section. 

Since this method for obtaining an absolute spac- 
ing was not reliable, an alternative means to calcu- 
late the area and spacing was tried. The second 
approach was to check for correspondence be- 
tween the area calculated from the ratio of theoret- 
ical current and current density calculations of 
section 2 and the area estimated from the Log 
I cv — S data with the use of table 1 or eq (2-70). 
From the results given in table 1 a conduction area 
of approximately 2.4 X 10" ' 3 m 1 , with a slight varia- 
tion with work function and electrode spacing, is 
obtained. If a spacing parameter S of 1.2 nm, a mea- 
sured current of 10~* A , and a work function of 3.0 
eV are assumed, then one may estimate a tunneling 
conduction area of approximately 5X 10" 16 m 2 by 
extrapolating the values given in table 1. This area 
corresponds to a region of 20 nm radius or about 
l/20 ,h that of the radius expected for perfect 
spheres. Such an area is possible and could be pro- 
duced by microscopic projections from, the spheri- 
cal surfaces. 

However, a second possible hypothesis to ex- 
plain this difference between measured and pre- 
dicted currents is that the assumed point of contact 
occurred when S was not zero. Several processes 
discussed in section 2 would produce this effect. 
First, the transition from equal Fermi energies to 
equal bulk electrostatic potentials, even though a 
small effect for the almost identical electrodes of 
this equipment, could produce a change in slope of 
the LogJ s „—S curve that would be interpreted as 
contact. Second, adsorbed or chemisorbed 
molecules would produce a nonzero 5 at contact. 
Third, a change in slope at some value of S will be 
produced by the onset of constriction resistance [7] 
dominated current flow. For this conduction pro- 
cess the electrical resistance, R, between two large 
conductors connected by a small conduction area 
of radius, a, is given by R =p/(4a), where p is the 
conductors' resistivity. Gold's resistivity is 2.3 
XlQ -8 fl-m. Thus with the area resulting from the 
first hypothesis, the expected conduction resistance 
would be approximately 0.3 A. The measured resis- 
tance at "contact" was, however, approximately 
10* A, Constriction resistance for the conditions en- 
countered in this experiment, therefore, does not 
appear to be a current limiting factor. 

No consistent interpretation of the data given in 
figures 32, 33, and 36 has been found which gives 
an absolute electrode spacing and a tunneling con- 
duction area. Future experiments will be required 



to resolve the unknown relationships between the 
barrier parameter <S and the true physical spacing, 
and the difficulties of unknown conduction area. 

Somewhat as a last resort the I*V characteristics 
of the tunneling device were explored at large elec- 
trode spacings. Figure 37 shows the breakdown be- 
havior which was typically observed for voltages 
greater than about 1.0 V, even with relatively large 
electrode spacings. 

Effects of Mechanical Stress on Tunneling 
Characteristics 

There are two likely sources of attractive forces 
between the tunneling electrodes which could, by 
producing surface deflections, affect both the Log 
/ vs. S and the I vs. V characteristics. These are the 
van der Waals type forces between two micro- 
scopic bodies and the force arising from the elec- 
trostatic field produced by the applied bias voltage. 

The van der Waals type forces between two mi- 
croscopic bodies spaced a distance D apart in a 
vacuum may be calculated from a general micro- 
scopic theory of dispersion forces proposed by Lif- 
shitz [68,69]. When D is less than the smallest 
characteristic absorption wavelength of the mate- 
rial of the bodies, the force per unit area, S VDVr , 
between two planar bodies is given by 



Svdw=W(8 7T 2 Z> 3 ) 



(4-1) 



ficu is known as the Lifshitz-van der Waals interac- 
tion constant, By using optical data to deduce the 
dielectric functions for metals and the Lifshitz the- 
ory, Krupp, Schnabel, and Walter [70J computed 
fia for the case when both bodies were gold; fiu> 
was found to be approximately 1.9 X 10""' J. With 
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Figure 37-Correnr-TOftage characteristic for electrode Spacing 
of approximately 3.D nm. The breakdown behavior was typi- 
cal of all attempts to obtain I-V measurements for bias 
voltages greater than about 1.0 V. 
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D expressed in nm, S VD w between gold planes in 
vacuum is therefore given by 



( S VD w Au =2X 10 7 /D J N amVm* 



(4-2) 



The stress, T, produced by an electric field, E, 
perpendicular to a metal surface is [4] 



r= \/2<i<£\ 



(4-3) 



e =8.85xl0" 1! FD/m 1 . For a voltage V sad elec- 
trode spacing D (in nm), T is therefore given by 

7-=4.42x 10 6 {V/Df N nm7(V> m 2 ). (4-4) 

Lifshitz's theory is based on the existence of a fluc- 
tuating electromagnetic field in the region between 
the two bodies and the resultant component of the 
Maxwell stress tensor. Thus, both types of these 
forces may be viewed as being the static and dy- 
namic manifestation of the same force. 

A complete treatment of the strains generated in 
the gold electrodes by these stresses would be an 
extremely difficult problem. One would need to ac- 
count for the nonplanarity of the two electrodes in 
terms of its effect on the van der Waals forces and 
the effect of the electric field not being parallel to 
the axis of the electrodes. However, since the radii 
of the effective area is small compared to the 
sphere radius, a reasonable approximation is that all 
the forces act parallel to the axis of the two spheri- 
cal electrodes. Further assume that the equation 
derived by Timoshenko and Goodier [72] for the 
deflection, W(r,Q), at a point ir,Q) on a boundary 
plane of a semi-infinite solid, produced by a stress, 
q{r',&), over the area r'dr'd& holds for the sur- 
faces of the spherical electrodes. Timoshenko and 
Goodier's relation between W and q is 

W{r,Q)={{\ -p'ViTrY)] g(r',&) dr'd&, (4-5) 

p is Poisson's ratio for the solid and Y is the elastic 
modulus. Equation (4-5) is based on the assumption 
that the solid's elastic properties are isotropic. 

Let A) be the minimum spacing between the 
electrodes assuming no mechanical stresses were 
present. Then, for the small deviations from the 
axis of the electrodes, the spacing, D , at a distance 
r from the axis is given by D=D +r 2 /R; R is the 
radius of the electrodes. With the assumption about 
the force directions made earlier, the maximum de- 
flection will occur on the axis. Here eq (4-5) gives 

W =[2(1 -p 2 )/Y] J g(r')dr\ (4-6) 



Substituting the sum of the stresses described by 
eqs (4-2 and 4-4) for q into eq (4-6) gives 

W=a \ [bD{ r T l +CV 2 D(rT 2 \^\ (4-7) 

a =2(1 -p'-yY, b =2X 10 T N nmVm 2 , and 
C=4.42X 10 6 N nmW 2 m 2 ). 

Let y ^D (r')/D then dr = \/l (D^R ) yi 

O-ir^rfyand, 

J } trMt 
W>/Do)y- 3 

(y-ty-w + Clfiy-ity-iy^y. (4-8) 

Both parts of the integrand are standard forms. The 
integral evaluates to: 

W(D )=(a/2) (R/Dj) m [b(y - \f n /QI>&*) 

+ (3b /4D a + CV 2 )(y - l) v '/y 

+ tan- 1 (y-l) l '' 2 K«". (4-9) 

Let y^ go to infinity, then 

W(D Q )=(a/2) (R/DoY 1 rr/2[(3b/4D Q ) 

-fCK 2 ]. (4-10) 

Y for Gold is 8.3 X 10' N/m 1 , p =0.33 and R = 1.5 
X 10* nm; with these values, eq (4-10) reduces to 



W(£>o)=2.1xlO- 7 -Do-^ 

(3&/4A,+C^)nm. 



(4-11) 



For the voltages used to obtain the Log I„ vs. 5 
data, 10 and 100 mv, the dominant term is that 
from the van der Waals force. If D is 1 nm, then 
W=3.l nm; if D =2 nm, then W=0.56 nm. Taking 
into account the relatively gradual and controlled 
behavior of most of the Log I&, vs. S data, the 
"avalanche" magnitudes obtained from this calcu- 
lation are very puzzling. The history of the elec- 
trodes offers a possible explanation. During the 
process to obtain the first sets of data for Log I av 
vs. S, typified by figure 31, an abrupt transition 
such as one would expect from the predictions of 
eq (4-1 1) was observed. However, after the acci- 
dental and purposeful electrode cJeaning by capac- 
itive discharges, the abrupt or almost step function 
behavior was no longer observed, e.g., figure 32. If 
the capachive discharges produced localized re- 
gions of smaller curvature, the magnitude of the 
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deflections, from both type forces, would be 
greatly reduced, i.e., if it =10 3 nm, then W would 
be reduced by a factor of approximately 35 from 
that calculated from eq (4-11). 

The density of electromagnetic modes for the 
spherical surfaces would be slightly different from 
that of the flat planes used by Lifshitz. Other mi- 
croscopic geometrical deviations such as plateaus 
raised above a relatively smooth background 
would also affect the mode density. However, the 
geometry in both cases so closely approximates 
that of two planes over the areas which would con- 
tribute, i.e., about 10 jim in radius, that the changes 
in the mode densities should not affect the result 
given in eq (4-11). 

The effects predicted by the voltage dependent 
part support the previous proposal that the anoma- 
lously large curvature in the I vs. fdata for a fixed 
spacing was due to mechanical stresses. The differ- 
ence between the theoretically calculated value of 

I for a particular voltage and spacing, e.g., figures 

II and 12, and that measured, e.g., figure 36 would 
correspond approximately to a voltage produced 
deflection of 0.08 nm for voltages of 0.35 V. These 
values were obtained from the s— 0.3 nm curve of 
figure 36 as estimates of deviations from a linear 
I-V characteristic. The change in W between V = 
and V=0.35 V with A>= 1-5 nm as obtained from 
eq (4-11) is 0.06 nm. 10 

Conclusions 

Albeit crudely, an experiment has been per- 
formed which demonstrates quantum mechanical 
tunneling of electrons through a vacuum region of 
continuously variable thickness. Within experimen- 
tal error the results indicate that an image potential 
reduced barrier with the modifications suggested 
by Lang and Kohn gives a close approximation to 
the true barrier for metal-vacuum-metal tunneling. 
For the first time, the work function for the elec- 
trodes in a tunneling experiment has been deduced 
from experimental parameters independent of the 
tunneling device. 

The device which was developed for this experi- 
ment demonstrates that vacuum tunneling experi- 
ments are feasible and can reliably be performed. 



'"Further exploration and utilization of the effects of van der 
Waals forces has recently been initiated to understand the oper- 
ation, of the scanning tunneling microscope (J. H. Coombs and 
J. B. Pethica, IBM Europe Institute, STM Workshop, Oberlech, 
Austria, July 1-5, 19*5, IBM J. Res. and Dev. (to tie published)) 
and to enable the investigation of insulator surfaces (G. Binnig, 
C. F. Quate, and Ch. Gerber, The Atomic Force Microscope, 
Phys. Rev. Lett. 56 930 (March 1986). 



With this demonstration many new areas of re- 
search may be explored. A particularly direct ap- 
plication of the vacuum tunneling technique is that 
of locating the true position of a surface to a preci- 
sion of 0.1 nm, 20 to 100 times greater than 
presently available techniques. Without the limita- 
tion to electrode materials which form oxides that 
are good barriers, vacuum tunneling experiments 
will be a very versatile tool for performing inelastic 
electron tunneling spectroscopy. 

Some contributions toward accomplishing the 
objectives given in section 1 have therefore been 
made, but much remains to be done. The problems 
which have been discovered, e.g., strains generated 
by field generated stresses and the unknown rela- 
tionship between physical and theoretical spacing 
parameters, will need resolution before further 
work can be performed. However, the work here 
does give direction for solving the problems and 
once they are solved the full potential of vacuum 
tunneling experiments can be realized, 
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Appendices 



Appendix I. NUMERICAL CALCULATIONS OF TUNNELING CURRENT 



The computer programs, all in FORTRAN, 
which were used to calculate the data for the 
graphs and tables given in section 2, are contained 
in this appendix. The programs were executed on 
the National Bureau of Standards 1108A Comput- 
ing Facility. The main programs of these computa- 
tions are CURDEN, to calculate the tunneling 
current density; CURNT, to calculate the tunnel- 
ing current between spherical electrodes and a pro- 
gram to calculate the tunneling current density 
according to Shmnins' theory. 



Aside from the problems of roundoff and over- 
flow errors, the major difficulty which had to be 
overcome was long computational times. In an ear- 
lier version of these programs, the Newton-Raph- 
son method was used to find the classical turning 
points for the required energy values and a stan- 
dard iterative Simpson's 1/3 rule integration was 
used for all integrations. With these algorithms, the 
computational time for one current density value 
ranged from 20 to 30 seconds, with the longer 
times corresponding to larger electrode spacings, 
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The process of determining the turning points is 
deeply nested within the chain of integrals to cal- 
culate the current and current density. Its computa- 
tional time is therefore a major component of the 
overall computational time. 

By changing to the indicated tanh rule for inte- 
grating over K(x) and using an explicit solution for 
the roots of a cubic equation, the computational 
times were cut by a factor of 10 to 15 [73]. The 
iterative method of increasing the number of parti- 
tions until the difference between the result for n 
and 2n partitions was less than a specified error 
was also dropped. After determining what number 
was required for the needed accuracy, a fixed num- 
ber of 64 or 128 was used in all integrations. The 
particular number is specified in the programs in- 
volved. 

With the programs as indicated, the time to com- 
pute one current density value was two seconds; 



the time for one current value was 1 5 seconds. The 
resulting values, within the limitations of the physi- 
cal model, should be accurate to five and four sig- 
nificant digits, respectively, since doubling the 
number of partitions did not change the least signif- 
icant digit. 

The fundamental data unit of the 1 108A central 
processing unit is a 36-bit word. Average process- 
ing rate is about 900,000 instructions per second; 
0.75 microseconds for integer addition to 17.25 mi- 
croseconds for double precision floating-point divi- 
sion. Single precision floating-point numbers 
carried about 8 significant decimal digits (27 bits) in 
the range: 2xlO~ J9 to lXlO +ai . Double precision 
floating-point numbers carried about 18 significant 
decimal digits (60 bits) in the range: 3xl0 _3D9 to 
8X10 4307 . 




30C 
10 



PROGRAM SIMMON 



COMPUTE APFPOXIMATE VALUES OF C'JBPENT DENSITY 
ACCORCSIJS TO SIMMONS THEORY, SEE THE FCLLOWIWG REFERENCE 
BURSTEIN F.LIA5 *M5 LUNDQVIST STIG TLMajtJG PHENOMENA 
IN SOLIDS" CK 13 3V J G SIMMONS. 
EQUATION 11. p. 141 HAS BEOI CORRECTS} TO REAO AS FOLLOWS 

S2 - SM-e3Ka/|3WF*|10/KS!-2VJ*S1 

NOTE THAT FOR < - 1. VACUUM BAFFIER, W - S2 - SI REBUTS. 
VBI IS F.QN 21 AND JO IS EON IS IN ~HE REFERENCE CITED. 

REAL JOJT 
DIMEMSICK JTfll: 

W(S,V,WF] - S*[ 1 .-23.ra.'WF'S+ 1 0.-2."V"SB 
XBKS.V.WS > .S2-WF -tV1(2.-S)].'iS US2i 
I -<2.8S;lS2-5l>]"ALOG((S2 , (5-S1)MS1 '(S-S2DI 
DATA WFAC'5.0,-1025.6.JEIJ.' 
OPEtfeFIUE-'CUFCUT-;- 

DO 10 I - 3.3D 
S - I 

.Td> . sno. 

00 2D J - 2 11 

V - (J-1.K2D. 
SI - 3.,'WF 

S? . W(S.V.WF1*S1 
LM I - WSS*J WF) 

VI - VBKS.V.WF.S1.S21 
V2 - VI * V 

JO - CfVl , EXPlA-WV'SQRTIV1W(W1"2) 
JO - JO , C'|V2'EXPIA«W1>SQRT(V2)II'|W1"2I 
JTUI - JQ'l1.t«.E-9 , i[W"300.|»2)/Vll) 
CONTINUE 

wffirB's.sao) (jTtU.L-i.nj 

=CBMATf ! 2X FJ. 1. 1X.1DES.3I 

CCNTIHUE 

STOP 

EKD 



Listing of CURDEN 

C PROGRAM TUNNEL. CURRENT DENSITY (CURDEH] 

C 

C WRITTEN BY E. CLAYTOH T£AGL=E MARCH 1978 

C 

C SUBPBOQUknS CALLED: RS iBV3 , PHOD 



EX - EV)/KT| I 



this 13 a prqgrah to calculate the tunneling cubr3nt as a 
function of electrode spacing and sbubc voltage, the insul- 
ating barrier is assumed to be vacuum as modified by the image 
potential, the calculations are based on a free £lectroh model 
for the electrodes with a free electron mass hi both electrodes 
and in vacuum . a symmetrical structure hith identical 
electrodes is assumed. 

the basic tunneling equations assuming a wkb solution 

to the scp r0dingkr equation are! 

(c.3. duke, tunneling in solids, academic press, new toilk, 1969, 

PP. -49-53) 

CURRENT CENSITYIV, D, T,)= ISTEGKAHO TO INFIlDEX' D| EX, V, D, ) 
M2X, V, T)l 

R'HERE ; 

NIEX, V, TJ = (■UPDMEKTl/tl*"! * LNtXPOHTl 

WHEHE XPOBT = 1 + EX?[(FE-EX1/KT1/[1 +EIJtE > [ f FE 

4[PI)MEK/1I**3 - 1.335E10 AM PS/< E V, K } 

DIE>, V, D( = (i + EXP[-J(EX, V, D>])**-1 

WHERE 3(EX, M, Dl * I HI ['111BX, VJ TO H21EH, U)l 

aP[4[EI)(2W|*"0.5!/H]*lU[X)- EXJ*»0.5 

THE; BOTTOM OF THE C0NDUCTS3N BAUD OF THE "LEFT" 

ELECTO0D1 IS ASSUMED TO BE E •> . A VOLTAGE V IS APPLIED TO 

LOWER THE ELECTRON ENEUGIES OF THE "RIGHT" SLECTRODE, 

THUS TH2 TOTAL POTENTIAL ENERGY FUNCrtOtJ FOR THE BARRIER IS 

GIVES BY: 

U(X) = 0, X [jBSS THAN Oil EQUAL TO 

I![X) » 3ARRIEP H3IGHT - ELECTRIC POTENTIAL -IKAGE POTENTIAL 
< X < L. 

ilex) = - ev >: greater than or equal to d. 

THE BARRIER 15 THEREFORE GIVEN BY: 

U[X) = FE + WE - (EV/DJK - 

[3*<J) D ' LI 2 ' i.lS / ((8 PI * EPPO * X(D -X))** -1J 

If 11 19 III HH KHD V , FE, AHDWF AHE In VOLTS , 
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THIS REDUCES TO: 



U(X) = FE + WF - (V/D) 



0.2867 D/ (X[D-X) ) . 



THUS FOR CALCULATION ONE MUST FIRST DETERMINE FOR EACH VOLTAGE 
THE CLASSICAL TURNING POINTS XKEX, V) AND X2(EX, V) FOR EACH 
RELAVENT ENERGY BETWEEN AND FE + WF. 

THE TUNNELING PROBABILITY IS THEN CALCULATED AS A FUNCTION 

OF ENERGY, VOLTAGE, AND BARRIER, THICKNESS . AFTER 

CALCULATING THE SUPPLY FUNCTION "N" AS A FUNCTION 

OF ENERGY, VOLTAGE, AND TEMPERATURE, THE CURRENT DENSITY FLOW 

IS FINALLY CALCULATED AS A FUNCTION OF THE VOLTAGE 

BARRIER THICKNESS, AND TEMPERATURE. 

THE USER IS REQUESTED TO SUPPLY THE TEMPERATURE OF THE 
TUNNELING DEVICE, AND THE WORK FUNCTION AND FERMI ENERGY 
OF THE ELECTRODES. 

REAL NTGL1,NTGL2,NTGL3 
COMMON T,FE,WF,V,D,FK 
EXTERNAL PROD 
DATA NP /128/ 

NP2=NP/2 

GET TEMPERATURE 

WRITE (6,90) 
READ (5,80) T 



GET FERMI ENERGY 



WRITE (6,100) 
BEAD (5,80) FE 



Listing of DNSITY 



C 

10 

C 



FUNCTION DNSITY (Z) 
WRITTEN BY E. CLAYTON TEAGUE MARCH 1978 
SEE ANNOTATIONS IN CURRENT DENSITY PROGRAM. 
SUBPROGRAMS CALLED: RSIM13, PROD 

REAL NTGL1,NTGL2,NTGL3 

COMMON T,FE,WF,V,Z1,FK 
EXTERNAL PROD 
DATA NP /64/ 

NP2=NP/2 
Z1=Z 

BLZCTT=8.6E-5*T 
DLTAM0 . *BLZCTT+V 
ETOP=FE+WF-l. 1468/Z 

SET MAXIMUM ENERGY BARRIER HEIGHT AS WF + DLTA 

A=ETOP-WF-DLTA 

IF (A.LT.0.1 A-0 

B=FE-DLTA 

FK=1. 

NTGL1=-RSIM13(FR0D,A,B,NP2) 

A=FE-DLTA 

DIFF»ETOP-(FE+80.*BLZCTT) 

B=ETOP 

IF (DIFF.GT.0.1) B-FE+80.*BLZCTT 
FK=2. 

NTGL2»RSIM13(PROD,A, B,NP) 

IF (DIFF.GT.0.1) NTGL3-0. 
IF (DIFF.GT.0.1) GO TO 10 
A=ETOP 

fl=FE+80,*BLZCTT 
FK=3. 

NTGL3-RSIM13(PROD,A,B,NP2) 

DNS ITY=NTGL1+NTGL2+NTGL3 

RETURN 
END 



Listing of PROD 



C 
C 
C 

10 



FUNCTION PROD (EM) 
WRITTEN BY E. CLAYTON TEAGUE APRIL 1978 
SUBPROGRAMS CALLED: RNTOP, TUNNPB 

COMMON T,FE,WF,V,D,FK 

K=FK 

GO TO (10,20,30), K 

NOTE THAT T IN CONSTANT OF RNTOP WAS CANCELLED BY T IN EV/KT. 

PROD-(1.385E10)*(V/8.6E-5)*TUNNPB(EM> 
RETURN 



C 
20 

C 
30 
C 



PROD-RNTOP ( EM ) *TUNNPB ( EM ) 
RETURN 



PROD= RNTOP ( EMI * ( l + ( D/2 . 293 6 ) * ( EM-FE-WF ) ) 



RETURN 
END 



Listing of RNTOP 



C 
10 

C 

c 

20 
C 



FUNCTION RNTOP (EM) 

WRITTEN BY E. CLAYTON TEAGUE APRIL 1978 

DOUBLE PRECISION FERMIL, FERMIR, ALPHA, BETA, BLZCTT 
COMMON T,FE,WF,V 

BLZCTT-8.6D-5*DBLE(T) 
ALPHA=DEXPI DDLE(FE-EM)/BLZCTT) 
BETA=DEXP(DBLE(FE-EM-V)/BLZCTT) 

AVOID ROUNDOFF ERROR IN DOUBLE PRECISION BY USING APPROXIMATION 
THAT LN|(1+A)/(1+B)] - A-B WHEN A AND B ARE MUCH LESS THAN 1. 

IF (ALPHA. LT.1D-6) GO TO 10 

AVOID OVERFLOW ERRORS BY USING APPROXIMATION THAT 
LN[(1+EXP(A))/(1+EXP(A+B) )) » -B WHEN 
A IS MUCH. GREATER THAN 1. 

IF (BETA.GT.1D6) GO TO 20 

FERMIL-1.0D0+ALPHA 

FERMIR=1.0D0+BETA 

RNTOP-1.385Dl0*DBLE(T)*(DLOG(FERMIL)-DLOG(FERMIR)) 

RETURN 

RNTOP=1.385D10*DBLE(T)*( ALPHA-BETA) 

RETURN 

RNTOP-1.385D10*DBLE(T)*DBLE(V)/BLZCTT 

RETURN; END 



Listing of TUNNPB 

FUNCTION TUNNPB (EM) 
C 

C WRITTEN BY E. CLAYTON TEAGUE APRIL 1978 

C 

C SUBPROGRAMS CALLED: DCUBIC, RNIHF, RTPOT 



C 

10 



DOUBLE PRECISION A3,A2,A1 ,A0 , Z1R, Z1I, Z2R, Z2I, Z3R, Z3I 
EXTERNAL RTPOT 
COMMON T,FE,WP,V,D,FK 

CALCULATE X1(EM,V,D) AND X2(EM,V,D) TURNING POINTS. 

CUBIC COEFFICIENTS CALCULATED FROM POT (IN SUBPROGRAM RTPOT). 

A3=V/D 

A2— FE-WF+EM-V 

A1=(FE+WF-EM)*D 

A0=-.2867*D 

CALL DCUBIC ( A3 ,A2,A1 ,A0, Z1R, Z1I ,Z2R, Z2I ,Z3R,Z3I ) 

X1EM=Z2R 

X2EM»Z3R 

CALCULATE THE TUNNELING PROBABILITY AS A FUNCTION OF ENERGY 

IF (X1EM.LT.0.0.OR.X2EM.GT.D) GO TO 10 
CALCULATE TUNNELING INTEGRAL 

Y=RNIHF(RTP0T,X1EM,X2EM,16) 
USE MILLER AND GOOD'S FORM OF WKB APPROXIMATION 

TUNNPB-1/l 1+EXP (Y) ) 

RETURN 

WRITE (6,20) X1EM,X2EM 

TUNNPB=0.5 

RETURN 



C 
C 
20 



FORMAT (IX, 'TURNING POINTS ' ,1X,F8 . 5, IX, ' AND' ,1X, F8. 5, IX, 
2 'ARE OUT OF BOUNDS TUNPB SET = 0.5") 



Listing of RTPOT 



FUNCTION RTPOT (XEM) 
WRITTEN BY E. CLAYTON TEAGUE 
COMMON T,FE,WF,V,D,FK,EM 



APRIL 1978 



POT=PE+WP-EM-(V/D)*XEM-0.2867*D/(XEM*(D-XEM)) 
IF (POT.LT.0.) POT=0. 
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RTPOT=10.25*SQRT(POT) 



RETURN 
END 



Listing of RSIM13 



FUNCTION RSIM13(F,A,B,N) 



C RSIM13 WRITTEN BY CHARLES P. REEVE, STATISTICAL ENGINEERING 

C DIVISION, NATIONAL BUREAU OF STANDARDS 

C FOR: NUMERICAL INTEGRATION BY SIMPSON'S 1/3 RULE 

C SUBPROGRAMS CALLED: USER DEFINED FUNCTION F 

C CURRENT VERSION COMPLETED APRIL 28, 1978 

c - " 

COMMON T,FE,WF,V,D,FK,X 

N2=N/2 

IF(2*N2.NE.N) GO TO 5 

IF(A.GT.B) GO TO 4 

N21-N2-1 

H=(B-AI/N 

H2=2*H 

X=A 

RSIM13=F(X) 

X=B 

RSIM13=RSIM13+F(X) 

X=A-H 

DO 10 1=1, N2 

X=X+H2 

RSIM13=RSIM13+4.*F(X) 
10 CONTINUE 

X=A 

DO 12 1=1, N21 

X=X+H2 

RSIM13=RSIM13+2.*F(X) 
12 CONTINUE 

RSIM13=H*RSIM13/3. 

RETURN 
4 WRITE(6,601) 
601 FORMAT! IX, '*** LOWER LIMIT OF INTEGRATION EXCEEDS UPPER LIMIT* 



Listing of RSIM13 



S4X, 'VALUE OF INTEGRAL SET =0.') 
RSIM13-0. 
RETURN 
5 WRITE (6, 60 2) 
602 FORMAT(lX,'***THE NUMBER OF PARTITIONS MUST BE EVEN**«'/4X, 
S'VALUE OF INTEGRAL SET - 0.') 
RSIM13=0. 
RETURN 
END 



HSJH13 WAS A DUPLICATE COPY OP RSIM13 



Listing of RNIHF 



FUNCTION RNIHF(F,A,B,N) 



RNIHF WRITTEN BY CHARLES P. REEVE, STATISTICAL ENGINEERING 
GET WORK FUNCTION 



C 
C 

C 

10 

C 

c 
c 



WRITE (6,110) 
READ (5,80) WF 

GET MINIMUM ELECTRODE SPACING 

WRITE (6,120) 
READ (5,80) IMIN 

GET MAXIMUM ELECTRODE SPACING 

WRITE (6,130) 
READ (5,80) IMAX 
IDIF=IMAX-IKIN 

GET INCREMENTS OF ELECTRODE SPACING 

WRITE (6,60) 
READ (5,80) NIES 

IF (IDIF.LT.NIES) IINC»1 
IF (IDIF.LT.NIES) GO TO 10 
IINC=IDIF/NIES 

GET MINIMUM VOLTAGE 

WRITE (6,140) 
READ (5,80) JMIN 

GET MAXIMUM VOLTAGE 

WRITE (6,150) 
READ (5,80) JMAX 
JDIF=JMAX-JHIN 

GET INCREMENTS OF VOLTAGE 

WRITE (6,70) 
READ (5,80) NV 



C 
20 



C 

30 



C 

c 

c 
c 

c 

40 
50 



IF (JDIF.LT.NV) JINC=1 
IF (JDIF.LT.NV) GO TO 20 
JINC-JDIF/NV 

WRITE (6,170) T,WF,FE 
WRITE (6,180) 

DO 50 ID=IMIN,IMAX,IINC 
D=ID/1000. 

DO 40 JV=JMIN,JMAX,JINC 
V=JV/1000. 
BLZCTT=8.6E-5*T 
DLTA-40 . *BLZCTT+V 

IF (E-FE) IS LESS THAN 40KT THEN: LN(XPONT) = EV/KT 

ETOP=FE+WF-1.1466/D 

ETOP = TOP OF IMAGE POTENTIAL WITH NO BIAS VOLTAGE. 

SET MAXIMUM TUNNELING BARRIER HEIGHT AS WF + DLTA 

1ST INTEGRAL COVERS ENERGY RANGE WHERE LN(XPONT) MAY BE 
APPROXIMATED BY EV/KT 

A=ETOP-WF-DLTA 

IF (A.LT.0.) A="0 

B=FE-DLTA 

FK=1. 

NTGL1=RSIM13(PR0D,A,B,NP2) 

2ND INTEGRAL COVERS ENERGY RANGE WHERE FULL VARIATION OF 
LN(XPONT) AND TUNNPB MUST BE TAKEN INTO ACCOUNT. NOTE THAT 
SINCE TUNNPB INCREASES AS E INCREASES UPPER LIMIT IS EXTENDED 
TO 80KT 

A=FE-DLTA 

DIFF=ETOP-(FE+80.*BL2CTT) 

B=ETOP 

IF (DIFF.GT.0.1) B=FE+80.*BLZCTT 

FK-2. 

NTGL2=RSIM13(PROD,A,B,NP) 

3RD INTEGRAL COVERS ENERGY RANGE IN WHICH TUNNPB IS ASSUMED 
TO VARY LINEARLY WITH ENERGY, AND WHEN ETOP-FE IS LESS THAN 
80KT. 

IF (DIFF.GT.0.1) NTGL3=0. 

IF (DIFF.GT.0.1) GO TO 30 

A=ETOP 

B=FE+80.*BLZCTT 

FK=3. 

NTGL3-RSIM13(PROD,A,B,NP2) 

RJ=NTGL1 +NTGL2+NTGL3 
WRITE (6,160) D,V,RJ 

AT THE NATIONAL BUREAU OF STANDARDS COMPUTER FACILITY THE FOLLOWING 
STATEMENT WRITES INTO AN AUTOMATICALLY CREATED TEMPORARY FILE CALLED 
7. UPON TERMINATION OP THE RUN FILE 7 MAY BE COPIED INTO THE USERS 
OWN FILE FOR PERMANENT STORAGE. 

WRITE (7,160) D,V,RJ 

CONTINUE 
CONTINUE 



(REAL*)') 
(REALt)' ) 
(REAL*)') 
( INTEGER I ) ' ) 



60 FORMAT (1X,'110W MANY INCREMENTS OF ELECTRODE SPACING?') 

70 FORMAT (lX,'HOW MANY INCREMENTS OF VOLTAGE?') 

80 FORMAT ( ) 

90 FORMAT ( * ENTER TEMPERATURE IN DEGREES KELVIN 

100 FORMAT (' ENTER FERMI ENERGY IN ELECTRON VOLTS 

110 FORMAT (' ENTER WORKFUNCTION IN ELECTRON VOLTS 

120 FORMAT (' ENTER MINIMUM ELECTRODE SPACING IN PM. 

130 FORMAT (' ENTER MAXIMUM ELECTRODE SPACING IN PM. (INTEGER #)' 

140 FORMAT (' ENTER MINIMUM VOLTAGE IN MILLIVOLTS. (INTEGER #)') 

150 FORMAT (' ENTER MAXIMUM VOLTAGE IN MILLIVOLTS (INTEGER!) 1 ) 

160 FORMAT (2F8 . 3 , E12. 5 ) 

170 FORMAT (lHl,/,/,8X, ' TEMPERATURE - MX.FS.l,' DEGREES KELVIN ' / 

2 8X,' WORK FUNCTION = ',1X,F4.1,' ELECTRON VOLTS .' ,/,SX, " " ' 

3 ' FERMI ENERGY = ',1X,F4.1,' ELECTRON VOLTS.',////) 
180 FORMAT (8X,' COLUMN 1 = ELECTRODE SPACING IN NM.',/,8X, 

2 ' COLUMN 2 = BIAS VOLTAGE IN VOLTS. ',/, 8X, 

3 ' COLUMN 3 = CURRENT DENSITY FLOW IN AMPERES PER SQ.M 1 ,////) 



Listing of CURNT 



PROGRAM TUNNEL CURRENT (CURNT) 

WRITTEN BY E. CLAYTON TEAGUE MARCH 1978 

SUBPROGRAMS CALLED: RSJM13, DNSITY 

COMMON T,FE,WF,V 
EXTERNAL DNSITY 
DATA NPART /16/ 

GET TEMPERATURE 

WRITE (6,80) 
READ (5,70) T 



GET FERMI ENERGY 
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WRITE (6,90) 
READ (5,70) FE 

GET WORK FUNCTION 

WRITE (6,100) 
READ (5,70) WF 

GET MINIMUM ELECTRODE SPACING 

WRITE (6,110) 
READ (5,70) IMIN 

GET MAXIMUM ELECTRODE SPACING 

WRITE (6,120) 
READ (5,70) IMAX 
IDIF=IMAX-IMIN 



IF(A.GT.O) GO TO 4 
1I=3.141S92654/SQRT(2.*N)-1/(2.*N) 
X=(B+A>/2 
RNIHP=F(X) 
DO 10 1 = 1, N 
1II = I*H 

Y=((B-A)/2)«TANH(HI) 
D=COSH (til 1**2 

RNIUF=RNIHF+(F(X+Y)+FIX-Y))/D 
10 CONTINUE 

RNIHP=((B-A)/2)*H*RNIIIF 
RETURN 
4 RNIHP=0 

WRITE(6,600) 
600 FORMATl IX, '""LOWER LIMIT OF INTEGRATION EXCEEDS UPPER LIMIT***'/ 
£4X, 'VALUE OF INTEGRAL SET ■ 0. IN SUBPROGRAM RNIHF') 
RETURN 
END 



c 
c 

c 

10 



c 

20 



GET INCREMENTS OF ELECTRODE SPACING 

WRITE (6,50) 
READ (5,70) NIES 

IF (IDIF.LT.NIES) IINC=1 
IF (IDIF.LT.NIES) GO TO 10 
IINC=IDIF/NIES 

GET MINIMUM VOLTAGE 



Listing Of DCUBIC 



WRITE (6,130) 
READ (5,70) JMIN 

GET MAXIMUM VOLTAGE 

WRITE (6,140) 
READ (5,70) JMAX 
JDIF=JMAX-JMIN 

GET INCREMENTS OF VOLTAGE 

WRITE (6,60) 
READ (5,70) NV 

IF (JDIF.LT.NV) JINC=1 
IF (JDIF.LT.NV) GO TO 20 
JINC-JDIF/HV 

WRITE (6,160) T,WF,FE 
WRITE (6,170) 

DO 40 ID-IMIN,IMAX,IINC 
D=ID/1000. 

DO 30 JV=JMIN,JMAX,JINC 
V-JV/1000. 

INTEGRATE THE CURRENT DENSITV OVER SPHERICAL 
SURFACES TO OBTAIN THE NET CURRENT PLOW. 

A-D 

B=D+0.4 

CURNT=RSJM13<DNSITY,A,B,NPART) 

CURHT=CURNT*3 .14159265*1 . 5E-1 2 

WRITE (6,150) D,V,CURNT 



C 

C 
C 

c 

C 

C 
30 
40 
C 

C 

C 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 



AT THE NATIONAL BUREAU OF STANDARDS COMPUTER FACILITY THE FOLLOWING 
STATEMENT WRITES INTO AN AUTOMATICALLY CREATED TEMPORARY FILE CALLED 
7. UPON TERMINATION OF THE RUN FILE 7 MAY BE COPIED INTO THE USERS 
OWN FILE FOR PERMANENT STORAGE. 



WRITE (7,150) D,V,CURNT 



CONTINUE 
CONTINUE 



FORMAT (IX 

FORMAT (IX 

FORMAT ( ) 

FORMAT ( ' 

FORMAT ( ' 

FORMAT ( 



'HOW MANV INCREMENTS OF ELECTRODE SPACING?') 
'HOW MANY INCREMENTS OF VOLTAGE?") 



ENTER TEMPERATURE IN DEGREES KELVIN 
ENTER FERMI ENERGY IN ELECTRON VOLTS 
ENTER WORKFUNCTION IN ELECTRON VOLTS 



FORMAT (' ENTER MINIMUM ELECTRODE SPACING IN PM 



FORMAT ( ' 
FORMAT ( ■ 
FORMAT ( ' 



2 ■ COLUMN 2 

3 ' COLUMN 3 



(REAL!)') 
(REAL!)') 
(REAL!)') 
(INTEGER I) 

( I NTEGER 



SUBROUTINE DCUBIC ( A3,A2, Al ,A0, Z1R.Z1I , Z2R, 221 , 23R, Z3I ) 

DCUBIC WRITTEN BY CHARLES P. REEVE, STATISTICAL ENGINEERING 

DIVISION, NATIONAL BUREAU OF STANDARDS 
FOR: COMPUTING THE ROOTS OF THE CUBIC EQUATION 
A3 * X**3 + A2 * X**2 + Al * X + A0 = 

WHERE THE COEFFICIENTS A3, A2, Al, AO ARE 

DOUBLE PRECISION. 

THE METHOD OF SOLUTION IS TAKEN FROM PAGE 17 OF 

"HANDBOOK OF MATHEMATICAL STATISTICS" BY 

H. ABRAMOWITZ AND I. A. STEGUN, NBS-AMS-55. 
SUBPROGRAMS CALLED: KJD 
CURRENT VERSION COMPLETED MAY 3, 1978 

IMPLICIT DOUBLE PRECISION ( A-Z ) 

IF(DABS(A3).LE.lD-8) GO TO 4 

PI-314159265 35897 93D-15 

A2=A2/A3 

Al=Al/A3 

A0=A0/A3 

Q=Al/3-(A2/3)**2 

R=(Al*A2-3*A0)/6-(A2/3)**3 

T=0**3+R*R 

IFIT.LT.0. ) GO TO 2 

T-DSQRT(T) 

S1C=DCBRT(R+T) 

S2C=DCBRT(R-T) 

GO TO 3 

2 TCR=0D0 
TCI=DSQRT(-T) 
CR-R+TCR 
CI=TCI 

A=DATAN(CI/CR)+PI*KJD(CI,CR) 
A=A/3 

RR=DSQRT(CR*CR+CI*CI) 
RR-DCBRT(RR) 
SlR-RR*DCOSlA) 
S1I=RR*DSIN(A) 

CR=R-TCR 

CI=-TCI 

A-DATAN(CI/CR)+PI*KJD(CI,CR) 

A-A/3 

RR=DSQRT( CR*CR+CI *CI ) 

RR=DCBRT(RR) 

S2R-RR*DCOS(A) 

S2I=RR*DSIN(A) 

3 UCR=-(SlR+S2R)/2-A2/3 
UCI=-(SlI+S2I)/2 
VCR=-(DSQRT(3D0)/2)*(S1I-S2I) 
VCI-(DSQRT(3D0)/2)*(S1R-S2RI 
ZlR=SlR+S2R-A2/3 
Z1I=S1I+S2I 

Z2R=UCR+VCR 
Z2I-UCI+VCI 
Z3R=UCR-VCR 
Z3I=UCI-VCI 
RETURN 

4 WRITE(6,600) 

600 FORMAT(/1X,'NO SOLUTION — COEFFICIENT OF X**3 IS LESS THAN 1E-8 1 
RETURN 
END 



Listing of KJD 



) 



ENTER MAXIMUM ELECTRODE SPACING IN PM. (INTEGER!)') 
ENTER MINIMUM VOLTAGE IN MILLIVOLTS. ( INTEGER! )' ) 
ENTER MAXIMUM VOLTAGE IN MILLIVOLTS (INTEGER!) 1 ! 

FORMAT (2F8.3,E12.5) 

FORMAT (1111,/,/, SX, ' TEMPERATURE » ■,1X,F5.1," DEGREES KELVIN. 1 

2 8X,' WORK FUNCTION = ',1X,F4.1,' ELECTRON VOLTS. ',/ ,BX, 

3 " FERMI ENERGY = , ,1X,F4.1, 1 ELECTRON VOLTS.',////) 
FORMAT <8X,' COLUMN 1 = ELECTRODE SPACING IN NM.',/,8X, 

= BIAS VOLTAGE IN VOLTS .',/ ,8X, 
CURRENT FLOW IN AMPERES ' ,////) 



./, 



FUNCTION KJD(T.B) 



DIVISION, NATIONAL BUREAU OF STANDARDS 
FOR: NUMERICAL INTEGRATION USING THE TANH RULE 

REFERENCE: HABER, SEYMOUR, THE TANH RULE FOR NUMERICAL 
INTEGRATION, SIAM J. NUMER. ANAL., VOL. 14, NO. 4, 
SEPTEMBER 1977 
SUBPROGRAMS CALLED: USER DEFINED FUNCTION F 
CURRENT VERSION COMPLETED APRIL 27, 1978 



KJD WRITTEN BY CHARLES P. REEVE, STATISTICAL ENGINEERING 

DIVISION, NATIONAL BUREAU OP STANDARDS 
FOR: DETERMINING THE QUADRANT FOR THE INVERSE TANGENT FUNCTION. 
THE PROPER MULTIPLE OF PI IS ADDED TO THE PRINCIPLE VALUE 
OF THE INVERSE TANGENT. THE MODE IS DOUBLE PRECISION 
SUBPROGRAMS CALLED: -NONE- 
CURRENT VERSION COMPLETED JUNE 8, 1976 

DOUBLE PRECISION T,B 
IF(T.LT.O) GO TO 1 
IF(B.LT.O) GO TO 2 
KJD=0 
RETURN 

1 IF(B.GE.O) GO TO 3 

2 KJD=1 
RETURN 

3 KJD=2 
RETURN 
END 
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Appendix II. STATISTICAL FLUCTUATIONS IN THE LENGTH OF A BAR 



The fluctuations to be calculated here are those 
arising from the thermodynamic properties of mat- 
ter and do not account for any additional quantum 
fluctuation In the length of the bar. At the tempera- 
tures used for this experiment, the quantum fluctua- 
tions wouid be negligible. The calculations are 
based on the development given by Landau and 
Lifshitz [74]. Starting from the argument that the 
probability for a quantity .v to have a value in the 
range x to x+dx is proportional to e six)/>: , where 
S(x) is the entropy of the system involved and k is 
Boltzmann's constant, they derive the result 



(AfoWrOg) 



(ii-i) 



Here (A Vf is the mean fluctuation of the system 
volume at the temperature T, 

This may be related to the properties of a solid 
by the bulk modulus, B, which is defined by the 
relation 



B 



~KS 



(II-2) 



where V\s the volume of the solid. Combining eqs 
(.11—1 and II-2) gives the result: 



(AVf=kTV/B. 



(II-3) 



If one now assumes that a bar of length, I, and 
cross-sectional area, A , is constrained such that all 



fluctuations are manifest in length changes, then 
the result 



(Mf=(kT/B)(l/A), 



(H-4) 



is obtained, where (A/) 2 is the mean square fluctua- 
tion in the length of the bar. 

Values for the parameters of eq (II-4) which rep- 
resent the gold electrode leads are 



r=300K 

B = I.OX 10 10 dynes/cm 2 

/=1 cm 

A=\0~ 2 cm 2 



Thus, [(A0 2 ] 1/2 =2X 1CT 13 m. Calculated values for 
the fluctuations of other parts of the electrode as- 
sembly, i.e., the aluminum housing and the copper 
holders, are also close to this value for the gold 
electrodes. Since this fluctuations in all the parts of 
the electrode assembly would be uncorrected, the 
net results would be the square root of the sum of 
the squares of the individual components. The re- 
sultant mean square fluctuations in the electrode 
spacing would then be approximately 4x 10~ 15 m. 
Note that these estimates include all frequency 
components. 
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Appendix III. VACUUM CONDITIONS TO ACHIEVE AND MAINTAIN A CLEAN 
SOLID SURFACE 



Let T be defined as the time in which an initially 
clean surface becomes covered with a monolayer 
of adsorbate when the surface is exposed to a gas at 
a pressure P Torr, then the purpose of this approxi- 
mate expression for T is 



r = 2xlO _ V/ > seconds. 



(UI-D 



The following derivation assumes that the surface 
contains 10" adsorption sites/m 2 and that the stick- 
ing coefficient for adsorption is one, i.e., every ad- 
sorbate molecule which strikes the surface is 
adsorbed. 

According to the kinetic theory of gases, the 
number of molecules striking unit area per unit 
time, Z, is given by 



Z=N A P(2TrMR g T)~ 1! \ 



(111-2) 



where AT A is Avagadro's number, M is the molecu- 
lar weight of the gas, R s is the universal gas con- 



stant, and T is the gas temperature. For Nitrogen, 
N 2 , at 300 °K, 

Z = 2. 873 X 1 12 P molecules/(m 2 s). (III-3) 

Converting the units of /> from N/m 2 to Torr gives 

Z = 3.83X 10 14 P molecules/(m 2 s). (III-4) 

Time for one monolayer is density of adsorption 
sites/Z, so 

7 , = (10 1 Vm 2 )/Z = (2.6x 10- 6 )/P seconds. (111-5) 

The general form of eq (III-l) for N adsorption 
sites/unit area is 



T=CVN(2vMR s T) y2 /(N A P), 



(IfI-6) 



where C is a proportionality constant determined 
by the adsorbate's sticking coefficient. 



222 



Journal of Research of the National Bureau of Standards 



Appendix IV. AN ELECTRONIC DEVICE UTILIZING QUANTUM MECHANCAL 
TUNNELING 



The idea for this electronic device resulted from 
attempts to design a metal-vacuum-metal tunneling 
structure with the use of microfabrication tech- 
niques commonly employed by the microelectron- 
ics industry. Conceptually, the device would 
combine the properties of piezoelectric ceramics 
with the high sensitivity of tunneling to changes in 
electrode spacing to produce an electronic device 
with a very high transconductance. 

Starting with a highly polished piezoelectric sub- 
strate which is polarized as a transverse expander, 
i.e., fields applied normal to substrate plane result 
in dimensional changes in the substrate plane, the 
device would be fabricated as shown in figure 38. 
After applying electrodes to both faces of the 
piezoelectric substrate, a thick layer of polycrys- 
talline metal or semiconductor, shown as the dot- 
ted layer, would be deposited over the anticipated 
device area. The next step would be to heavily oxi- 
dize the metal or semiconductor layer or by other 
means form an integral insulating layer, shown as 
the lightly cross-hatched strips. Following this op- 
eration, the tunneling electrodes would be formed 
on top of the insulating layer using standard evapo- 
ration and photomasking technology. The final 
step in the fabrication process would be to remove 
the insulating layer using a preferential plasma or 
chemical etching process to leave a structure simi- 
lar to that shown in figure 38. A gap within the 



range of the piezoelectric substrate would be 
formed by either highly precise photolithography 
or by melting a narrow bridge into two parts. 

Assuming that one hits produced a gap between 
the two electrodes which is shorter than the trans- 
verse dimensional capabilities of the piezoelectric 
material, the device characteristics should have the 
following qualitative properties. Let the piezoelec- 
tric electrodes be A and B and the tunneling elec- 
trodes be C and D. Then, if x } is the lateral distance 
between the bases of the tunneling electrodes and 
<f 3 „ is the lateral expansion coefficient of the 
piezoelectric ceramic, 



(for 1 = (rf 3[ Jc,/-yf / 'A 



(IV-1) 



where t is the ceramic substrate's thickness, and 
dV A B is the change in voltage across the substrate. 
For a fabricated gap of 500 nm, the voltage Vl B to 
place the electrodes within tunneling distances for 
available piezoelectric ceramics would be approxi- 
mately 50 V for a 0. 1 mm thickness. 

Considering the device input to be the voltage 
applied to the piezoelectric ceramic and the output 
as the tunneling current J C d for a Rxed bias voltage, 
the transconductance is defined as 



g = <dlcv/d Pah) Vcd =* const . 



(IV-2) 



\ 




Luy?r 



F-l.-zcelGc-rlc Substrate 
-I-r Transverse PvLzjrlz- 



Figure 38-Proc-ess to fabricate tunneling devices using pho- 
tolithographic techniques. 



Combining eq (IV-1) and (IV-2) with the rule-of- 
thumb tunneling equation 



/ CD =C 1 exp(C 1 S2J 1/2 )^F'cD, 



(IV-3) 



where C, = 3.16Xl0 13 A/m 2 /V, C 2 = 10.25 (ev)- 1/2 
nm" 1 , S is the electrode spacing, B the mean barrier 
height and ,4 is the tunneling area, gives 



g=C 1 C ! 8 lfl e qa " fl ^^^ 



(IV-4) 
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For tunneling areas of 10 
*, = 4 mm, r =0.1 mm, B 
F CD =0.1 V. 

g = \Q- 1 A/V{ = mho). 



i 2 , </ 3 i=0.25 nm/V, 
ev, S = l nm, and 

(IV-5) 



This value is about the maximum transconductance 
obtainable with conventional electronic devices 
[75,76]. Increases in tunneling area to 10~ 10 m 2 
would raise g to 1 mho. 

A photograph of a device resulting from the ex- 
ploration of this technique for fabricating tunneling 
structures is shown in figure 39. The approach was 
not pursued any further, because at the time of the 
investigation, only limited time was available and 
the preferential etching process to produce the re- 
quired electrode length of 2 mm was not suffi- 
ciently developed. Shorter lengths place demands 
on photomasking technology that could not be 
achieved, i.e., forming lines less than one microme- 
ter in width. However, since the time of these at- 
tempted fabrications, April 1976, the capabilities of 
both technologies have significantly advanced. 




Figure 39-SEM photomicrograph of prototype device fabri- 
cated with the use of photolithographic techniques. 
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Appendix V. ELECTRONICS FOR THE TEMPERATURE CONTROL SYSTEM 



The purpose of this appendix is to give a descrip- 
tion of the two-stage temperature control system 
and its associated electronics in terms of the trans- 
fer functions of each part of the system, to derive 
from these transfer functions the response of each 
stage's control system to setpoint changes and envi- 
ronmental temperature changes, and finally to de- 
scribe the procedure employed to choose the lead 
and lag times of the compensation network and to 
"tune" the network for an acceptable response to 
setpoint changes. Descriptions of both stage-con- 
trollers are given in terms of the block diagram 
shown in figure 23. The intent here is not to give a 
detailed analysis of the electronic circuitry, but 
rather to give a reasonably complete description of 
the basic properties of the control system in terms 
of the properties of its subunits. The transfer func- 
tion of each subunit will now be obtained. 

The transfer function of the thermistor with re- 
spect to ambient changes in temperature, T, can be 
obtained from the equation for conservation of en- 
ergy 



Q E = C^f+K(T z -T), 



(V-l) 



7* z = thermistor temperature, Qe= electrical power 
input, C= thermal capacity of the thermistor, and 
K= dissipation constant. From this equation, the 
effects of input power level may also be obtained. 
At a steady state temperature, r| s , T z s — T REF 
=Q E /K. With the input power levels and dissipa- 
tion constants described below, 7f s — r outer =0.1 
mK and 7f s — r inner = 1 ju,K. The Laplace transform 
of eq (V-l) gives the result 



T z (S)=[l/(l+ZS)]T(S) 

+Z/[C(\+ZS)]Q E (S), 



(V-2) 



Z=K/C is the time constant of the thermistor, and 
S is the Laplace transform variable in the transform 
domain. The thermistor transfer function with re- 
spect to ambient changes in temperature is there- 
fore 



The time constant, Z, and dissipation constant, K, 
for both controllers are strongly dependent on the 
manner in which the thermistor is mounted. The 
estimates given below for K were based on the 
manufacturer's recommendations and the thermal 
conductivities of the mounting materials. Z must 
account for the thermistor's response time and the 
coupling of the thermistor to the heater and system 
mass. For the thermistor's response, independent of 
this last factor, Z was approximately 1 to 10 sec- 
onds. Accounting for the imperfect coupling of the 
thermistor and heater, rough estimates for Z based 
on observations of system behavior are 

Z outer =slO-20s, 
and 
Z imKr -100-150s. 

The coupling of the thermistor and heating coils on 
the inner stage was across 6 mm of the inner alu- 
minum enclosure. 

Let Z a , Z b , Z c , and Z d represent the four 
impedances in a Wheatstone bridge. Then, if 
Z=Z a =Z h =Z 4 and Z c =Z+dZ, with dZ much 
less than Z, the error signal, e, output by the bridge 
is [77] 



e=(E/4)(dZ/Z), 



(V-4) 



E is the bridge excitation voltage. Thermistors 
have the property that Z=Z(T Z ), where T z is the 
device temperature. Over a quite large working 
range Z(T Z ) for the thermistors used in the con- 
troller bridges was given by Z(T Z )—Z(T KBF ) 
=a(r z -r REF ) with a = -0.04 Z(T KEF )/K. This 
along with the thermal transfer function of eq (V- 
3) gives 



e=A(S)(T REF -T z ), 



(V-5) 



Af(5)= T Z {S)/T(S)= 1/(1 +ZS). 



(V-3) 



where A(S)=k/(\+ZS), A:=0.01E V/K. 

To minimize the temperature difference between 
the thermistor and its surroundings, the self-heating 
of the thermistors by the bridge excitation voltage 
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must be kept to a minimum [78]. The manufactur- 
ers recommendation was that the self-heating be 
less than or equal to 0. 1 of the desired control reso- 
lution. The dissipation constants of the two ther- 
mistors as mounted on the two thermal enclosures 
were 90 mW/K and 42 mW/K for the outer and 
inner enclosures, respectively. At bridge balance, 
one-half of the exciting voltage will be across the 
thermistor. Thus, the powers which may be dissi- 
pated in the thermistors assuming 1 mK stability in 
the outer stage and 10 fiK stability in the inner 
stage are 9.0 ju,W and 42 nW, respectively. The rec- 
ommended excitation voltages are, in turn 0.85 V 
and 0.13 V for the outer and inner stages respec- 
tively. Since any variations in these voltages would 
result in an additional variation of the control sig- 
nal, stable sources are required. As a compromise 
between the recommended excitation voltage and 
the need for stability, a 1.35 V mercury cell was 
used for the outer-stage bridge. The inner-state 
bridge was driven by a stabilized oscillator with an 
output of 100 mV rms. 

The circuit diagrams for both stages are shown 
in figures 40 and 41. The lock-in amplifier em- 
ployed for the sensor amplifier was an instrument 
designed and built in-house at the National Bureau 
of Standards [79]. It had gain settings from 500 up 
to 10 6 , time constants available from 0.5 s to 20 s, 
and operated at a fixed frequency of 1 kilohertz. To 
a first approximation, the transfer functions of both 
sensor amplifiers may, for small deviation from 
bridge balance, be represented by the equation 



H,(S)= V,(S)/e(S)=A/{\ +qS), 



(V-6) 



A is the amplifier's dc gain or for the lock-in ampli- 
fier its gain when the rms value of the input is not 
changing, and q is the time constant of the respec- 
tive amplifier. This transfer function for the lock-in 
amplifier represents the relationship between the 
dc output voltage and the rms value of the input ac 
signal with the selected phase. For the lock-in am- 
plifier, q is simply the time constant setting of the 
amplifier's low pass filter. Equation (V-6) assumes 
the lock-in filter is equivalent to one RC stage. For 
the sensor amplifier of the outer-stage-controller, q 
was set by the 0.004 u-fd capacitor at 4xl0~~ 2 s 
since this was the approximate value expected from 
the 12 dB/octave gain characteristic and antici- 
pated gain setting. 

For both stages, the bridge-sensor amplifier com- 
bination was configured to obtain the required neg- 
ative feedback. With the outer stage controller, this 
only amounted to choosing the correct polarity of 
the bridge excitation voltage. With the lock-in am- 
plifier, this was accomplished by choosing the ap- 
propriate phase to detect. 

A brief note about the particular bridge circuit 
used for the outer sensor amplifier. This circuit al- 
lows the use of a single-ended, chopper-stabilized 
amplifier and its attendant low drift. Another ad- 
vantage of the circuit is that the gain is not depen- 
dent on the absolute value of the bridge resistors 
and can be simply adjusted by changing the ratio of 
the 10 nfl and 100X1 resistors. 



SUPPLY 



Figure 40-Schematic of electron- 
ics for outer temperature con- 
trol stage. 
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SUPPLY 



Figure 41-Schematic of electron- 
ics for inner temperature con- 
trol stage. 
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The transfer function of the compensation net- 
work shown in figure 42 is given by [80] 



(1+fflSKl + fcS) 



(V-7) 



where P is the proportional gain setting, / is the 
product of the integration channel's gain and the 
inverse of the integration time and D is the product 



of the differentiation time constant and this chan- 
nel's gain. The time constant, a, was 0.2 s for the 
outer-stage controller and 2.0 s for the inner-stage 
controller; b was 0.5 D for the outer-stage and 0.05 
D for the inner-stage. The open loop transfer func- 
tion of the operational amplifiers used in the com- 
pensation network is to first order 



M(S)=2X KT7(1 + 6X 10' 3 5). 
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Figure 42-Schematic of compen- 
sation network used for both 
temperature control stages. 
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Since this pole and the pole at <S = 100 in the pro- 
portional channel are so much higher frequency 
poles than all the other poles of the controllers, 
their effects, which are minor, will not be included 
in the following discussion. 

The response properties of the power amplifier 
for both control stages can be understood by not- 
ing that a fundamental property of an operational 
amplifier is to maintain the summing point, i.e., the 
inverting input, at a "virtual ground" condition 
[80]. The implication of this property for the power 
amplifier portion of the controllers is that all cur- 
rent flow to the inverting input of the operational 
amplifier which drives the MOSFET power tran- 
sistor, a VMP-12 or 22, must also flow through the 
10 Kfl feedback resistor, /? FB . connected to the 
device's source terminal. (This, of course, neglects 
the small amount flowing through the 100 Kfl re- 
sistor to the operational amplifier output.) This 
condition can be used to relate the voltage applied 
to the heater, V3, and the equivalent voltage out- 
put from the compensation network with an equiv- 
alent input resistance, R it of 100 Kfls. Thus, with 
these definitions; V2=I R„ where / is the summed 
current. The voltage at the source terminal will 
then be maintained at; y FB =f R FE . But, K FB =(T3/ 
RkWb, from a simple application of Ohm's law. So 



V3=(R fb /RWh/Rb) Y2- 



(V-9) 



Here R H is the heater resistance and R B is the bias 
resistance between common and the device source 
terminal. 

The transfer function of the power amplifier is 
determined by the dc gain obtained from eq (V-9) 
and the response of the operational amplifier in the 
feedback loop. The corner frequency of the MOS 
device is greater than 100 MHz. Thus, the overall 
transfer function of the power amplifier would be 



H$5)=V3/V2=A l> /(l+pS), 



(V-10) 



with Ap=(R P v/R-X r h/ r b) and/> = 6xl0~ 3 s. As 
indicated for the compensation network, the effects 
of this high frequency pole are not accounted for in 
the following discussion. 

To obtain a model of the thermal response of the 
heaters and their respective masses, a simple 
lumped parameter model of the heat conduction 
and radiation was utilized. Referring back to figure 
22, let R represent the thermal resistance of all con- 
duction paths between each stage's components 
and its environment. For the outer stage, these 
paths include the polyurethane foam, the stainless 
steel ports to the vacuum pumps, the electrical con- 



necting leads, and the epoxy supports used to ther- 
mally isolate the outer stage from the granite 
vibration isolation table. For the inner stage, these 
paths include the electrical connecting leads and 
the ceramic supports. By analogy with electrical 
circuit terminology, R is given by R =(R&—T^)/ 
Qci where T$ = temperature of stage, T B = temper- 
ature of the stage's environment and Qc = power 
loss due to conduction. 

Radiation transfer of energy is a major factor 
which must be considered in the design of the inner 
control stage. Zemansky [81] shows that the rate of 
energy transfer Q R , between a body with emissiv- 
ity, e, and temperature, T, and a surrounding sys- 
tem at temperature, T B , is given by 



Q K =A€<r(T 4 -T E t ). 



(V-ll) 



Here .4 is the surface area of the enclosed body and 
or is the Stefan-Boltzmann constant. 

With the use of this lumped parameter model, 
the thermal behavior of each stage can be de- 
scribed by the equation for the conservation of en- 
ergy 



Q=Qm-Qc~Q* 



(V-12) 



or 



C~ = V 2 /r - (T-T,)/R -A €cr(T*- T E 4 ), 

(V-13) 

where C is the heat capacity of each stage, V is the 
voltage applied to heater coils of resistance r, and 
the other symbols are as defined previously. If eq 
(V-13) is linearized by expanding T and V about 
their operating point values, Tref and F , at which 
d7Yd/ = 0, the following equation holds for small 
deviations about these values 



1 



6T 1 
dt^R'C ~RC ' R'C 



r=^F+, 



(V-14) 



where l/R' = l/R +4AecrT REP \ Let t T =R'C and 
t E = rC, then the Laplace transform of eq (V-14) 
gives the result 



nS) = G ] (S)V(S)+G 1 (S)T E (S) 



(V-15) 



where Gi(S)=2Fto-/(f,(l-KpS)) and G 2 (5)=l/ 
(1 + tf). 

The transfer function of the cascaded sensor am- 
plifier compensation network and power amplifier 
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is given by the product of their individual transfer 
functions [82]. Therefore, let 

#(s)= V3/e =HiQS)XXdS)XHdS). (V-16) 

The controllers' feedback responses are then deter- 
mined by the equations 

e^AiSXTxzp-T) 

V=V 3 =H(S)e 

T=G ) (S)V+G 2 (S)T E . (V-17) 

The transform variable is not indicated on the 
dependent variables. Solving this set of equations 
for T gives 

T=[GIS)T R +A {S)G,{S)H (S)T REF \ F($l 

F(S)=[\+A(S)G l (S)H(S)]-\ (V-18) 

The following discussion is presented in terms of 
changes m 7* REF . While the transient response of the 
system is different for changes in 7r EF and T E , the 
stability conditions are the same for both system 
changes. Calculations were only made in terms of 
changes in 7* REF since step changes in T^? were 
anticipated. 



Let K(S)=T(S)/T KW (S), then 
K(S)=G l (S)H(S)A (S)F(S). 



(V-19) 



K(S) is the overall transfer function describing the 
response of the closed loop system to changes in 
the setpoint, T REP . Substitution of the various 
parameter and function definitions yields the result 

G X (S)H(S)A(S) 

g(S+S l )(S+S 2 )(S+S 2 ) 

S(l + f„SXl+«S)(l+6SXl+ZS)(l+ffS), 

(V-20) 

where g =2kAA p V ( jt T /t e . 

The zeros S lt S 2 , and 5 3 are the roots of the equa- 
tion 

S' + tf/p + l/b + l/a +D/(abP)]S 2 

+ [I/bP + I/aP + 1 /ab ]S +I/abP = 0. (V-2 1) 

This equation reduces to the following equations 
for the two stages; 



outer stage: 

S 3 + 15S 2 + [51+ W/D]iS + 10 1/D =0, (V-22) 

and inner stage: 

S>+ 1 \S 2 + [20/+ 10/D]jS + 10 1/D =0. (V-23) 

For eqs (V-22 and V-23), the roots are approxi- 
mately [83]: 

outer stage: 

S liBi -15 > 5 1 ,S 3 «-(2S,)- 1 [] 1 ±Sf l []i 1/ ', 

and inner stage: 

S l3 - 1 1, S lt S^-{2S t )- ] [ h±Sf i [ V n . 

Parameters of the outer-stage controller were 

heat capacity= C = 1 .42 X 10* J/"C 
heater resistance = r = 60 fi 
thermal resistance^ =0.8 °C/W 
steady state voltage = V =20 V 
bridge constant = fc=0.01 V/°C 
stage time constant =R C = t T = 1 . 1 X 10 4 s 
power amplifier gain =A P = 6 

The poles, other than the high frequency ones not 
accounted for in eq (V-20), were a = 0.25s, 
6=0.5Ds, Z = 10-20s, and gr=4xI0 -2 s. 
Parameters of the inner-stage controller were 

heat capacity =C = 2.9 XlO 3 J/°C 
heater resistance = r = 1.05X 10 3 £2 
thermal resistance =R =8 °C/W 
steady state voltage = T = 10 V 
bridge constant — k = 1 x 10~ 5 V/°C 
stage time constant=i?C=2.3X 10 4 s 
power amplifier gam=A P = 3 

The poles were located at a=2.0s, 6=0.05Ds, 
Zsl00-150s, and ?sQ.5-2s, The thermal resis- 
tance R of this stage represents the parallel resis- 
tance of conduction and radiation. 

The poles and zeros of K(S) in the 5-plane com- 
pletely determine the stability of the temperature- 
control stages and their transient responses to 
changes in setpoint temperature and environmental 
temperature [82]. Some guidelines for the design of 
a feedback system which are usually scattered 
throughout most texts on control systems are 

A. The system will be unstable if K (S) contains 
any poles in the right half-plane. 

B. The system's transient response is increas- 
ingly oscillatory as the poles of K(S) move 
toward the imaginary axis. 
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C. The system's speed of response is propor- 
tional to the distance of the dominant poles 
from the S -plane origin. 

D. Introduction of negative real zeros in K(S) in 
general increases system stability. 

E. The relative weighting, R ir of the k !h mode of 
a system, i.e., a pole or a complex pair of 
poles, is determined by 



\RA=- 



n \P,-n,\ 

i— ] 



ia i rr ip.— p. i 



(V-24) 



where n h I— I, m are the zeros of K(S) and 
Pj,j—l, n, are the poles of K (S). Thus, plac- 
ing zeros near troublesome or oscillatory 
poles minimizes their effect. 
The response of a feedback system (as in fig. 
23) to a constant reference input has a zero 
steady-state error if and only if the open loop 
transfer function G\(S)H(S)A(S) has at least 
one pole at S =0. For zero steady-state error 
with reference inputs of the form, at", the 
requirement is that Gi{S}H(S)A (S) have a 
pole at S = with multiplicity m greater than 
n. 

For disturbance inputs, a similar property 
holds. Here the pole at S = must be in the 
controller transfer function H (S). Therefore, 
a controELer with a pole at S = insures that a 
zero steady-state error may be achieved in 
both cases. 



An attempt was made to use these guidelines in 
adjusting the lead and lag times of the compensa- 
tion networks. The procedure for adjustment was 
to set I and D to zero and P = 1 . The overall gain of 
the feedback loop was then increased by increasing 
A , the sensor amplifier gain, until the stage-con- 
troller broke into oscillation. D and / were then 
increased to the values calculated below; D to 
eliminate the instability by introducing zeros near 
the origin and / to eliminate or minimize any offset 
by introducing a pole at the origin. The effects of 
these changes can be visualized in terms of a root- 
locus diagram. This type of diagram is a graph of 
the paths traced by the poles of K(S) as g of eq 
(V-20) is varied from zero to infinity. A root-locus 
diagram, constructed with the aid of the root-locus 
properties described by Fortmann and Hitz [E4] 
and Kuo [85], for K(S) of the outer-stage con- 
troller is sketched in figure 43. Only the dominant 
poles and zeros are drawn since the effects of the 
other ones are principally to change the angles of 
the asymptotes of the root-loci. 

An examination of the diagram in figure 43 will 
show that before compensation, the two dominant 
poles are the ones due to the thermal lag of the 
stage mass and to the time constant of transferring 
energy from the heater to the thermistor. These are 
labelled \/t T and 1/Z, respectively, in the diagram. 
The compensation network introduces three zeros 
and three new poles. The zeros tend to straighten 
the root-locus branch that causes instability so that 
the same gain may be used with no oscillations. An 
equal margin of stability without compensation 
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would have required significantly reduced loop 
gain and, thus, less attenuation of environmental 
disturbances. 

D and I were initially chosen to bring the zeros, 
S 2 and 5 j, to be adjacent to the pole at 1/Z to ap- 
proximately cancel its effect. From the roots given 
for eqs (V-22 and V-23), this condition requires 
that for the outer stage 



30 



[5/+10/£>]=0.1 



and for the inner stage 



22 



[20I + 10/D] = 0m 



l=0AS-\D = 3s 
and for inner stage: 
1=0.01 S~\D = 50s. 

After some "tweaking" of I, D, and P to get 
maximum gain without oscillations, the responses 
to setpoint changes shown in figures 44 and 45 
were obtained. These responses are certainly not 
optimum, but are acceptable. The sensor amplifier 
gains and compensation network settings to obtain 
these responses were: outer stage — A = 5000 V/V, 
7 = 0.013 r 1 , Z> = 4 s, 7^ = 1, and inner stage — 
A =5X 10 5 V/V, 7=0.015 s l , D = 200 s, i>=0,75. 

With these sensor amplifier gains, the closed 
loop gain, including the bridge and thermistor re- 
sponse, were 



Initial setting for D and I were 
for outer stage: 



outer stage: A CL =A A p k= 3,000 V/'C 
inner stage: A CL =A A p k = 15,000 V/"C. 
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